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Uber stabilisierte Netzgeraite mit Transistoren 


E- Von Ernst BALDINGER und WoLFGANG Czaja, Basel*) 


Einleitung 


Die Zahl der Anwendungen stabilisierter Netzgerate als Gleichspannungs- 

quellen niederen Innenwiderstandes ist durch die Méglichkeit, Transistor- 
: Netzgerate zu bauen, wesentlich vergréssert worden. 

Mit verhaltnismassig kleinem Aufwand kann man heute Probleme lésen, 
die bisher die Verwendung von Akkumulatoren notwendig machten. Beriick- 
sichtigt man, dass sich die Ausgangsspannung der stabilisierten Netzgerate sehr 
einfach iiber grosse Bereiche variieren lasst, so erkennt man, dass Transistor- 
Netzgerate in vielen Anwendungen Akkumulatoren iiberlegen sind. Ferner gibt 

_ der Temperaturkoeffizient der Ausgangsspannung solcher Netzgerate zu keinen 
Schwierigkeiten Anlass. 

Die hier behandelten stabilisierten Netzgerate sind ausserdem zur Sonate 
komplizierterer Transistorschaltungen, Stabilisierung der Heizung von Rohren- 
schaltungen usw. interessant und diirften die Lésung vieler Probleme auch bei 
industriellen Geraten erleichtern. 

Im folgenden wird gezeigt, dass man die Figenschaften dieser Netzgerate 
beherrscht und mit den heute kauflichen Transistoren bereits sehr brauchbare 
Resultate erreichen kann. 

Fiir den Transistor beniitzen wir durchwegs das Ersatzschema nach Figur 1, 
das in der Literatur zum Teil als natiirliches Ersatzschema bezeichnet wird — 
[1-4]2). 

€ 


Figur 1 


Ersatzschema des Transistors. 


1) Physikalisches Institut der Universitat Basel, Abteilung fiir angewandte Physik. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 25. 
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2 ERNST BALDINGER und WOLFGANG CZAJA — 
_ Die Verwendung des Ersatzschemas nach Figur 1 bietet drei wesentliche 
Vorteile. 

a) Die Behandlung von Transistorschaltungen wird der von Réhren- 
schaltungen ahnlich. Durch passende Grenziibergange erhalt man Verhiltnisse, 
- wie sie bei der Verwendung von Réhren vorliegen. 

b) Die Arbeitspunktabhangigkeit (I,- bzw. I,-Abhangigkeit) der Parameter 
lasst sich in guter Naherung durch einfache Funktionen darstellen [2—4]. Dies 
_ gestattet zum Beispiel, die Lastabhangigkeit des Stabilisierungsfaktors bei 
Stabilisierungsschaltungen zu berechnen. 

c) Die Aufstellung eines derartigen Ersatzschemas lasst sich physikalisch 
begriinden [1, 2]. Die Parameter haben eine einfache physikalische Bedeutung. 

Typische Werte der Parameter eines Niederfrequenztransistors kleiner 
Leistung OC 71 sind zum Beispiel R, = 3002, R, proportional 1/J,, wobei 
R, = 1,3 kQ bei I, = 1 mA. Die differentielle Stromverstarkung SR, des Tran- 
sistors in Emitterschaltung bei kurzgeschlossenem Ausgang ist nahezu konstant, 
mit dem Mittelwert SR, ~ 55%). Weiter erweisen sich R, und R, ebenfalls in 
guter Naherung als proportional zu 1/Z, und nehmen bei J, = 1 mA die Werte 
R, + 2,5 MQ R, » 40 kQ an. 

Den Widerstand der Emitterzuleitung werden wir in allen Rechnungen als 
geniigend klein voraussetzen und von seinem Einfluss absehen. Schliesslich 
k6énnen wir die Kapazitaten C, und C; so lange vernachlassigen, als wir uns mit 
gentigend langsamen Vorgangen befassen. Zur Diskussion von Stabilitatsfragen 
sind hingegen diese Kapazitaten wesentlich. Von allgemeinem Interesse ist die 
Feststellung, dass die bei etwa 1kHz gemessenen Parameter auch bei sehr 
langsamen Veranderungen (« Gleichstromschwankungen ») ihre Werte praktisch 
beibehalten. 


1. Einfache Stabilisierungsschaltungen 


1.1. Stabilisierungsschaltungen 


Stabilisierte Netzgerdte stellen im allgemeinen als Schaltungen Sechspole 
dar (Figur 2). 

Die Ausgangsspannung ist eine Funktion der Eingangs- und der Referenz- 
spannung. 


Us = U4 (Uy, Up) . 


3) Aus Messungen an einer Vielzahl verschiedener Transistortypen ergab sich, dass bei Nieder- 
leistungstransistoren (bis einige hundert Milliwatt) SR, verglichen mit R, nur schwach von J, ab- 
hangig ist und im Hinblick auf die Berechnung von Schaltungen in vielen Fallen in ausreichender 
Naherung als konstant angenommen werden kann, Bei Leistungstransistoren scheint hingegen all- 
gemein die I,-Abhangigkeit von SR, starker zu sein. (Die Arbeitspunktabhangigkeit yon SR 
entspricht im wesentlichen der Websterschen Theorie iiber hohe Injektionsdichten. Vel. nin 
Beispiel [20].) 
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Figur 2 


Stabilisierungsschaltung als Sechspol. 
u,, Eingangsspannung; wu, Ausgangsspannung; #, Referenzspannung; I; Laststrom. 


Wir fragen nach dem Einfluss von AUg und AUp auf uw, und entwickeln u, 
dazu in eine Taylor-Reihe*) 


= Ou, Ous 
Be Malavp-0 eS on pee ee IE ek as ee AUp rs 
sete AE R=? AUR=0 (1) 
(“4)aug-0 = Us (Ug, Up) - 
AUR=0 


Eine derartige Entwicklung ist immer dann méglich, wenn wir geniigend kleine 
Abweichungen AU, und AUp von einem Bezugszustand betrachten. Unsere 
Schaltung soll in dem Sinne als linear zu betrachten sein, als das Abbrechen 


der Entwicklung (1) mit den linearen Gliedern auch bei endlichen AUg und 


AUpz bereits eine gute Naherung darstellt. Sowohl U, als auch die Koeffizienten 
von AU,g und AUg in (1) werden im allgemeinen vom Laststrom, das heisst 
vom Arbeitspunkt der Transistoren abhangen. 

Weitergehende Aussagen lassen sich erst dann machen, wenn wir die All- 
gemeinheit einschranken und an die Schaltung zusatzliche Anforderungen 
stellen. Dazu definieren wir den Stabilisierungsfaktor 7 als Verhaltnis 


$8 AU 1, 
int anata (2) 
und den Innenwiderstand 
= Ous 
Se eis ieee y (3) 


Beziiglich der Ausgangsklemmen lassen sich bei bekannter Ausgangsspannung U4 
die stabilisierten Netzgerite eindeutig durch 7 und R; beschreiben, das heisst 
als Spannungsquelle #4, mit dem Innenwiderstand R, darstellen. Dabei gilt 


1 AU 
et aU, =U, (be - eae 


yj und R,; hangen vom Arbeitspunkt, das heisst von der Belastung ab. 


4) Wir bezeichnen eine konstante Spannung im folgenden mit U,, die tiberlagerte Rippelspan- 
nung bzw. Schwankung mit 4U,, so dass u, = U, + AU, gilt. 


’ - - 
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Die Stabilisierungsschaltung soll eine méglichst gute Approximation an 
eine ideale Spannungsquelle darstellen, das heisst eine méglichst konstante 
Gleichspannung liefern. Dies verlangt 

1. einen méglichst kleinen Innenwiderstand R, derart, dass J; R;/U,<1, 

2. einen moglichst grossen Stabilisierungsfaktor 7, so dass 4U4/AUp<1 
und AU,4/U,<1 gilt. 

3. Werden wir verlangen, dass sich 7 und vor allem FR; mit der Frequenz 
nicht stark andern, das heisst fiir schnelle Schwankungen der Eingangsspan- 
nung und des Laststromes ahnliche Werte aufweisen wie fiir extrem langsame. 

Die zweite Bedingung verbietet eine starke Abhangigkeit der Ausgangs- 
- spannung von wz. Es ist daher plausibel 


Uy = Uy(Up) + Duy, = <1 (4) 


zu schreiben. Diese Uberlegung gilt auch dann, wenn wir unsere Schaltung aus 
nichtlinearen Elementen aufbauen®). - ; 

Verwendet man keine Bauteile mit nichtlinearen Kennlinien, so ist die Ein- 
fihrung einer Referenzspannung und gesteuerter Spannungs- bzw. Strom- 
quellen (Réhren, Transistoren usw.) notwendig. Die Einfiihrung gesteuerter 
Quellen erlaubt, dass die abgegebene stabilisierte Leistung durch den Eingang 
geliefert und die Referenzspannung nur schwach belastet wird. 

Eine Stabilisierungsschaltung wird also entweder passend nichtlineare Ele- 
mente oder eine Referenzspannung®) und gesteuerte Quellen enthalten. Im 
folgenden werden wir die zweite Variante vorziehen, die sich mit normalen 
Bauelementen realisieren lasst und weitaus anpassungsfahiger ist. 

Es ist nicht erforderlich, dass die Referenzspannung eine extreme Konstanz 
aufweist. Man kann sich zum Beispiel eine ebenfalls vom Netz gespiesene, 
eigene Hilfsspannung als Referenzelement mit der Eigenschaft 


tpn = Up + eo ae AU are (5) 


denken. Setzen wir (5) zusammen mit (4) in (1) ein, so erhalten wir 


ot i oR 1 [( 04a Ou, OUR 
U4 (U4 (UR) } b Ur) if (ae4 | satay a lean) apes lose area | 
AUR=0 AUR=0 | (1a) 
xX AUZ +++, AIDE Ale 
A 


°) Derartige Elemente lassen sich stets durch eine Spannungsquelle und eine Serieimpedanz 
darstellen. 


6 eh West ; : ; . 
) oder ein als Referenzelement wirkendes, wesentlich nichtlineares Element. 
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EDurch Briickenschaltungen’) kann der eesumte Koeffizient von AU, in (1a) 
_oder ein Teil, zum Beispiel 


(Fre) 
OUR AUp=-0’ 


_zam Verschwinden gebracht werden. Da die Differentialquotienten, wie wir 

bereits betonten, ebenfalls Funktionen von J; sind, wird eine derartige Kom- 
pensation vor allem fiir Netzgerate von Vorteil sein, deren Belastung konstant 

ist. Praktische Kompensationsschaltungen werden aus linearen Elementen auf- 
gebaut, sie kompensieren also nur die linearen Terme in (5) bzw. (1a); héhere 
Terme sind stérend und werden nicht kompensiert. Die Tatsache, dass bei 
praktischen Geradten eine ausgezeichnete Kompensation auch fiir gréssere AU; 
méglich ist, kann als Beweis fiir die Zulassigkeit der Naherungen (1) und (la) 
gewertet werden. 

Wir wollen diese allgemeinen Betrachtungen nun verlassen und uns der 
Behandlung der Stabilisierungsschaltungen in Figur 3 zuwenden, wobei wir 
unseren Uberlegungen zunidchst die Schaltungen in der einfachsten Form 
zugrunde legen. 

a bezeichnen wir als Schaltung mit Quertransistor, b als Schaltung mit 
Langstransistor. Beide Varianten sind gebrauchliche, von den Réhrenschaltun- 


Figur 3 


Einfaché Stabilisierungsschaltungen: @ mit Quertransistor; b mit Langstransistor. 


gen her bekannte Anordnungen [6-8], die im Zusammenhang mit Transistoren 
bereits erwihnt wurden [9-11, 13, 22]. 


1.2. Die Schaltung mit Quertransistor 


Der Einfachheit halber nehmen wir zundchst die Referenzspannung als 
ideal an. Es soll also wz = Ug und R; zr = 0 gelten. Ausserdem vernachlassigen 
wir den Innenwiderstand des speisenden Netzgerates. Spater werden wir uns 
von diesen einschrankenden Voraussetzungen wieder befreien. 


7) Unter Briickenschaltung verstehen wir hier allgemeine Schaltungen im Sinne von Bope [5], 
S. 243. 
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Aus der Schaltung nach Figur 3a ergibt sich sofort, dass 
lg Up 3 Uns (6) 


wobei in den meisten praktischen Fallen die Basis-Emitter-Spannung U,, klein, 


das heisst U,,/Up<1 ist. Zur Berechnung von 7 bestimmen wir zunachst das — 


Verhaltnis AU;/AU,, das wir aus der Ersatzschaltung (Figur 4) entnehmen 
kénnen: AUg/AU, = (R+ R,)/R; Wegen R,/R<1 lasst' sich dies zu 


dU ===; 
AU; R A # 
dle 
Figur 4 Figur 5 
Ersatzschaltung zu Figur 3a. Zur Berechnung des Innenwiderstandes. 


AU;/AU, » R/R; vereinfachen. Ferner gilt Ug — Uy =I R8), wobei mit I 
die Summe aus Querstrom Jp, das heisst dem Strom durch den Quertransistor 
in Figur 3a und Laststrom I, bezeichnet wird. Im Leerlauf ist I = Ig. Mit der 
Definition von (2) erhalten wir 


Us 


nme (1-2), no Be a 


Aus diesen Beziehungen liest man die fiir die Dimensionierung von stabilisierten 
Netzgeradten geeigneten Zusammenhange 
U4 
Noo Ri = TH 7 Re= te (7a) 
ab. 7. bedeutet den optimalen Stabilisierungsfaktor, den man fiir Uz > oo 
erreicht. Aus (7) entnimmt man, dass eine Erhéhung der Eingangsspannung 
tiber den Wert von etwa 2 Uy nicht empfehlenswert ist. Eine grossere Ein- 
gangsspannung bedeutet nur eine bescheidene Verbesserung, die sich asymp- 
totisch dem Faktor 2 nahert, wihrend die Eingangsleistung unbeschrankt 
wachst. Weiter zeigt sich, dass das Produkt aus Innenwiderstand und Stabili- 
sierungsfaktor durch die Gleichstromgréssen allein festgelegt ist. Beachtet man 
ferner, dass U, und I von der Last weitgehend unabhangig sind, so ist sofort 
ersichtlich, dass die Bezichungen (7) und (7a) ihre Giiltigkeit auch bei Belastung 
des Netzgerates beibehalten, wenn nur die entsprechenden Werte fiir 7 und 
R, eingesetzt werden. 
Wir fragen nun nach dem Innenwiderstand R,. Da stets dI, = —adlo¢ gilt, 
erhalt man fiir den differentiellen Innenwiderstand aus (3) 


aE esc 8) 


8) Man beachte, dass dies eine Beziehung zwischen « Gleichstromgréssen » darstellt. 


‘ 


- 
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Die fiir die Berechnung des Innenwiderstandes massgebende Ersatzschaltung 
ist in Figur 5 dargestellt. Mit Hilfe des eingangs angegebenen Ersatzschemas 


fiir den Transistor erhalten wir im Grenzfall sehr langsamer Stérungen 


oe ee R> Ri. (9) 


Beschranken wir uns nicht nur auf langsame Schwankungen, so enthalt (9) 
-statt R, die entsprechende Impedanz (vgl. Figur 1). Besitzt ferner die Refe- 
renzspannung einen Innenwiderstand ungleich Null, so ist in (9) als R, die 
Summe aus Basiszuleitungswiderstand des Transistors und Innenwiderstand 
der Referenzspannung einzusetzen. 
Es bleibt nun noch die Diskussion der Lastabhangigkeit von R; und 7; 
wegen (7a) kénnen wir uns auf die von R; beschranken. Mit (9) und den in der 
Einleitung angegebenen Abhingigkeiten der Parameter vom Arbeitspunkt 


ergibt sich : R, LJt 


0 ail Ae See Tee iLO) 


R; ra (R;) 


Dieser Ausdruck vereinfacht sich, wenn der Querstrom gross gegen den Last- 
strom ist (bzw. R, sehr klein ist) zu 


Riv (Ry, -0 p= pr 1 Mrz-0 (1 — LID. (11) 


(11) hat den Vorteil, dass die Eigenschaften desTransistors nicht bekannt sein 
miissen. Sie erlaubt auch fiir unbekannte Transistoren eine gute Abschatzung, 
die etwas zu ungiinstige Werte liefert. R, und 7 sind in Figur 6 als Funktion der 


06 08 40 
— lf ~~ —- if > 
a b 
Figur 6 
Schaltung mit Quertransistor nach Figur 3a. 4 Lastabhaingigkeit von R,; 6 Lastabhangigkeit von 
Noo- Gleichung (11); —— — Gleichung (10) bzw. Gleichung (10) und (7a). Messpunkte: 


Schaltung mit OC 71; im Leerlauf erhalt man die Werte oo = 35,8 berechnet, 7o9 = 33,8 gemessen, 
R, = 10Q berechnet, Rj = 11Q gemessen. ; 


| 
| 
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Figur 7 : 
Schaltung mit Langstransistor nach Figur 3b. a Lastabhangigkeit von R,; b Lastabhangigkeit von 
Nopt: Gleichung (16) und (17); ——— Gleichung (17a). Fiir eine Schaltung mit einem 
OC 71 ergibt die Rechnung (Ri)7, +00 fm 5,5 Q, (Mopt)r, =0 © 1690. Beste Ubereinstimmung der 


= 6,02 und (%op4)7, _9 = 1700. 


Messwerte mit der Rechnung erhalt man fiir (R,) Tp->00 Tr 
Belastung aufgetragen. Baut man mit dem eingangs erwahnten Transistor eine 
Querschaltung auf, so stimmen die gemessenen Werte innerhalb der Mess- 
genauigkeit mit der Rechnung iiberein®). Es ist jedoch zu beachten, dass die 
hier berechneten funktionellen Abhangigkeiten von der Belastung nur unter 
der Voraussetzung SR, ~ const gelten (vgl. S. 2). Die gemessene Abhangigkeit | 
von 1/SR, ist von der Form A + BI, und hat die Tendenz, die Lastabhangigkeit 
von R; und 7 abzuschwachen, mindestens solange SR, im Vergleich zu Ry nur 
schwach von J, abhangt. 

Wir haben bisher den Innenwiderstand des speisenden Netzteiles null 
gesetzt. Man tiberlegt sich nun leicht, dass bei Schaltungen mit Quertransistor 
der Innenwiderstand R; y des speisenden Netzteiles in Serie zum Vorwiderstand 
R liegt und deshalb einfach beriicksichtigt werden kann. Im allgemeinen ist in 
dieser Schaltung der Einfluss von R;y jedoch klein. 

Bis auf die Berechnung des Innenwiderstandes und die Belastungsabhangig- 
keiten besitzen die Uberlegungen dieses Abschnittes insofern eine allgemeine 
Bedeutung, als sie auch auf die Verhiltnisse bei Stabilisatoren mit Glimmréhren 
und Zener-Dioden [11, 12] direkt anwendbar sind. 


1.3. Die Schaltung mit Léngstransistor 


Um zu einer einfachen Behandlung der Langsschaltung zu gelangen, be- 
trachten wir zunichst den Einfluss des Seriewiderstandes R; y des speisenden 


xt) Langsame Schwankungen, die der Transistor ausregelt, bewirken eine Arbeitspunktver- 
schiebung und damit eine Temperaturanderung des Transistors. Dieser zusatzliche Effekt, der hier 
nicht berticksichtigt wurde, lasst sich anhand von Abschnitt 3.1 angeben, ist jedoch nicht von 


grosser Bedeutung. Im erwihnten Beispiel ist die Korrektur stark lastabhangig und variiert von 
~ 3% bis ~ 10%. 


TEN rere eee aT 


#t 


Vol. IX, 1958 Uber stabilisierte Netzgerate mit Transistoren 9 


Netzteiles auf den Innenwiderstand und den Stabilisierungsfaktor. Eine strenge 
Behandlung der Schaltung nach Figur 3b mit dem Ersatzschema des Transi- 


_stors nach Figur 1 zeigt, dass, solange 


Ry, 
R@tsSR) 
Bs R, + By 


Rin < (12) 
ik 


oh 


_ gilt, der Innenwiderstand bei kurzgeschlossenem Eingang berechnet, das heisst 
_ Ry = 0 gesetzt werden kann. 


- sammenhang 


Der Einfluss von R;y auf den Stabilisierungsfaktor lasst sich wie folgt er- 
fassen: Zwischen den Schwankungen 4U# am Eingang?®) und den Spannungs- 
schwankungen AU; zwischen Kollektor und Basis in Figur 36 besteht der Zu- 
Zz 


mit der Eingangsimpedanz Z (vom Netzteil aus gesehen) der Schaltung in 


_ Figur 3d. Wir berechnen nun 


iN 
AUR Us Riy\ AUg Us Up pil: rane 
ele ee =o . =H : 
oo Au. Us SS eg Oe ei 1 LR 
Ug 


Da die gleichstrommissige Eingangsimpedanz U;/I <Z sein muss, wird 


We fom 
wenn die Bedingung 
ee (13) 
aesar <— 


erfiillt ist. Es bleibt nun noch die Bestimmung von Z. Dazu hat man zunachst 
zu beachten, dass in der Schaltung nach Figur 30 bei einer Anderung AJ des 
Eingangsstromes die entsprechende Anderung des Laststromes um den Faktor f 
grésser ist. Die Differenz fliesst durch die Basisleitung des Transistors. Der 
Wert von f liegt sehr nahe an eins. AU, kénnen wir nun einerseits aus dem 
Spannungsabfall iiber dem Lastwiderstand und andererseits aus dem Stabili- 
sierungsfaktor bestimmen. Setzen wir beide Ausdriicke gleich, so ergibt sich 


AU, U; U, S 
= z BR Gl - da Bal. 


Die Bedingung (13) ist die bei weitem einschrankendere, das heisst, der Einfluss 
von R;v ist auf den Stabilisierungsfaktor grésser als auf den Innenwiderstand. 


10) Ux bzw. AUZ = EMK bow. Anderung der EM K, Ug bzw. AUpz = Klemmenspannung 
bzw. Anderung der Klemmenspannung des Gleichrichters mit dem Innenwiderstand k; y. 


10 Ernst BALDINGER und WoLFGANG CzajJA ZAMP 


Die Ungleichungen (12) und (13) diirften wohl fiir viele praktisch auftre- — 


tende Falle erfiillt sein; wir werden daher im folgenden vom Innenwiderstand 
R;y des Netzteiles absehen, also R;y = 0 setzen. 
Die Ausgangsspannung der Schaltung mit Langstransistor nach Figur 3b 


bestimmt sich zu 
U, =e 0. (14) 


Es ldsst sich leicht einsehen, dass der Innenwiderstand R, dieser Schaltung 
durch (9) gegeben ist. Zur Bestimmung des Verhaltnisses AU;/AU, stellen 
wir unsere Schaltung formal durch eine Ersatzschaltung entsprechend zu 
Figur 4 dar. Die beiden Widerstaénde legen wir als Parameter dieser Ersatz- 
schaltung durch folgende Bedingungen fest: 
1. Vom Ausgang aus gesehen soll die Ersatzschaltung den gleichen Innen- 
widerstand wie Figur 30 besitzen. 
2. AU;/AU, soll fiir beide Schaltungen den gleichen Wert aufweisen. 
Aus 1. folgt, dass der zum Ausgang parallel liegende Widerstand gleich R; 
wird, wahrend man mit 2. den zweiten Widerstand R, streng zu 
R 
Ba ure R 
Rutsey t+ Rar, (+54 Re) 
R,+R, (1+ SR,) 4 Ry 
R, R; 


t+ 


berechnet. 

Mit diesem Vorgehen erreicht man, dass die Behandlung der Langsschaltung 
formal auf die gleiche einfache Art wie die der Querschaltung erfolgen kann. 
Man hat dabei jedoch die Lastabhangigkeit von R, zu beachten. Es gilt 


AUsz i? Ra + Ry Ra Rk; 
40 et Ra eh 


Damit wird der Stabilisierungsfaktor 


ST ge lant Te (15) 


Optimale Verhiltnisse liegen offenbar dann vor, wenn U,/Ug = 1. Man wird 
also bei der Schaltung mit Langstransistor bestrebt sein, U; méglichst klein 
zu wahlen, das heisst nur die minimal notwendige Kollektor-Emitter-Spannung 
an den Transistor anzulegen. Aus (15) erhalt man sofort 


Noot Ri = Ry, (15a) 


und dieses Produkt ist im Gegensatz zur Schaltung mit Quertransistor, vom 
Arbeitspunkt, das heisst von der Belastung des Netzgerates abhiangig. 


_ 
: 


— 
+3 
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i Wir untersuchen nun noch die Lastabhangigkeit von R; und y. Die Varia- 
_ tion von Rk; mit der Belastung ergibt sich leicht zu 

E 1 

q FoR) ie (1 a =) , (16) 
£ 


wobei wir die dimensionslose Grésse x = R,/R, = R, I;/Ry") eingefiihrt 


_ haben und (Rj);, 0 den Wert R;/(1+ SR,) besitzt. Der Stabilisierungsfaktor 
_ variiert fiir kleine x, das heisst kleine Belastungen, naherungsweise wie 


R, 1 
opts 2p (Novt)r, =0 Tan (17) 


7 


mit 
R,SR 


(Nope)t, +0 > x A val 


In (17) wurde dabei R, = (R,)7, +0 > R, eingesetzt. Fiir beliebige x erhalt man 
mit den immer zulassigen Vereinfachungen R,/R,<1, R,/R; <1, SR, > 1, 
R 


mS a — 
Hat 1 RaSRs % 
ee eta y 
und weiter 
1 
Nopt = (Nopi)r, -0 Pk eR SR (17a) 
Td ed 
Rgo 


Nach (16) geht R; > co fiir J, > 0. Betreibt man die betrachtete Schaltung 
im Leerlauf, bzw. bei stark wechselnder Belastung, so ist es vorteilhaft, eine 
geringe zusatzliche Belastung einzubauen, die R,; passend begrenzt. Weiter ist 
bemerkenswert, dass 7,; fiir A, > 0 lastunabhangig wird und (R,)7,.. > 0 
geht. Figur 7 gibt eine Darstellung der Arbeitspunktabhangigkeit von 7 und 
R,. Beziiglich des Einflusses der J,-Abhangigkeit von SR, und des Temperatur- 
effektes gilt das Entsprechende wie auf Seite 8. 


1.4. Erginzende Bemerkungen zum Verhalten der Schaltungen nach Figur 3 


a) Varianten der Schaltung der Referenzspannung. Fiir manche Anwendungen 
kann es zweckmassig sein, die in Figur 3 angegebene Schaltung wie folgt abzu- 
andern (Figur 8): 

Mit y bezeichnen wir das Spannungsteilerverhiltnis, dessen Definition aus 
Figur 8 entnommen werden kann. Wir bemerken, dass der Spannungsteiler 
eine zusatzliche konstante Belastung des Netzgerates darstellt. Der Span- 
nungsteiler soll ferner durch den Transistor praktisch nicht belastet werden. 


11) Fiir Parameter R, ~ 1/I, schreiben wit R, = Rjo/l, (A= 2,3, 4): 
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Zur Definition von ¥ 


Uy 


YU, 


om Y= 


Figur 8 


Stabilisierungsschaltungen mit Spannungsteiler. 


Eine allfallige Belastung des Spannungsteilers durch den Transistor kann durch 
den Faktor y beriicksichtigt werden. Bei den Schaltungen mit Quertransistor 
erhalten wir jetzt fiir die Ausgangsspannung 


Up + Uz 
U, = - j 
und es ist anschaulich verstandlich, dass der Stabilisierungsfaktor 7 um den 
Faktor y verkleinert, der Innenwiderstand um 1/y vergréssert wird, wobei das 
Produkt beider Gréssen erhalten bleibt, vorausgesetzt, dass die Gleichstrom- 
daten des Gerates nicht verandert werden. Bei den Schaltungen mit Langs- 
transistor erhalt man nach Figur 9 eine Ausgangsspannung von 


Ug+U,, 
Ley 


Us 2 

und sowohl fiir den Innenwiderstand als auch fiir den Stabilisierungsfaktor 
gilt das bereits oben erwahnte Verhalten, wobei y durch (1 — y) zu ersetzen ist. 
Besitzt ferner die in Figur 8 eingezeichnete Spannungsquelle u, den Innen- 
widerstand R;,p, so tritt dieser in beiden Schaltungen additiv zum Basiszulei- 
tungswiderstand hinzu. 

Verwendet man als Referenzelemente Zener-Dioden [11, 12], so erhalt man 
zum Beispiel die in Figur 9 dargestellten Schaltungen. 

Zener-Dioden kénnen naherungsweise als Serieschaltung einer Spannungs- 
quelle (Elektromotorische Gegenkraft) mit einem Widerstand dargestellt wer- 
den, womit die Behandlung der Schaltung analog zu dem soeben gegebenen 
Beispiel erfolgen kann. 


Figur 9 


Stabilisierungsschaltungen mit Zener-Dioden als Referenzelement. 


(CAPA T Ea 


a Te: 
va 
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b) Wirkungsgrad der Stabilisierungsschaltungen. Wir fragen nach dem Wir- 
kungsgrad & der beiden Varianten von Figur 3. € soll wie tiblich das Verhiltnis 
von abgegebener (stabilisierter) Leistung zu Eingangsleistung bedeuten. Fiir 
die Schaltung mit Quertransistor erhalt man 


ey Cae il Pe 
Unt e Pare 


Wie bereits erwahnt, wird U,/Uz; ~ 1/2, und ein Verhltnis von Ip/I ~ 20% 


__ ist erfahrungsgemass noch verniinftig. Damit wird im glinstigsten Fall é ~ 40%. 


In allen anderen Fallen wird ein entsprechend grésserer Bruchteil der hinein- 
gesteckten Leistung im Quertransistor vernichtet. Es kann von Vorteil sein, 
durch einen passend gewahlten Widerstand in der Kollektorleitung den Transi- 
stor vor Uberlastung zu schiitzen. &, R; und 7 werden jedoch von dieser Mass- 
nahme nicht beeinflusst. Die Schaltung mit Langstransistor ist demgegentiber 
wesentlich giinstiger, hier ist ndmlich 
b= Catt, Us te 
Cate U1, * 


und dieses Verhaltnis liegt unabhangig von der Belastung wohl in den meisten 
Fallen nahe an eins. 

Nimmt man bei der Langsschaltung in Kauf, dass die notwendige Vorbela- 
stung den Wirkungsgrad vor allem bei kleiner Last verschlechtert, so fallt der 
Vergleich eindeutig zugunsten der Langsschaltung aus. Dies um so mehr, wenn 
man beachtet, dass unter Verwendung gleicher Transistoren die Schaltung 
nach Figur 36 einen um Grdssenordnungen héheren Stabilisierungsfaktor auf- 
weist als die Schaltung nach Figur 3a, wahrend die Innenwiderstande in der 
gleichen Gréssenordnung liegen. 

Als Vorteil der Schaltung mit Quertransistor ist vor allem ihre Kurzschluss- 
sicherheit zu erwahnen, wahrend die maximale Ausgangsspannung durch die 
zulassige Emitter-Kollektor-Spannung begrenzt wird. Eine derartige Be- 
schrankung existiert bei den Langsschaltungen nicht. Die héchste zulassige 
Ausgangsspannung wird hier durch die Grdésse der zulassigen Schwankungen 
der Eingangsspannung bedingt. Unter Umstanden wird man allerdings zusatz- 
liche Schaltungsmassnahmen ergreifen miissen, um den Langstransistor im 
Falle eines Kurzschlusses, bzw. bei zu starker Belastung, zu schiitzen (vgl. 
zum Beispiel [13)). 

c) Serieschaltung stabilisierter und unstabilisterter Quellen. Der Spannungs- 
bereich von Transistor-Netzgeraten lasst sich erh6hen, wenn man zum Ausgang 
der Schaltungen in Figur 3 eine unstabilisierte Quelle in Serie schaltet (vgl. 
Figur 10). Diese Massnahme erlaubt, die gegeniiber den Anodenstromen der 
iiblichen Réhren grossen Kollektorstréme der Transistoren auch bei Netzge- 


raten mit hédheren Spannungen auszunttzen. 


" s ealiteeliel sell 
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Figur 10 
Serieschaltung von stabilisierten und unstabilisierten Quellen. 
Upo + Upi = Up Upg + U4i = Ua, Spannungsteilerverhaltnis y = Ral(Ry + Rp). 
Voraussetzungen tiber den Spannungsteiler vgl. Abschnitt 1.4a. 


Wir untersuchen zunachst das Verhalten der gemass Figur 10a veranderten 
Querschaltung von Figur 3a und sehen vorerst vom Innenwiderstand Rj, der 
unstabilisierten Quelle ab. Nehmen wir Ugg > Ug, an, so kann man AUg, 
gegen AUgy vernachlassigen, das heisst AUg ~ AUgg setzen. Man erkennt 
folgendes Verhalten der Schaltung: A Ugo verursacht Variationen AJ = AUgo/R 
von I. Diese Stromanderungen werden nun vom Transistor nahezu ausgeregelt, 
es bleibt eine Schwankung AU, = AI R; tibrig, wobei mit Ri der Innenwider- 
stand der Querschaltung mit dem Spannungsteiler nach Figur 10a aber mit 
Uso = 0 und Ro = 0 bezeichnet wird. Man erhalt somit 


AUs = AUson =e 


AUP AU ee 


und 


/ 
Nrotal — 1 = : 
ola. Uro 


n’ ist entsprechend zu Rj definiert. Diesen Uberlegungen liegen die normaler- 
weise erfiillbaren Voraussetzungen R’/R <1, Rj/R, <1 zugrunde. Der 
Innenwiderstand der Anordnung nach Figur 10a wird mit R;, =.0 zu 


eS , 
‘total to 


Der Einfluss eines R,) + 0 auf den Stabilisierungsfaktor lasst sich nun 
sofort angeben. R,) wird nur vom Laststrom durchflossen, wir haben also 
anstelle von AUgy jetzt 


AUX, ae AUz os Rig Bhs 


j 
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mu setzen. Es ergibt sich schliesslich 


2 Nica = Neotal R Pet 
- ; a0 
ve 


_ wobel 775,,, mit R;g < R, in das vorher abgeleitete Notar Ubergeht. Mit ent- 
_ sprechenden Uberlegungen erhalt man den Innenwiderstand zu 


R; 
Baal v4 Bes al (1 eee F 
Wir wollen zur Berechnung des Stabilisierungsfaktors noch zwei Bemer- 
_ kungen anfiigen: 
1. Der Einfluss der Schwankungen 4U;, kann mit Hilfe des Superposi- 
tionsprinzips leicht berticksichtigt werden. Man erhalt zum Beispiel 


ke’. 
AU, = (AUg_ + AU; 1) » usw. 


2. Da mit Ugg > Ug, die absolute Grésse der Schwankungen A Ug betracht- 
lich werden kann, ist die Rechnung unter Umstanden durch die Beriicksich- 
tigung der Arbeitspunktabhangigkeit von Rj zu verfeinern. Man erhalt dann 


UA, T, 
y AU, =7fat=/R dl, 
Va, ty 


und dieses Integral lasst sich mit (11) bzw. (10) streng lésen. 

Wird zum Ausgang eines stabilisierten Netzgerates in Langsschaltung 
(Figur 3b) eine unstabilisierte Quelle in Serie geschaltet (Figur 10d), so ergibt 
sich als neue Schaltung wieder ein normales Netzgerat in Langsschaltung,. 
Kann der Innenwiderstand der unstabilisierten Quelle vernachlassigt werden, 
so entnimmt man y und R, aus Abschnitt 1.4a. Der Einfluss eines von Null 
verschiedenen Widerstandes FR; lasst sich entsprechend zum Abschnitt 1.3 
beriicksichtigen. Man erhalt 


Be 
Niotal = 7 i Ty ee 
ifmeteres es 
U; 


und in einer hdheren Naherung als der in 1.3 mE ANS: 


R? * aR, Sg 
Yotal oc opt 


mit Rip < Ry, Rip < Ry (1+ S Rz). Rj, n/ usw. bedeuten dabei die Gréssen 
der Schaltung unter Beriicksichtigung des Spannungsteilers, aber mit Ry) = 0. 


TE 
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Voraussichtlich ist es in vielen Fallen vorteilhaft, in derartige «aufgestockte » 
Netzgerate einen Regeltransformator mit Servomotor-Antrieb einzubauen. Die— 
unstabilisierte Spannung lasst sich dann jeweils so nachstellen, dass die Stabi- 
lisierungsschaltung stets im giinstigsten Regelbereich arbeitet. Ferner verbes- 
sert man mit dieser Massnahme die Eigenschaften der Schaltung fiir geniigend 
langsame Schwankungen erheblich, wahrend fiir kurzzeitige Eingangsstorungen 
die soeben abgeleiteten Beziehungen giiltig bleiben. 


2. Methoden zur Verbesserung der Stabilisierungsschaltung 
2.1. Einbau zusdatzlicher V orverstarker 


Es erscheint selbstverstandlich, dass mit dem Einbau von Vorverstarkern 
in die bisher diskutierten einfachen Stabilisierungsschaltungen ihre Eigen- 
schaften wesentlich verbessert werden kénnen. Derartige Gleichstromverstar- 
ker sind zwischen der Basis des stabilisierenden Transistors und der Referenz- 
spannung anzubringen. Die Auswirkungen dieser Massnahme kann man bei der 
Querschaltung sofort tiberblicken. Bezeichnen wir mit « > 1 die Verstarkung des 
Vorverstarkers, so bleiben unsere Rechnungen unverdndert, wenn wir dem 
stabilisierenden Transistor nun die Steilheit S’ = w S zuordnen. Es ergeben 
sich somit folgende Zusammenhange: 

Kom =, epithelia 

Treten im Transistor selbst wesentliche Riickkopplungseffekte auf, wie dies 
zum Beispiel in der Langsschaltung tiber R, und R, der Fall ist, so versagt das 
soeben bentitzte einfache Verfahren. Eine genauere Analyse der Schaltung 
zeigt, dass unter den Voraussetzungen 


R, R, RR, 
ae Sb. Sly \Raeea ep os pepe ee ae 


ebenfalls 


Ri = = , er =a Nopt) 
gilt. R, wird, wie zu erwarten war, nur sehr schwach von u beeinflusst, und zwar 
im Sinne einer Verminderung der Lastabhangigkeit. Die Rechnung zeigt, dass 
im Ausdruck fiir R, (vgl. Abschnitt 1.3) statt R; nun w R, auftritt. 

Ks soll noch bemerkt werden, dass die Schaltungen nun kompliziertere, zum 
Teil mehrfach riickgekoppelte Systeme sind und es notwendig ist, die Stabilitats- 
kriterien sorgfaltig zu untersuchen. Eine theoretische Behandlung derartiger 
Systeme findet man in [5] und [8]. 


) Die mit ’ gekennzeichneten Gréssen beziehen sich auf eine Schaltung mit der Verstarkung 


4 > 1, die tibrigen auf die gleiche Schaltung, aber mit esis 


3 
2 
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2.2. Kompensation des Einflusses von Schwankungen der Eingangsspannung — 


Wir haben bereits im Abschnitt 1.1 auf die Méglichkeit einer Kompensation 
_ des Einflusses der Schwankungen der Eingangsspannung auf die Ausgangs- 
'spannung hingewiesen und untersuchen nun eine derartige Kompensations- 
_schaltung an einem Netzgerat mit Quertransistor nach Figur 3a, jedoch im 
- Leerlauf. Die Schaltung ist dann wie folgt zu variieren: 


Figur 11 


Kompensationsschaltung fiir ein Netzgerat nach Figur 3a. 


_ Dabei bedeuten R; p den Innenwiderstand der Batterie bzw. der Referenzspan- 
nung, R, einen Seriewiderstand in der Basiszuleitung des Transistors, R, den 
Kompensationswiderstand, und R, ist durch (9) gegeben. Man berechnet AU, zu 


Rat FR, 1 
oe se ORS Re le 
TE ee ALE R,+Rt+Rrpt kh esR, | ; x) 
A E ie Reo : (n+ = 
re Bp a L+Sik, 
Fir 
Rir thi ( SR; 7) 1 
R; 1+ SR R 
etn ss et 


erhalt man eine abgeglichene Briickenschaltung mit AU, = 0. Man beachte, 
dass R;, > 0 sein muss und somit ein Abgleich der Briicke fiir R, die Giiltigkeit 
der Ungleichung ES eee 
R; 
bedingt. a 
Der Innenwiderstand R; der Schaltung nach Figur 11 wird durch R, nur 


schwach beeinflusst: 


oe Re 
R; =a (R;) pp +00 ( 5 Ry . Ry ) R; > 1. 
Rirt ky Rk, +k, 


Dabei bedeutet 
= Rint te R : 
(Rap ued = Ri a mis Ay” >i. 


ZAMP IX/2 


ae 
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Im Falle der exakten Kompensation erhalt man 


D D Riz + Fy 
R; a (Ri) g,r00 ( Hi : (= ire 1) ) 
a Re 


z 


das heisst praktisch den Wert fiir R, > oo. 

Der Stabilisierungsfaktor der Schaltungen in Figur 3 kann also mit zwei 
Methoden vergréssert werden: ae 

1. Durch den Einbau zusatzlicher Vorverstarker (vgl. Abschnitt 2.1). 

2. Durch Kompensation (vgl. Figur 11); dabei sind die Innenwiderstande 
der so modifizierten Schaltungen natiirlich verschieden. Wir vergleichen die 
beiden Varianten 1 und 2 der Querschaltung nach Figur 3a, fordern in beiden 
Fallen das gleiche 7 und erhalten mit den Voraussetzungen 


(R,) (Ri)p 


Rp-oo ko 
ee a RS 


. 


eine Aussage iiber die Lastabhangigkeit des Stabilisierungsfaktors: 


(arr) R 
paperless 
Oly Ju >1 


k 

Die Lastabhangigkeit von 7 ist demnach bei einer kompensierten Querschal- 
tung ungefahr um den Faktor yw grésser als die Lastabhangigkeit des gleichen 
Stabilisierungsfaktors einer entsprechenden Querschaltung mitVorverstarker. 

Es sei noch bemerkt, dass der Einfluss von Lastschwankungen auf die Aus- 
gangsspannung ebenfalls kompensiert werden kann, wobei mit dieser Methode 
auch negative Innenwiderstande realisierbar sind. Eine interessante Anordnung 
ist ktirlzich von PEMBERTON [14] angegeben worden. Vergleiche auch die Arbeit 
von MIDDLEBROOK [22]. 


3. Temperaturkoeffizient und Beispiele ausgefiihrter Geridte 
3.1. Temperaturkoeffizient der Ausgangsspannung 


Man tiberzeugt sich zunachst anhand der Uberlegungen des Abschnitts 1.4b 
leicht davon, dass die Ubertragung von Schwankungen der Referenzspannung 
auf die Ausgangsspannung bei Schaltungen mit Quertransistor im Verhaltnis 
1:y bei Schaltungen mit Langstransistor im Verhaltnis 1:(1 — y) erfolgt. Be- 
schranken wir uns auf die einfachen Schaltungen in Figur 3, so gilt mit y =k 
bzw. y = 0 


AU An 
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Be eteactitet man weiter nur die Schaltung mit Quertransistor (Figur 3a), so 
_ erhailt man aus (6) die Aussage 


Wy Aves dey 
a ar aT (18) 


- BES lasst sich zeigen, dass die Anderungen des Innenwiderstandes und des Quer- 
_stroms mit der Temperatur nur einen von hdherer Ordnung kleinen Einfluss 
auf U, haben. 

; Fir dU,,/dT kénnen wir schreiben 


7 reat (245) a9 


_wobei U.; die Emitter-Junction-Spannung (@ U,, — R, 4, da R, praktisch 
unabhangig von J, und T) bedeutet. Eine nahere Begriindung fiir (19) geben 
_ wir im Anhang. Mit (19) erhalt man aus (18) fiir langsame Anderungen 


Se e+ BYR), (<0. eo 


Aus (20) liest man die Méglichkeit der Kompensation des Temperaturganges 
der Ausgangsspannung ab. Der Temperaturkoeffizient der Referenzspannung 
(zum Beispiel einer Zener-Diode) kann den des Transistors mindestens teilweise 
aufheben. 

Der gemessene Temperaturkoeffizient der Ausgangsspannung einer Schal- 
tung nach Figur 3a (OC 71), bei dUp/dT = 0, ergab sich zu etwa — 3,5 mV/°C 
bis zu 40°C. Oberhalb 40°C erhalt man einen etwa dreimal grésseren Wert. 
Nach (20) und Figur 17 berechnet man den Temperaturkoeffizienten in befrie- 
digender Ubereinstimmung zu — 3,0 mV/°C bis zu etwa 45°C. 

Die Uberlegungen lassen sich ebenfalls auf Schaltungen mit Langstransi- 
storen anwenden, wobei aber nun zu beachten ist, dass wegen (14) der Tempe- 
raturkoeffizient der Referenzspannung mit demselben Vorzeichen eingeht wie 
der Temperaturkoeffizient des Langstransistors, das heisst, eine Kompensation 
des Temperatureinflusses der Ausgangsspannung entsprechend (20) nicht ein- 
tritt. Abgesehen von diesem Unterschied bleibt die Rechnung dieselbe.’ 

Weist die zu untersuchende Schaltung einen Vorverstarker auf, so ist die 
erste Stufe wesentlich fiir den Temperaturkoeffizient der Ausgangsspannung 
und bestimmt zusammen mit dem Temperaturkoeffizient der Referenzspan- 
nung den der Ausgangsspannung. 

Der Temperaturkoeffizient des Verstarkers lasst sich nun durch geeignete 
Schaltungsmassnahmen (symmetrische Stufen) weitgehend beherrschen (vgl. 
Abschnitt 3.2a und [15]), und durch passende Unsymmetrien die Temperatur- 
koeffizienten anderer Bauelemente (Referenzspannung usw.) kompensieren. 
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Eine genauere Analyse der Verhiltnisse ist auf Grund der in den Figuren 
16 und 17 dargestellten Zusammenhange prinzipiell méglich, wiirde hier jedoch — 


zu weit ftihren. 


3.2. Beispiele ausgefiihrter stabilisierter Netzgerate mit Transistoren 


a) Schaltung mit Quertransistor (Figur 12). Dieses Gerat, das bereits kurz 
beschrieben wurde (vgl. BALDINGER und CzaJA[17]), ist als Laborgerat ent- 


7 (60W) 


@) Intermetall 28 
@ Philips 0C 76 
@® Philips OC 16 
@ Delco 2N277 
©) _termetall FD7 


Figur 12 


Stabilisierungsschaltung mit Quertransistor. Der Schalter S reduziert bei dauernder geringer 
Belastung den Strom durch den Quertransistor zur Vermeidung unnG6tiger Warmeentwicklung. 


wickelt worden und weist in einer neueren Ausfiihrung eine Ausgangsspannung 
von 1+ 17 V auf. Der maximale Ausgangsstrom betragt 2,3 A bei 1 V und 
0,6 A bei 17 V. Der Stabilisierungsfaktor ist rund 1000 und wird im wesentlichen 
durch den Stabilisierungsfaktor der Referenzspannung bestimmt. Die Referenz- 
spannung selber ist ein kompensiertes Netzgerat mit Quertransistor. Wie die 
Erfahrung beim Bau von drei gleichartigen Schaltungen lehrte, ist der zwei- 
stufige symmetrische Verstarker geniigend unempfindlich auf die Toleranzen der 
Elemente, und ‘seine Temperaturabhingigkeit gibt zu keinen Schwierigkeiten 
Anlass. Als Leistungstransistor wird ein DELCO 2N277 verwendet. 

Der Einbau der Diode FD7 erweist sich als notwendig, um unter allen Um- 
standen eine Ubersteuerung der Endstufe des Verstirkers und damit Kipp- 
erscheinungen zu verhindern. Die eingebauten RC-Glieder dienen in der in [5] 
und [8] beschriebenen Art zur Korrektur des Frequenzganges und verhindern 
damit eine Instabilitat der Schaltung. 


at 
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Der Temperaturkoeffizient der Anseanesspannuie betragt rund + 0,5mV/°C. 
In Figur 13 ist der Drift der Ausgangsspannung in der ersten Stunde nach dem 
BPinschalten des Netzgerates gezeigt. Die statistischen Schwankungen konnten 
zum Teil als Schwankungen der Temperatur innerhalb des Gerites identifiziert 
_ werden. 

i Fur sehr niedrige Frequenzen betragt der Innenwiderstand R, ~ 1,5- 10-3Q 
und steigt bei etwa 100 kHz auf 0,2 Q. Im Gebiet 0 < # < 20 kHz wird der 
_Innenwiderstand durch die Schaltung bedingt. Bei héheren Frequenzen ist die 


a 
20 
10 

0 b 

ar | 3 
Kh th hk hh hk 2 rT 
Stunden —=— t 
500 ys 
Figur 13 Figur 14 
Drift der Ausgangsspannung des Netzgerates nach a dein Ausgang der Schaltung nach 
Figur 12 nach dem Einschalten. Figur 12 aufgepragter Stromimpuls; 


b Spannung iiber dem Innenwider- 

stand, das heisst an den Ausgangs- 

klemmen. Maximalwert: 0,28 V Spit- 
zenspannung. 


Impedanz der am Ausgang liegenden Kondensatoren massgebend, und R; hangt 
dann wesentlich vom Typ und Fabrikat der verwendeten Kondensatoren ab. 
Das Impulsverhalten des Ausganges entnimmt man aus Figur 14. 

Der Netzteil wurde als Graetz-Schaltung mit Drosseleingang dimensioniert 
(vgl. zum Beispiel BALDINGER [16]). Bei einer Rippelspannung von rund 0,3 V 
Spitzenwert an den Eingangsklemmen der Stabilisierungsschaltung betragt die 
Rippelspannung am Ausgang etwa 0,3 mV Spitzenwert?*). Der Netzteil weist 
einen Innenwiderstand von etwa 2 2 auf. Die gesamte Schaltung arbeitet bei 
einer Umgebungstemperatur von 45°C noch einwandfrei. 

Abschliessend soll erwahnt werden, dass sich die Schaltung nach Figur 12 
weiter verbessern liesse, wobei in erster Linie an eine Verbesserung der Konstanz 


der Referenzspannung gedacht werden muss. 


13) Da die Selengleichrichter im Netzteil den Strom zeitweilig unterbrechen [16], wird die 
Drossel zu Schwingungen angeregt, die sich auf den Ausgang des Netzgerates iibertragen. Im ange- 
gebenen Beispiel treten kurze Stérungen von etwa 5 mV Spitzenwert auf. Diese Schwingungen 
lassen sich prinzipiell beseitigen, wenn man statt Selengleichrichter zum Beispiel Si-Dioden ver- 


wendet. 
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b) Schaltung mit Léngstransistor..Die in Figur 15 dargestellte Stabilisie- 


rungsschaltung mit Langstransistor liefert eine Ausgangsspannung von 10,5.V — 
bis 15,5 V und ist bis 0,5 A belastbar, wobei eine konstante Vorbelastung von — 


50 mA eingebaut wurde. Die Referenzspannung ist wie im Beispiel a) ausge- 


bildet. Der Stabilisierungsfaktor des kompensierten Gerdtes betragt etwa — 


100 150. Ferner weist der Innenwiderstand den praktisch lastunabhangigen 
Wert von 0,6 2 auf, und sein Frequenzgang ist durch einen starken Anstieg bei 


Figur 15 


Stabilisierungsschaltung mit Langstransistor: 
@ 2 Intermetall Z 8 in Serie; @ Philips OC 76; @ Philips OC 16; @ 2 Philips OC 76 parallel. 
Netzteil: Graetz-Gleichrichter mit Drosseleingang, Innenwiderstand rund 20. Rippelspannung 
am Eingang der Stabilisierungsschaltung rund 0,6 V Spitzenwert und am Ausgang rund 3mV 
Spitzenwert. 


etwa 20 kHz charakterisiert. Die Ausgangsspannung zeigt einen Temperatur- 
koeffizienten von rund + 3,5 mV/°C, und dieser liegt damit in der abgeschatzten 
Gréssenordnung. Die maximal zulassige Umgebungstemperatur betragt etwa 
isles 

Dieses Netzgerat eignet sich zum Beispiel fiir die Anwendung in industriel- 
den Schaltungen. Auf die KurzschluBsicherheit wurde verzichtet und als 
Uberlastungsschutz sowohl fiir den Selengleichrichter als auch den Leistungs- 
transistor eine Feinsicherung eingebaut. Uber die Verbesserungsméglichkeiten 
gilt das beim vorherigen Beispiel Erwihnte. 


Anhang 


Wir geben hier Messungen iiber das Temperaturverhalten eines OC 71 
wieder, ohne jedoch niher auf die theoretische Deutung einzugehen. Als Vor- 
aussetzung erwdhnen wir zunachst, dass nach [2] und [3] die Parameter der 
Ersatzschaltung héchstens mit T bzw. 1 /T variieren, im Gegensatz zur expo- 
nentiellen Abhangigkeit des Arbeitspunktes also annahernd als konstant ange- 
nommen werden kénnen (vgl. auch [18] fiir Angaben tiber k,). Wir beschaftigen 
uns also mit der Temperaturabhangigkeit des Arbeitspunktes, 
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_ a) Der Basisstrom variiert bei konstantem Kollektorstrom wie in Figur 16 
dargestellt mit der (absoluten) Temperatur. 

Als Temperatur ist hier die innere Temperatur des Transistors aufgetragen, 

die sich zum Beispiel aus dem Kollektorsperrstrom ermitteln lasst. Der Kollek- 


_ torsperrstrom zeigt in sehr guter Naherung eine exponentielle Abhangigkeit 


36 
102 i" 
233 253 273 293° (313 

—Tr——* | fee oe 

Figur 16 ‘ Figur 17 
Basisstrom J, als Funktion der inneren Tem- Emitter- Junction-Spannung U,; als Funktion 
peratur mit dem Kollektorstrom als Para- der inneren Temperatur mit dem Kollektor- 
meter. Ge-p-n-p-Flachentransistor OC 71. strom als Parameter. Ge-p-n-p-Flachentran- 


sistor OC71. Die Kurven nahern sich mit ab- 

nehmender Temperatur den gestrichelt ein- 

gezeichneten Geraden, welche sich bei T= 0°K 
im Punkt U,; = 760 mV schneiden. 


_von der (absoluten) Temperatur. Der Schnittpunkt der Kurven mit der Ab- 


szisse ist in grober Naherung durch das Uberhandnehmen der Ladungstrager- 
erzeugung in der Basis zu deuten. Bei steigenden Temperaturen kehrt der 


Basisstrom sein Vorzeichen um. 
b) Emitter- Junction-Spannung: In Figur 17 ist die Emitter- Junction-Span- 


nung U,; U,;= Up — KG 


als Funktion der absoluten Temperatur bei konstantem Kollektorstrom auf- 


getragen. 
Den Zusammenhang zwischen U,; und T kann man qualitativ verstehen. 


Fiir gentigend grosse Kollektorspannung ist J, im wesentlichen durch die Eigen- 
schaften des Emitteriiberganges bestimmt, das heisst 


Toate ABIES = const, 7 U,, AE. 
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Mit AE bezeichnen wir die Breite der verbotenen Zone zwischen Valenz- und 
Leitungsband. Man erhalt 


qU.j—-4E _ _ sac ln ty 58 ; 
kT ae Sar Lies q ; 


wobei 6'von J, abhiingen wird. Wie Figur 17 zeigt, gibt diese Beziehung die 
Verhiltnisse qualitativ gut wieder, einschliesslich des extrapolierten Schnitt- 
punktes fiir T = 0, fiir den U,;= AE/q gilt. Die von GARTNER [19] unter gewissen 
einschrankenden Voraussetzungen, aber mit Beriicksichtigung der Temperatur- 
abhingigkeit der Beweglichkeit usw. gerechneten Kurven zeigen dasselbe Ver- 
halten. 

Es lasst sich ferner anschaulich erklaren, dass es Zustande des Transistors 
gibt, bei denen U,,; = 0, aber J, + 0 ist. Nehmen wir zum Beispiel an, dass ein 
Ge-p-n-p-Transistor bei einigen Volt Kollektorspannung einen Kollektorstrom 
von 2mA fiihrt. Zwischen der Lécherkonzentration in der Basis an der Emitter- 


Junction und der Spannung U,,; besteht dann ein wohldefinierter Zusammen- 


hang [vgl. zum Beispiel [21], S. 276, Gl. (8-14)]. Die Gleichgewichtskonzen- 
tration in der Basis wird nun mit wachsender Temperatur zunehmen, das heisst 
bei konstantem J, (konstantem Gradient der Lécher-Konzentration in der Basis), 
die Spannung U,;abnehmen. Bei geniigend hoher Temperatur wird U, ; schliess- 
lich Null werden und sogar sein Vorzeichen umkehren kénnen. Ahnliche Kurven 
wie in den Figuren 16 und 17 sind letzthin von GUGGENBUHL und SCHNEIDER [18! 
publiziert worden. 

Abschliessend sollen noch zwei Bemerkungen gemacht werden: 

1. Unterwirft man einen Transistor nicht zu grossen Temperaturanderun- 
gen, so hat dies Schwankungen des Arbeitspunktes um einen festen Wert zur 
Folge. Diese Schwankungen kann man auf einfache Art erfassen, solange die 
Temperaturabhangigkeit der Parameter des Ersatzschemas nach Figur 1 ver- 
nachlassigbar ist. Man erhalt so zum Beispiel Gleichung (19). 

2. In grober Naherung sind Leistungstransistoren darstellbar durch Parallel- 
schaltung vieler Niederleistungstransistoren™). Insofern beanspruchen vor 
allem die in Figur 17 dargestellten Zusammenhange eine gewisse allgemeinere 
Giiltigkeit; dies um so mehr, als die Darstellung unabhangig vom speziellen 
Basiszuleitungswiderstand R, ist. 


Der Stiftung Hasler-Werke, Bern, haben wir sehr fiir die finanzielle Unter- 


stitzung zu danken, die es erméglichte, die vorstehenden Untersuchungen aus- 
zufiihren. 


14) Allerdings bestimmt bei stark ausgesteuerten Leistungstransistoren im wesentlichen R, die 
Eingangsimpedanz (Ry > Rs); zum Beispiel betragt bei 1 A fiir einen OC 16 Ry ~ 6 Q, wahrend 
Ry < 0,52 ist. 


1 OR 
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Summary 


This paper deals with stabilized voltage supplies employing transistors. A 
theory is given which can readily be applied, and examples of existing devices 
are added. 

With the transistors now available, stabilized voltage supplies may be designed 

to display highly satisfactory characteristics — especially a small temperature 
dependence of output voltage. 

Some ideas concerning the influence of temperature on transistor-action are 
given in the appendix. 


(Eingegangen: 6. Juli 1957.) 


15) Anmerkung bei der Korrektur: Diese Arbeit erschien erst kiirzlich und wurde der Voll- 
standigkeit halber in das Literaturverzeichnis aufgenommen. 
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Finite Difference Methods of Solution of the von Mises 
Boundary Layer Equation With Special Reference to 
Conditions Near a Singularity 


By ANpDREw R. MircHELt and Joun Y. THomson, St. Andrews, Scotland?) 


Introduction 


The von Mises transformation of the boundary layer equations was first 
published in 1927 [1]?). A discussion of this form of the equations was given in 
1938 by PRANDTL [2], who had deduced the transformation before voN MISES 
without publishing it. GORTLER [3] considered the possibility of using this form 
of the equations to solve the problem of incompressible boundary layer flow 
along a semi-infinite flat plate against an adverse pressure gradient. The 
unpleasant singularity of the equation at the plate, however, caused GORTLER to 
abandon the transformation as a means of solving the problem. The singularity 
was avoided by LucKERT [4], who used the von Mises equation to work out the 
flow in the wake of a finite plate lying in the direction of the stream. 

In the present paper it is intended to show that the singularity of the von 
Mises equation at the plate in GORTLER’s problem can be overcome by using an 
appropriate expansion for the velocity in the vicinity of the plate. As a result, 
conditions at the plate can be deduced from the values of the velocity in the 
boundary layer obtained from a step-by-step finite difference procedure. Once 
the difficulty of the singularity has been overcome in the incompressible case, 
it is intended to use the method to obtain a solution of the corresponding 
problem in compressible flow. 


The von Mises Equation and Its Finite Difference Replacement 


The equations of motion and continuity for the incompressible boundary 
layer are in the usual notation 


Ou Ou du, 0?u 
Mon Pay Raper ope 
and 
Ou ov 
Pema, Foerare 


1) The University, Department of Mathematics. 
) Numbers in brackets refer to References, page 36. 
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_Tespectively. In terms of the von Mises variables x, y the above equations 
_ become 


FF Ou du, 0 Ou 
! Ox 1 dx Ow (u ae 
or in non-dimensional notation 
0z 022 
an og (1) 


~ where 


and u, z, x, py now stand for u/U, 2/U2, x/L, yR'?/UL, respectively. The non-' 
_dimensional free stream velocity distribution outside the boundary layer is 
taken to be 


u=1—2%, 


where the semi-infinite flat plate lies along the positive part of the x-axis. 
The four point forward difference replacement of (1) is 
ee 25k 6 [(uz),; << a (% aa ae % 4-1) , (2) 
where 6 = Ax/(Ay)?, and Ax and Ay are the non-dimensional mesh lengths in 
the x- and y-directions respectively. The calculation commences at x = 0-05 
with values of z taken from Howartu [5], and proceeds step-by-step in the 
x-direction. The problem is illustrated in Figure 1. 


0 0:05 


Figure 1 
Illustration of problem. 
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Stability of the Finite Difference Equation 


In a step-by-step calculation using a finite difference equation such as 
(2), care must be taken that round-off errors introduced at any stage do not 
grow with increasing x. This is usually accomplished by chosing a suitable © 
value of the mesh ratio 6. In the present problem the growth of error is governed ~ 
by the equation 


Gara = a + 9[(04); — al? (6 4a — 2&0 + &, 4-1) | ; 
ik Esk (3) 
— 2 [i — % PP (Zner— 2 %4 + pa) + 


1 


The latter is obtained by replacing z in (2) by z + «, and neglecting higher powers — 
of e. In the manner of O’BriEN, HYMAN, and KAPLAN [6], the substitution 


a7 Ax eibk Ay 


ei k = 6 
is made in (3), and the result 
\ 0) : A 
Oe a are [8 uf, sin? Fenn rs Oak tit “i-1)| 
ds 


obtained. The condition for stability is |e*4*| < 1, which leads to 


j< au 4 (4) 
= 2 2 
10 0), — (4) aaa + % aa) 
and 
4s pact 4 apne &y pecs Oe (5) 


Provided conditions (4) and (5) are satisfied at every step during the calculation, 
round-off errors introduced at any stage do not grow and the step-by-step 
procedure using (2) is stable. 


The Finite Difference Calculation 


The values of Ay and 6 are chosen to be 0-1 and 0-25 respectively, and the 
calculation carried out using (2) starts at x = 0-05 with values of z taken from 
HowartH, and proceeds step-by-step with 4x = 0-0025 as far as x = 0-125. 
During the calculation, stability condition (4) is satisfied at every point but 
condition (5) is violated slightly at nodes on y = 0-1 beyond x = 0-08. With 
Ay = 0-1, z is evaluated to four places of decimals and selected values are 
shown in Table 1. The boundary layer thickens from eleven mesh lengths at 
x = 0-05 to eighteen mesh lengths at x = 0-125. 


Te 
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Table 1 


| 0:07 | 0:09 | 0-10 | 0-11 | 0:12 


0-8649 0-8281 0-8100 0-7921 0:7744 
0-6801 0-6800 0-6770 0-6728 0-6675 
0-5004 0-5259 0-5337 0-5393 0:5430 
0-3478 0-3879 0-4025 0:4144 0-4241 
0-2283 0-2736 0-2913 0-3065 0-3195 
0-1419 0-1845 0-2025 0-2183 0-2325 
0-0829 0-1188 0-1350 0-1498 0-1634 
0-0454 0-0729 0-0862 0-0988 0-1108 
0-0232 0:0426 0:0526 0-0625 0-0724 
0-0110 0-0235 0-0306 0-0379 0-0455 
0-0047 0-0122 0-0169 0-0220 0-0275 
0-0018 0-0059 0-0088 0-0122 0-0159 
0-0006 0-0027 0-0043 0-0064 0-0087 
0-0002 0-0011 0-0019 0-0030 0-0046 
0-0000 0-0004 0-0008 0-0015 0-0022 
0-0001 0-0003 0-0006 0-0010 

0-0000 0-0000 0-0002 0-0004 

0-0000 0-0001 

0-0000 


The violation of condition (5) at nodes next to the plate can be explained 
by examining the original equation (1). At any station x, 02/0x is zero in the 
stream outside the boundary layer, negative on the plate, and has a maximum 
inside the boundary layer. As a result, 02/0x is negative in a region adjacent to 
the plate. In this region, however, w is positive and so from (1) it follows that 
02z/dw? is negative. Since difference replacement (2) is of the same order in *% 
and y as (1), it possesses the same inherent properties. Thus in a region next to 
the plate, ‘ : 

; ee > Zep ner — 24,4 + % 2-1) SO, 
and so condition (5) is violated. It should be pointed out that a difference 
equation of order higher than the original differential equation is also a mini- 
mum order replacement of a higher order differential equation, and so the 
higher order difference replacement need not possess the same properties as 


the original differential equation. 


The Unstable Region 


In the previous section, it is shown that a finite difference step-by-step 
procedure with Ay = 0-1 is mildly unstable at nodes beyond x = 0-08 adjacent 
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to the plate. If the mesh length Ay is chosen to be smaller than 0-1 in the 
original calculation, condition (5) is violated at a correspondingly greater — 
number of nodes. In fact, the smaller the mesh size Ay the more precisely is — 
the region of instability determined and the more accurate are the values of 
z at nodes where condition (5) is not violated. Thus the initial step-by-step ~ 
calculation using (2) with prescribed values of Ay and 6 divides the field of the 
boundary layer into two distinct sections, a large region covering most of the field 
where round-off errors do not grow, and a narrow region adjacent to the plate, 
where round-off errors may grow appreciably. The locus 02z/dy? = 0, which 
separates the two regions, diverges slightly from the plate as x approaches © 
the value 0-125. It should be emphasized that this region of instability 
is not a result of 02z/dy2 being infinite, but of 0z/0x being negative, on the plate. 
For example,-if the free stream velocity distribution outside the boundary 
layer is taken to be u,=1-+ x, then 02z/0x is positive on the plate and the 
stable region covers the entire field, although 0?z/dy? is still infinite on the 
plate. In general, the difference calculation is stable or unstable next to the 
plate depending on whether the flow is accelerating or decelerating in the 
main stream. This splitting-up of the boundary layer in the case of an adverse 
pressure gradient into two distinct regions can be compared with the method 
used in 1934 by von KARMAN and MILLIKAN [7], who obtained ‘outer’ and 
‘inner’ solutions depending on whether 02z/dy? is > or < 0. 

As explained in the previous section, the unstable nature of the difference 
calculation in the region next to the plate for the case of retarded flow is a 
direct consequence of the differential equation, and so there seems little hope 
of obtaining a finite difference replacement of (1) for which round-off errors 
do not grow in this region. In the present problem, therefore, an alternative 


procedure is required for finding z at nodes in the unstable region next to the 
plate. 


The Singularity at the Plate 


At the plate, w = 0 and 0z/0x = — 2 (1 — x), and so from (1) it follows that 
0*z/dy? is infinite. This singularity caused GORTLER to abandon the von Mises 
transformation as a means of solving the present problem. Before proceeding 
further with the evaluation of z at nodes near the plate it is desireable to ex- 
amine this singularity in detail. 

Consider at any station x, the expansion 


w= ayh + byt op? + dyrt on, (6) 


where the coefficients a,b,c, d,... are functions of x. This expansion gives 
0°z/d~p? infinite and u (du/dy) = a2/2 at p= 0. The coefficients are not all. arbi- 
trary, but must be chosen to satisfy the following conditions of compatability 
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at the plate; 


a user = — (55), = 20-4) 
Be - “ay (% dy (Hap)} =o. 


It can easily be verified that the above conditions lead to the results 


2 Oz 
Dis 3 a* fore 


See rate 
18 a alta 


and 


C= 
respectively, and so (6) becomes 


2 


, 2 0z a 0z 
= IE EY pe aly (i La lima i Bae ts) 3/2 Dy ea 
die 3a? eee 18 a ike shi (7) 
The corresponding expansions 
4 0z 3/2 1 Oz \2 2 
ae td eo ee Vedios! yo? age oe a 


ae gf" yp = gt yp =, rf ye sce. 


and 
where d’,e’, f’,..., a", e", f”,... functions of a, b,c, d,..., are obtained from 


(7) using z= v2 — uw? and (1) respectively. It should be noted that if, for 
example, u is terminated at p?, the comparable expansions for z and 02z/dx 
‘contain terms up to and including y*. The expansions (7), (8), (9) can be used 
for interpolating w, z, 0z/dx respectively in the vicinity of the plate. 

The basic formula (7) is now used to find du/dy on the plate at a station x. 
LEIGH [8] using the boundary layer equations with u as a function of y, ob- 
tained accurate results for the velocity in the problem of the present paper for 
a small range of x upstream of and including the separation point x = 0-1198. 
LEicu’s results at x = 0-1198 are converted to give w in terms of py and are 
shown in Table 2. From any two pairs of corresponding values of « and y, and 
expansion (7) terminated at y?, a sextic equation is obtained for a. This equa- 

tion in all cases yields two real roots, one positive and one negative, the positive 
root being the one required. The values of a?/2 obtained from different pairs of 
corresponding values of « und y are shown in Table 2, where in each case the 
smaller of the two ‘values of y is used in the tabulation. The results are also 
illustrated in Figure 2 where it is seen that a?/2 (= 0u/dy) extrapolates to zero as 


i 
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yw tends to zero. Thus provided the expansion (7) is used, the singularity does . 
not prevent an accurate value of 0u/dy being obtained on the plate, the accuracy 
depending only on the values of the velocity in the field. Of course, the closer 
to the plate the values of the velocity are known, the easier is the extrapolation 
and so the more accurate is the value obtained for du/dy on the plate. It can 
be stated that if a set of values of the velocity is given in the field in terms of 
both wy and y at any station x, then u (du/dy) (= 0u/dy) on the plate can be 
obtained just as easily and with as much accuracy using expansion (7), as using 


\ 


Table 2 

y u a?/2 

0:2 0-4780 
0-1 0-3200 0-254 
0-05 0-2075 0-187 
0-025 0-1325 0-164 
0:0125 0-0830 0-126 
0-006 25 0:0548 0-118 
0-003125 0-0345 0-088 
0-0015625 0-0210 0-065 
| 0-000 78125 0-0125 0-045 


Figure 2 
Separation point at x = 0-1198 from Lr1Gn’s results. 
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_ the standard expansion for w in terms of y, in which no singularity exists on 
the plate. 


Evaluation of the Velocity in the Unstable Region 


As shown previously, an initial step-by-step calculation, using (2) with 

_ prescribed values of 6 and Ay, produces a narrow region adjacent to the plate 

_ where round-off errors may grow. Accordingly an alternative method is required 
for finding the values of the velocity in this unstable region. 

In the actual calculation carried out with Ay and 6 equal to 0-1 and 0.25 
respectively, the stability condition is violated slightly at nodes on wy = 0-1 
beyond x = 0-08. The method used to find revised values of z at these nodes 

_ is now explained. 

At x = 0-0825, 02/0 is known at y = 0-2, 0-3, 0-4, and so 0z/0x at y= 0-1 
can be obtained by using (9) terminated at w?. However, z is known on y = 0-1 
at x = 0-0775, 0-08, and so z can be obtained at x = 0-0825 from the three point 
backward difference formula 


0z 
Foleo = 200 (2.0775 — 4 20.08 + 3 20.0825) « 


This procedure is repeated to find z at x = 0-085 and so on up to x = 0-125. 
The new values of z obtained at the nodes on y = 0-1 beyond x = 0-08 differ 
slightly from the original values (see Table 3), showing that only mild insta- 


Table 3 
| z 
x 

Original | Corrected 
0-080 0-6812 0-6813 
0-085 0-6809 0-6810 
0-090 0-6800 0-6803 
0-095 0:6787 0-:6793 
0-100 0-:6770 0:6776 
0-105 0-6751 0:6767 
0-110 0-6728 0-6750 
OATS 0-6703 0:6725 
0-120 0:6675 0-6700 


bility exists in the original step-by-step calculation. It should be mentioned 
that the changes in the values of z at nodes on py = 0-1 cause minor alterations 
in z at nodes on py = 0-2. 

The method described can be continued to find the velocity within the 
region 0-05 < x <0-125,0< p < 0-1. At « = 0-0525, 0z/0x is known at yp = 0, 
0-1, 0-2, 0-3, and so a formula for 02/0 is obtained from (9) terminated at yp. 


ZAMP IX/8 


ee ee eee 
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From Howartw’s results, z is known at « = 0-0475, 0-05, and so z is evaluated © 


at nodes within the range 0 < y < 0-1 on x = 0.0525 from the three point 


backward difference formula for (02/0%)o.9525. This procedure is repeated to find — 


z for 0 << y<0-1 at x = 0-055 and so on up to x = 0-125. Of course, once the 
separation point (du/dy = 0 on the plate) is obtained by the method described 
in the previous section, the values of the velocity at stations x downstream of 
the separation point can be discarded. Although it is considerably simpler to 
obtain the velocity in the region 0 < w < 0-1 directly from (7), the values are 
in fact less accurate. This is due to the direct method not making use of Ho- 
WARTH’S values of z at x = 0-0475, 0-05. 

An alternative method for finding the velocity at nodes within the region 
0<y< 0-1 is now described. At a station x, 0z/0x is known at y= 0, 0-1, 
0-2, 0-3, and so from (9) terminated at w?, an expansion is obtained for 0z/0x 
which can be used within the region 0 << wy < 0-1 at the appropriate station x. 
The velocity at y = 0-05 at a station x is now obtained from 

(Fo )o.ng == 400 (2 — A008)” (oa — 2 Ao.05 + 20)» 
which is the central difference replacement of (1) at y= 0-05. In a similar 
manner, the velocity at y = 0-025 is got from 


Oz 2 2 y 
hora: = 1600 (uy a ee (20.05 Es 20.025 = Zo) = 


and so on. This process is carried out at the section x = 0-12, and the values of 
the velocity obtained at y = 0-05, 0-025, .... , 0-:0015625, are shown in Figure 3 


0:20 


0-15 


0-10 


—— LEIGH'S Values 
e Calculated values 


Figure 3 
Comparison of results at x = 0:12. 
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along with the corresponding values obtained by LeicH at x = 0-1198. The 
agreement is surprisingly good considering the large increments Ay used in the 
early stages of the present calculation. 
_ More accurate results can be obtained by this alternative method if 
(1) the initial step-by-step calculation is carried out using a value of Ay 
less than 0-1. 
(2) at each station x, the range of y within the unstable region is divided 
up into more than two equal parts. For example, if the range 0 < y < 0-1 is 
divided up into four equal parts, a central difference replacement of (1) with 
» Ay = 0-025 is used at wy = 0-025, 0-05, 0-075, and so it is required to solve three 
equations for the three unknowns 29.995, 20.05, 20.075: Lhe range 0 < y < 0-025 is 
then divided up into four equal parts, and the above process used to find 
2.00625» 20.0125» Z0.01875, aNd so on. 

(3) expansion (9) is used with values of d”, e”, ¢” revised at each stage. For 
example, if the range 0 < wy < 0-1 is divided up into four equal parts at three 
neighbouring stations x, and z is obtained at each node, then 0z/dx at the 
station x furthest downstream can be calculated from a three point backward 
difference formula at y = 0-025, 0-05, 0-075. A revised expansion for 0z/0x, 
which can be used in the range 0 < wy < 0-025, at the appropriate station x is 
obtained from (9) using the values of 02z/0x calculated above. 


Concluding Remarks 


(1) The primary object of the present paper is not to obtain an accurate 
value for the separation point in the problem of boundary layer flow along a 
flat plate against an adverse pressure gradient, but to explore generally the 
possibility of using finite difference methods, based on the comparatively 
simple von Mises equation, to solve boundary layer problems. The values of 
velocity computed are generally in good agreement with those obtained by 
Lr1IGH, while separation, according to our calculation, occurs slightly down- 
stream of the accepted position, « = 0-1198. The accuracy of the solution 
obtained from the numerical calculation depends on the mesh size used, which 
in turn is governed by the computing facilities available. 

(2) Previous investigators did not consider the von Mises transformation 
as a means of solving the present problem because 0z/dy? is infinite at the 
plate. However, provided an expansion like (7) is used, which incorporates the 
conditions of compatability at the plate, there is no difficulty in obtaining a 
value for 0u/dy on the plate from the computed values of the velocity in the 
boundary layer. oe 

(3) The errors in the initial step-by-step calculation in the region next to 
the plate are caused not by 07z/dy® being infinite on the plate, but by 0°2/dp? 
being negative in the vicinity of the plate. This violates the stability condition 
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for the initial step-by-step calculation, and so an alternative procedure is 
adopted in order to find the velocity in the region next to the plate. 

(4) An apparent disadvantage of all step-by-step techniques based on finite 
differences for solving boundary layer problems is that values of the velocity 
are required at some station x in order to start the calculation. In the present 
example, Howartu’s values are used at x = 0-05. This disadvantage can be — 
overcome to some extent by commencing the calculation at the nose of the © 
plate (x = 0) with the values z= 1 at y=0 and z=0 for py > 0. This was 
tried out in the present problem with Ay and 6 equal to 0-1 and 0-25 respec- 
tively. At each station x, the computed values are compared with the theoretical 
values of HowarTH, the agreement between the two sets of values improving 
as the calculation proceeds. In fact, by the time the station x = 0-05 is reached, 
the values of z are within three per cent of HOWARTH’s values. This agreement — 
is of course improved if a smaller value of Ay is used in the calculation. 

(5) The velocity is obtained as a function of y by evaluating numerically 


ae 


¥ 
y= f dy/uw at each station x. The integrand is infinite at the lower limit and 
0 


so it is necessary to evaluate y in the form 


Uy 
ee: 
0 


where u, is the velocity corresponding to y,, an arbitrary small value of y, and 
u (du/dw) takes the value a?/2 on the plate. 

(6) Finally, it is hoped to use the methods outlined in the present paper to 
obtain a solution for the corresponding problem in compressible flow. The 
governing equations are of course much more complicated, and as a result a 
greater amount of numerical work is necessary in order to obtain a solution. 
The expansions (6), (7), (8), which are based on incompressible boundary layer 
flow, are sufficient provided that the coefficients in the expansions are obtained 
from compressible flow calculations in a region near the plate where the 
velocity is sufficiently small. 


a 


J 


du 
at 
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Zusammenfassung 


Friihere Forscher haben die Misessche Gleichung als Mittel zur Lésung des 
Problems der Grenzschicht an der langsangestr6mten ebenen Platte mit Druck- 
anstieg nicht verwendet wegen ihrer unangenehmen Singularitaét an der Platte. 
In der vorliegenden Untersuchung wird eine Geschwindigkeitsverteilung berech- 
net, welche die Plattenbindungen beriicksichtigt und folglich den Hautreibungs- 
wert durch die Geschwindigkeitswerte in der Grenzschicht berechnen lasst. 
Diese Werte werden ermittelt durch ein Fortsetzungsdifferenzenverfahren auf 
der Grundlage der Misesschen Gleichung sowie durch ein modifiziertes Differen- 
zenverfahren in der Nahe der Platte. Dieses letzte Verfahren ist nétig wegen der 
Instabilitat der urspriinglichen Fortsetzungsberechnung in der Nahe der Platte. 
Angesichts der verhaltnismassig geringen Rechnungsarbeit stimmen die Ergeb- 


mnisse gut tiberein mit den Lésungen anderer Forscher, die andere Formen der 


Grenzschichtgleichung benutzten. Die Genauigkeit der in dieser Untersuchung 
erreichten Ergebnisse lasst sich natiirlich noch vermehren durch Anwendung 
einer kleineren Schrittweite bei der ersten Fortsetzungsberechnung. 
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1. Introduction 


' 
+ 


‘ 


Although there is some difference of opinion on the approach to the theore- — 
tical and experimerital phases in studying the hail control problem (cf. e. g. — 


WEICKMANN [1]?) and DESSENs [2]), most will agree that proof is a matter of 
theory and verification is a matter of experimentation and these must go hand 
in hand in the development of science. While theoretical development in mete- 
orology presents no peculiar difficulties, experimentation presents a major 
difficulty resulting from the fact that nature ultimately controls all experiments 
in the atmosphere. WEICKMANN [3] aptly points out the main factor in this 
difficulty as the necessity of accounting for the natural variability in time of 
meteorological phenomena. With proper understanding of meteorological data, 
however, we believe this difficulty can be overcome to an extent that allows 
valid conclusions to be drawn from experiments now being conducted in the 
atmosphere. Certainly this can always be accomplished with properly designed 
experiments. 

Although it is not our present objective to evaluate the results from a 
specific hail project, it would be advantageous to work within practical limits 
of hail occurence as might be associated with an actual hail suppression project. 
In order to develop and apply the methods we have in mind, the suppression 
project target should contain at least one but preferably more than one long 
record meteorological station with good hail records. Examination of our records 
showed no United States project meeting this requirement. The situation in 
Europe, however, was better, for hail suppression has been practiced much more 
extensively there. 

There are, at present, at least two experimental operations to suppress hail 
underway in Europe: One is in France and has been described by H. DESSENS [4] 
and the other is in the Magadino Plain in Switzerland under direction of 
R. SANGER, H. Staus and J. C. Tuams [5]. The latter operation is especially 
well controlled and seemed most appropriate as a site around which we could 
develop our methodology. 

A review of the Swiss reports [5] showed that there are three good meteoro- 
logical stations in the Magadino Plain observing hail occurence. Although, as 
pointed out to us by Prof. SANGER in correspondence, these stations are not 
properly located to sample the hail suppression operations as conducted 
adequately, this was not critical since our main objective was to develop 
methods, and the station arrangement and records were good for this purpose. 


Long, carefully kept records were readily available in the Swiss Meteorological 
Office publications [6]. 


*) Numbers in brackets refer to References, page 63, 
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2. The General Problem 


It would appear that hail suppression efforts might be evaluated either by 
using some meteorological variable or some related variable such as damage to 
crops. In this paper we discuss only the use of a meteorological variable, the 
related variable being left to another study. 

The only hail variable on which historical data are ordinarily available is 


_ the frequency of days with hail, and it suits our present purposes best to use the 


number of hail-days per year as our variable. From a consideration of the 
stochastic aspects of hail occurence, it is seen that the number of hail-days will 


_ be proportional to the total area of hailstorms. This depends on the area of 


individual storms as well as on the frequency of storms. If it is assumed that 
hail suppression attempts could either prevent a storm from producing hail 


‘altogether or reduce the area over which hail occurs, then the total area of 


hailstorms is in either case reduced and hence the frequency is reduced. Pre- 
sumably, then, if the conventional frequency of hail-days is reduced, the damage 
from hail would also be reduced. 

- It must not be forgotten that hail damage might be reduced by the sup- 
pression operation having affected some physical properties of the hail such as 
size, intensity, and hardness. Hence, while it might appear that damage is 
related to frequency, it might also be affected by other factors so that frequency 
would seem to provide a conservative measure of hail damage. 

Since the covariance method of using a control area in conjunction with the 
target proved so useful in evaluating seeding for. precipitation augmentation, 
it was natural to give it first consideration for hail suppression. There were, 
however, some already existing questions on its applicability. In the first place, 
hail is associated with summer convectional storms for which the correlation 
with distance on precipitation drops off rapidly with increasing distance 
because of the sporadic nature of convectional precipitation. Second, since the 
hail pattern is usually a small part of the rain pattern, it tends to be even more 
sporadic than rain and hence have a more rapidly decreasing correlation with 
distance. This is borne out by correlations found between the Magadino-Plain 
stations Locarno-Muralto, Bellinzona, and Lugano. These stations are located 
approximately on an isosceles triangle with Locarno-Muralto and Bellinzona 
on the base which runs roughly West-East and Lugano to the South at the 
apex. Locarno-Muralto and Bellinzona are about 10 miles apart and each is 
about 14 miles from Lugano. Even with these short distances between stations 
the correlations coefficient between Locarno and Bellinzona is only 0:5318, 
between Locarno and Lugano 0-3085, and between Bellinzona and Locarno 
0-3971. These stations are all in the target area. If we went outside of the 
target area, as we must do for a control, the correlations would be still lower. 
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Since 1 — v2 is a measure [7] of the advantage gained from a control, it can 
be seen that for 7 = 0-50 the gain is only 13%. Consequently, it is nearly 
useless to employ a control where correlations are near 0-50 or less. There may 
be some instances where a satisfactory control may be found, but this is much 
more unlikely for hail than it is for rainfall and hence we must look to other 
tools for effective evaluation methods: 

As with precipitation [7] if the seeding effect had been very large, it would 
have been very noticeable and hence there would have been little need for a 
sensitive evaluation method. With hail, however, the situation is perhaps even 
more difficult than rain, for the effect is not as easily observed and tends to be 
more obscured by the natural variability of weather. This is no doubt largely 
due to the lack of a more complete understanding of the theory of hail formation 
as well as to a lack of engineering techniques necessary to be able to transmit 
effective quantities of nuclei to proper regions of the atmosphere at will. As is 
often the case where benefits from an operation could be large compared to 
costs, there is a great urge to bypass difficult, and therefore slow to develop, 
theory with application. Immediately there is also then an urge to evaluate 
the effectiveness of such operations. This is the situation at the present time. 

In such evaluation one is not only faced with the variability in the skill of 
the operator but with the variability due to the uncertainty of the theoretical 
effect. These, together with the great weather variability in atmospheric opera- 
tions, make the evaluation of hail suppression at present largely one of statistical 
analysis based on whatever physical knowledge is available. Of course, since 
evaluation must inevitably use the results of observation and the atmosphere 
has a certain uncontrollable. random vatiability, evaluations must always be 
partly statistical in nature. This inherent random variability at the same time 
limits the sensitivity of purely statistical methods and so we must look ulti- 
mately to improved physical models or theory for major improvements in the 
sensitivity of evaluations. For the present we must make best use of the data 
available to us. 

Since the covariance method using the classical method of testing hypo- 
theses is of little aid, and we need a method of testing which makes the best use 
of observations on the target alone, we decided to devise a so-called sequential 
method of testing the hypothesis that hail suppression was effective. The 
sequential test of hypothesis can be shown to reach a decision on the average 
with fewer observations than the classical test. See WALD [8], p. 54. This partly 
compensates for the lack of a control area. 


3. Tests of Hypotheses 


_Every experiment, if carried to a proper conclusion, ends with a decision 
being made as to whether to reject or accept the hypothesis which the experi- 
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‘ment was meant to test. In many fields, experiments can be controlled or 

‘replicated to such an extent that the probability of committing errors in 
making decisions is so small that the correct decision seems quite obvious at 
‘the end of the experiment. This type of experiment is so prevalent in the 
physical sciences that many do not realize that a probability of coming to a 
wrong decision even exists; and so the need is not recognized, nor in fact, is 
there much need for systematizing the decision-making procedure. 

In some fields, such as meteorology, replication is costly and time consuming 
because the trials of an experiment must be made in a time sequence of nature’s 
choosing and the possibilities of control are severely limited. Here the proba- 
bility of making a wrong decision is not often negligibly small. Hence a knowl- 

edge of this probability becomes important, and a systematized decision- 
making procedure is indispensable to the making of enlightened decisions. This 
is provided by the tests of hypothesis developed by statisticians. 

In simplest terms, the statistical test of hypothesis consists of an hypothesis 
and a rule for accepting or rejecting that hypothesis on the basis of a set of 
experimental trials called a sample. A sample clearly can only lead to making 
a decision which is correct with a certain probability. There are chances then 

that two kinds of errors may be committed: Type 1, the hypothesis may be 
rejected when it is actually true or type 2, it may be accepted when it is false. 
In addition, sample size must be involved since this must affect the confidence 
we place in our decision. Thus the basic elements of the decision-making 
procedure may be said to be the two types of errors and the sample size or 
number. 

The two forms of test we have under consideration here involve fixing two . 
of these elements and allowing the third to vary. For the classical test the 
sample size and the first type error are fixed, and the second type error is 
minimized [7]. One test is then made on the entire sample. For the sequential 
test the two types of errors are fixed and the sample size is allowed to vary, 
i. é., the sample is tested progressively as it accumulates. After each trial a 
computation is made which provides one of three results: accept the hypothesis, 
reject the hypothesis, or make another trial. Owing to the fact that the proba- 
bility of accepting or rejecting the hypothesis approaches one for a proper test, 
the process always terminates in one of two decisions—acceptance or rejection 
of the hypothesis. 

In order to arrive at a test procedure the hypothesis we have been discus- 
sing must be expressed in numerical terms. Suppose we are concerned with a 
parameter « being less than or greater than C. Then one hypothesis of interest 
might be expressed as follows: yw is less than C. In the hail problem if C is the 
mean hail frequency before seeding, we would like to be able to accept or reject 
ye < C on the basis of a series of trials. To do this we need a rule for accepting 


or rejecting our hypothesis. 
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‘ 


Since, as we have seen, the acceptance or rejection of an hypothesis, say H, 
involves probabilities of being correct and in error, the rule must entail a 
population of values of w. Thus in general we accept the hypothesis H if aq 
sample value mu is a member of the population associated with H, and we reject — 
H if wis not a member of the population of H and hence is a member of some 
other or alternative population, say H’. 

Except in very special instances the populations H and H’ overlap, so for 
any sample value of « we merely know its position with respect to H and H’ 
in terms of probability. Hence a reasonable rule would be to favor the accept- 
ance of H when the probability that ~ came from population H is high and to 
favor rejection of H when the probability that ~ came from H is low. To 
determine such probabilities we must have some knowledge of the frequency 
distributions on H and H’. These must be obtained from the climatological 
record. 

In a previous paper [9] we developed the application of the Poisson and ne- | 
gative binomial distributions (n. b.) to hail frequency series. We reasoned 
there that for low average frequency of hail occurrence successive hailstorms 
are independent and hail frequency tends to follow the Poisson distribution in 
which the variance is equal to the mean. With higher average hail frequency 
the hailstorms occur closer together resulting in some pairing or correlation 
which makes the variance of annual frequencies. exceed the mean annual 
frequency—a departure from requirements for a Poisson distribution. In these 
circumstances we found the n. b. or Polya distribution, which has a variance 
exceeding the mean, fit the hail series very well. Here, however, our interest is 
. In the total hail frequency for a network of several stations in the target area 
since this would seem to be a more efficient indicator of hail suppression. 

It is well known that the sum of independent Poisson series, such as those 
from several hail observing stations, would also be distributed in a Poisson 
distribution with mean equal to the sum of the means of the several series. 
This implies also that the variance of this series would be equal to the sum of 
the means. If the stations are in close proximity, the station hail series will not 
be independent but will tend to be positively correlated. Hence, while the 
means will again add to produce the mean total, the variance of the total will 
be greater than the sum of the variances by a covariance which is proportional 
to eg Sons between the stations. Here the"departure from the Poisson 
‘series takes place in a simil it di 
low hail a aie to ane ite oe eh a eo Se aa —— ae 
we have substituted a space variable for a ti neste ae : 

au me variable. Although this is a 
special case and convolution does not take place generally in this manner 
ceealeles rie sti ae combine to form an n. b. distribution except 

; perience shows that n. b. fits the distributions 
of totals quite well. We shall therefore assume that the principles applied 
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through time also apply over space (area) and employ the same distributions 
for network totals of hail frequency as we have for point totals. 


4. Estimates of Distribution Parameters 


Before we can undertake development of specific tests, we shall need 
estimates of the distribution parameters both on the hypothesis and its alter- 
native. Our test series will consist of the total frequency of hail-days at the 
three stations: Lugano, Locarno-Muralto, and Bellinzona. Although Locarno- 
Monti is a main station in the hail-suppression project for the Magadino Plain, 
we could not use it because there was only a short record available to us. It 
also appeared to us that the effects of hail suppression, if real, might be less 
pronounced at Monti although this did not affect our choice of stations. The 
three stations were chosen for no other reasons than that they had the only 
records that met our requirements in the target area. The choice was made 
before any computation or testing was done and so it was absolutely inde- 
pendent of the methods or results later obtained. 

It is well known that the sequential test is sensitive to non-randomness in 
the climatological series. Such non-randomness we reasoned could come from 
two causes: (a) long term fluctuations in the climate which most agree would 
be relatively gradual and (b) changes.in observational technique (regulated or 
unregulated) which could either be gradual or abrupt. The hail-occurrence 
observation being aperiodic is a difficult observation and in some instances is 
subject to unregulated change from observer to observer. It is also subject to 
some regulated variations caused by changes in rules, e. g., ruling that a hail- 
storm extending past midnight or the end of the day is two hail-days as against 
one. Clearly, this should be one hail-day since only one hailstorm is involved. 
If two occurrences are assigned to such an event the hail-days are erroneously 
inflated and the number pairs increased which artificially emphasizes the Polya 
effect. Here is an excellent example of the difference between the requirements 
for climatological and synoptic observations. Notwithstanding the possible 
weakness in the hail observation, every indication we have is that the hail 
records at these Swiss stations are of good quality and have been little affected 
either by regulated or unregulated changes. 

Henceforth we treat only the total annual hail frequency at the three 
target stations Lugano, Locarno-Muralto, and Bellinzona. Complete record for 
their total frequency of hail-days was available to us back through 1888 with 
the exception of 1896 for which the total was incomplete. Hail-suppression 
operations began in the Magadino Plain in 1948 and have been conducted in 
two phases: The first five year’s operations were with rockets 1948-1952; the 
next three year’s operations and still continuing are with ground based Agl 
generators ([5], No. 6). We therefore had an undisturbed record which begins in 
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1888 and extends through 1947 excluding 1896. Of this, use was also made of 


the period 1925-1947 as we shall see presently. 


: 


To support our contention that the two series are essentially unaffected by 


climatic and artificial change and are therefore random, we have made simple 
run and pair tests. The entire record 1888-1947 gave 29 runs above and below 
the median which is well enclosed by the 0-90 acceptance region (24, 37). The 


same period showed 31 pairs above and below the median which is below the ~ 


0-05 limit 35. The shorter period 1925-1947 had 10 runs which is on the middle 
of 0:90 acceptance region (6, 15). For the pair test this period gave 13 pairs 
which is below the 0-05 limit 15. Though we do not maintain these are powerful 
tests for randomness, they do confirm the opinion stated above. 

We view these results as justification for proceeding with a sequential test 
although we took one further precaution against trend in the record. This was 
to use the parameters for the period ending in 1947 with greater than 20 years 
record and having the least mean. This was the 1925-47 period which has a 
mean of 4-78 hail-days as compared to a mean of 5-31 for the full record 
1888-1947. It was also chosen without knowledge of the outcome, and before 
any tests of hail suppression were made with a view to forestalling arguments 
that the mean hail-days had decreased naturally up to the time the experiments 
began. Thus we have chosen a mean on the low side with respect to the long 
time experience and if suppression operations gave some indication of even 
lowering this, it would seem worthwhile to at least continue the experiments 
with greater emphasis on developing the physical theory and the engineering 
aspects of the operations. 

N. b. distributions were fitted to the periods according to methods described 
in [9]. In both instances the method of moments was adequate. Results are 
shown in Table 1: 


Table 1 
Distribution Statistics 


Period | x | # | p | k 
1925-1947 4-78 6-632 0-3875 12-33 
1888-1947 5-31 7-630 0-4381 12-11 


Here x is the arithmetic mean, s? the variance, and p und k the parameters 
of the negative binomial distribution. For the longer record SUKHATME’s 7? 
test showed the variance definitely to exceed the mean so that the n.b. gives 
the better fit. For the shorter period the results are indifferent and the Poisson 
distribution gives as good a fit as the n.b. Since our test is based on the para- 
meters of the shorter period, we use this lack of preference for n.b. against 
Poisson distribution as a basis for developing tests on both distributions. 
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5. The Sequential Test 


The sequential test as we have seen is a type of test of hypothesis suited 
to deciding whether to accept or reject the hypothesis that a parameter esti- 
mated from a sample is different from some preassigned value of that para- 
_ meter. Suppose C is the value of a parameter for some.undisturbed hail regime 
and D is a value of that parameter which would be considered a successful 
reduction from C if attained by hail-suppression operations. A sequential test 
_might then be made during hail-suppression operations to decide whether or 
not the parameter had been reduced to less than D. Such a test would be made 
in annual stages on hail frequency as we propose to do later. 

Inasmuch as many meteorologists will not have encountered sequential 
methods before, we shall present the developments in some detail even though 
_ much of this is already in the literature. After we had completed the analysis 
presented here, we found that G. B. OAKLAND [10] had performed a similar 
analysis in a quite different field. With some modification of the hypothesis 
this served as an excellent check of our own work. We give the analysis for 
both the Poisson and the n.b. distributions realizing that the n.b. may have 
much wider although somewhat more difficult application. 

The basic tool for the development of sequential as well as the classical 
tests is the likelihood ratio ([11], pp. 249 and 365) defined by 


f1(41) fi(%2) ++: Am) 
camer ACANACA EAC ” 
where /,(x) is the probability density function when Hy is true, and /,(x) the 
probability density when H, is false. (x, %2, ..., %m) is a random sample. (1) is 
called the likelihood ratio because the products in the numerator and denomi- 
nator are likelihoods as defined by R. A. FisHer. All the common tests of 
hypothesis developed by FIsHER can be derived by means of equation (1). 

For the classical or fixed sample-size test we choose an interval (0 < 4 < A) 
called the critical region so that the probability of A being on (0, A) is «. This 
is the probability of a type 1 error or the error of rejecting Hy when it is true. 
The likelihood ratio ensures choice of the best critical region which is one for 
which the probability of rejecting Hy when it is false is a maximum or the 
type 2 error f, the probability of accepting Ho when it is jalse is a minimum. 
The probability of rejecting Hy when it is false is called the power of the test 
and is 1 — f. The test is made by drawing a sample (%, Xap vey Xm) and com- 
puting 4 by (1). If Ais on (0, A), Hy is rejected; if A > A, Hy is accepted. 

For the sequential test we have seen that « and f are fixed, but the sample 
size is allowed to vary. Thus /,, is computed for each term of the sequence 
n =1,2,..., m. An interval (B <A, < A) is established on the basis of pre- 
scribed values of « and , such that for each m beginning with m= 1 if-A, 18 
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in this interval J, ,, is computed, and if /,, is not on the interval the sequence — 
is terminated. If further /,, << B, Hy is accepted; if A,, = A, Ho is rejected. 
Thus at each term of the sequence 1, 2, ..., m a test is made which gives either 
the decision to accept or reject H, or to take another sample. Since we require ~ 
a two decision test the process must terminate. WALD [8] has demonstrated 
that whatever the form of /(x), the sequence must terminate with the accept- 
ance or rejection of Hy. He has also given much other ingenious analysis which 
will enable us to measure the effectiveness of the test as well as the expected 
value of m for termination of the sequence or reaching a decision. 

The sequential test is in some respects simple1 to develop than the classical 
test. For example, establishment of the limits on /,, is done very easily. For 
convenience the logarithm of the likelihood ratio is usually employed because it 
reduces products to sums which are easier to handle. Taking logarithms of (1) 


we have 
m 


dD, % = 1OgAm, (2) 


t=1 


where z; represents the individual ratio for each sample. The inequality for 
the test correspondingly becomes 


log B < Siz; <logA. (3) 
i= 


The essential simplicity of the test arises from the fact that as WALD [8] has 
shown A and B to be closely approximated by the relations 


i= p | 
Awe and Base, (4) 
For all practical purposes these may be accepted as equalities as they will be 
here. These four equations as we shall see are the four basic relationships for 
devising sequential tests. 

As we have agreed previously, our measure of the effectiveness of hail sup- 
pression will be based on whether the hypothesis that the mean total hail 
frequency at Lugano, Locarno-Muralto, and Bellinzona has been decreased 
can be accepted and with what risks. When the decision to accept or reject 
such an hypothesis is based on the methods we have been discussing, it is 
necessary to employ the device of setting up some plausible alternative hypo- 
thesis which, if accepted, is equivalent to rejecting the hypothesis that there 
was no decrease. The problem is relatively simple if the distribution has a 
single parameter like the Poisson, for then one needs only to assume some 
value for the mean hail frequency based on past experience or on the desired 
magnitude of the effect. If, however, the distribution has two parameters like 
the negative binomial, it is necessary to consider a composite hypothesis or 
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reduce the latter to a simple hypothesis by assuming a constant value for one 
of the parameters. We shall attempt the latter. 
_ We have observed in connection with the evaluation of precipitation 
augmentation that the correlation in storm precipitation amounts between two 
stations A and B in the same climatic regime is very much lower in convec- 
tional storms than in cyclonic storms. This is due to the fact that the cyclonic 
precipitation pattern being extensive is much more likely to cover both sta- . 
tions than the small convection pattern. In fact, the stations do not have to 
be very far apart before a single convection cell cannot possibly cover both 
stations. When we consider the totals of storms for the convectional season or 
a fixed period in it, we found the correlation higher. Further examination 
showed that this increased correlation was due mainly to the correlation in 
frequency of occurrence. This is not caused by the same cell covering both 
stations, but by the fact that increased frequency tends to be associated with 
weather situations that affect the whole climatic region. Consequently weather 
situations with a large number of convection cells tend to raise the frequency 
at individual stations. The higher the areal frequency in such situations the 
greater will be the chance that stations A and B will be covered by cells and 
hence the higher the correlation between A and B. It follows that we would 
expect that the correlation between A and B would be proportional to the 
mean frequency at A and B. Since hail is a typical form of precipitation from 
a convective cell, hail frequency could be expected to behave like precipitation 
in general. We were able to test the proportionality of frequency and corre- 
lation for Denver (Colorado) and Amarillo (Texas) where the correlation between 
annual hail frequencies was 0-45. When only those pairs of frequency were 
selected whose Amarillo frequency was below a certain value, the total fre- 
quency was reduced to 0-62 of its former value and the correlation to 0-67 of 
its completed series value. Thus a reduction in the hail frequency caused a 
proportionate reduction in the correlation coefficient so the ratio w/v remained 
nearly constant. 

Consider now two stations A and B in the same climatic regime whose hail- 
day series are each distributed in Poisson distributions with mean p. Since 
the variance at A and B will also be y, if there is a correlation 7 between A 
and B, the variance of the total frequency will be 


O72 te ZF Me. 
If we fit a n.b. distribution to the total frequency of A and B, we find by 


substitution in equation (9) of [9] that 


2 
i ee Meee 
te 


i 
In general for / stations tatath’siprp, 


j 
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where 7 is the mean correlation for the / stations. This yields by a similar 
substitution 
a ee aa 
Sl (Me a 


Thus if as hypothesized above 7 is proportional to mw, reduction in frequency 
by hail suppression should leave the ratio u/r constant and hence & will not 
change. This will be our hypothesis in using the n.b. in the sequential test. 
From the relationship 
b=kp (5) 


for the n.b., it follows that if & is constant any change in yw results from a 
change in ~. Hence the composite hypothesis is reduced to a simple hypothesis 
on ~ by assuming k constant. This is the same assumption OAKLAND [10] made 
on the basis of quite different observational material. It makes the test on p 
equivalent to a test on w which is what we desire. 

Referring back to the table of Distribution Statistics we see that the period 
1925-1947 chosen for its low mean has nearly the same value of & as for the 
full record. This agreement has little significance because the standard error 
of k is large, although there is some satisfaction in knowing that k for the 
shorter record is very close to a better estimate obtained from a longer record. 
The test on the Poisson distribution will also be made, of course, using the 
mean directly. We proceed now to the development of the tests. 

We shall speak generally at first in terms of a parameter yw which could 
be either # of the n.b. or the mean of the Poisson distribution. Suppose that 
if wis near wu’ we are only indifferent whether we accept or reject a hail-sup- 
pression operation as effective. If however wu < y’, we would tend to want to 
accept the hypothesis Hy either that seeding was effective or that it would be 
worthwhile to conduct more stringent tests based on more extensive data. If, 
on the other hand jy > yu’ we would tend to accept the hypothesis H, that the 
operation produced no worthwhile effect and we would want to discontinue 
operations or change to another system which had greater possibilities of 
being effective. More specifically our preference for accepting Ho, or hail sup- 
pression, would increase the farther 4 was below yw’ and our preference for 
rejecting hail suppression the farther « was above py’. Hence it will be possible 
in general to find two values wy and jy (Wo < pa’ and 4, > pw’) such that rejection 
of seeding effect is considered an important error when [eS wy and acceptance 
of seeding an important error when pu => [;. We may say that the zone ot prefer- 
ence for accepting seeding as effective is 4 < po, and the zone of preference 
for rejecting seeding as being effective is | 2 wy. The zone of indifference 
would be all values of « between jy and fr- 

Clearly the choice of jg and jy is not a statistical problem; neither is the 
setting of the risks of making the errors referred to above. If the cost of seeding 
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_ is cheap compared to the potential benefit, or if we are making an exploratory 
experiment on a highly desired effect, we should want a high probability of 
finding the effect when it is real and would more readily tolerate a larger 

probability of accepting the effect when it is not real. Thus in general for an 
exploratory operation, the probability of rejecting seeding should not exceed 

_ some small value « when mS po, and the probability of accepting seeding 

when uw < p, should not exceed a value f with 6 = a. 

Without any knowledge of possible results we chose wu, equal to the mean 
or # for the period 1925-1947, py) = 0-5 £1, My = 0-5 wy, a = 0-01, and f = 0-20. 
4 Equations (2), (3), and (4) then lead directly to the desired sequential tests. 


6. The Negative Binomial Test 


In reference [9] the n.b. distribution was expressed in the form 


T(k + #) p* 
Cie P(e +1) r(k) (1+ pete” (6) 


When #, and #, are substituted in (6) to obtain f, and /, for the likelihood 
ratio, we note that the factor involving the gamma functions cancels in the 
ratio since it does not contain p. Substituting in the 7-th ratio factor of (1) and 
taking logarithms yields 


ss es! : 1+ Po 
a= x; log (5) + (k + x,) log( pe oak 
. For simplification let 
Pr 
b= —+ 7 
op (7) 
and 
ae it Pr 8 
S T+ Po.” : 
Substitution of these values gives 
a 
z,= x, logb + (k + %,) log—. (9) 
‘Summing (9) over the sample (41, %2, --+» %m) yields 
Sus (s7x,) log-— + k mlog—. (10) 


By (3) we set (10) successively equal to logA and then log B which gives 
the acceptance and rejection lines for the sequential test in the form of the 


pair of equations 
m m 
Dua Smt ho, Di Sm+ hy. (11) 


ZAMP IX/4 
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The values of the slope S and the intercepts ty and h, are easily found from the — 
substitution of the two values of equation (3) into equation (10) to be 


loge 
S= Tog (yey @ 
log B 
O log (bJe):? a 
log A 
n= Tog (be) es 


Substituting «= 0-01 and 6 = 0-20 in (4) and taking logarithms gives 
log A = 4-382 and log B = — 1-599. These, together with our assumed values 
1 = 0-3875, k = 12-33, b = 2, and c = 1-1622, give the slope and intercepts 
for the critical lines from equations (12), (13), and (14) and hence the equa- 
tions (11) in the form 


DEF = 3-4195 m — 2-9475 , Doe = 3-4195 m + 8-0774 . (15) 
These are shown in Figure 1 as the two parallel lines. 
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Sequential tests, negative binomial distribution. 
«= 0-01; B = 0-20, 
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To make the test, the accumulation of hail-days by years x; is plotted 
against the year number m. As soon as the accumulation falls outside of the 
band formed by equations (15) the test is terminated. If the accumulation is 
above the band, the evidence from the test is that we should reject the hypo- 
thesis that there was a suppression effect. If the accumulation is below the 
band, the evidence from the test is that suppression had some effect. If the 
accumulation remains within the band, the trials are continued another year. 

As mentioned above, there have been two experiments on the Magadino: 

_ The experiment with rockets which was begun in 1948 and terminated in 1952, 
and the experiment using ground based AgI generators begun in 1953 and 
still continuing. Since the sampling by our network of three stations was 

especially poor for the rocket series, we shall employ this period only in a 
test of heterogeneity. This means that since the stations did not sample a 
1ocket effect, their records were essentially undisturbed by the rocket experi- 

-ment and hence the record should be homogeneous through the rocket test 
period. This is verified by the dotted line of Figure 1 which gives the accumula- 

tion of annual total hail-days through 1948, 1949, and 1950. In 1950, the 
dotted line exceeds the upper critical line indicating that we should. reject the 
hypothesis that the period beginning in 1948 was different from the previous 
record, and therefore accept the hypothesis of homogeneity for the entire 
record including the rocket-trial period. If our stations had been properly 
sampling the effects of the rocket trials and the sequential tests were being 
conducted during the rocket period, it might have been used to make a decision 
about the trials. As it was, this decision had to be made on the basis of other 

- evidence [12]. 

Since data from our three stations were available through 1955, we may, 
in a similar fashion, test whether the period beginning with 1951 is also homo- 
geneous with respect to the historical record. The test is shown by the solid 
line giving the accumulated annual total hail-days for the years 1953, 1954, 
and 1955. In 1955, the accumulation falls below the lower critical line indicat- 
ing that we should reject the hypothesis that the period beginning with 1953 
is a sample from the population which produced the record before that date. 
Since the test was arranged to evaluate decreases in hail frequency, the deci- 
sion given by the test is to accept the hypothesis that the hail frequency was 
lower beginning in 1953. 

As we have seen above, the test has been arranged to readily accept the 
hypothesis that there was a decrease in hail frequency, but at the same time, 
making the acceptance of the hypothesis that there was no decrease 20 times 
more difficult to accept. Thus, even though there was considerable difficulty 
in accepting the hypothesis of no increase or homogeneity for the first period, 
this hypothesis was accepted at average sample size 3 which is near the value 4 


predicted by Figure 4. 


. 
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The test period beginning in 1953 as we saw includes the AgI ground gener- f 
ator trials. Due to the nature of ground generator operations, the sampling 
of the target by our three-station network is somewhat better but was not 
considered adequate according to correspondence with Prof. SANGER. However, © 
there seems to be little doubt that the three-station network does sample the 
suppression effect of the AgI generators at least part of the time. If this is 
true, the results given by the test are quite interesting. They tell us to accept 
the hypothesis that the total hail frequency was reduced somewhat during AgI 
ground-generator operations. This, of course, can be considered only as explora- 
tory information, for the probability of coming to this decision, when actually 
the hail frequency had not been decreased, is high, i.e., 0-20. If we reduce this 
risk of error to 0-10, a more common value used in precipitation augmentation 
evaluations, no decision can be reached with sample size 3 and at least another 
year is required. This lowering of the risk does not, however, change the con- 
clusion reached concerning the first experimental period. 

Clearly the deficiency in the sampling of the network and the high risk of 
error we have chosen for this exploratory test makes it impossible to say any- 
thing conclusive about the experiment. It does seem, however, to show that (1) 
the method of evaluation has some merit, especially if experiments were 
designed for its application, and (2) that the specific application to the AgI 
generator trials gives at least a suspicion of effectiveness which is worthwhile 
following up by further experimentation. 


7. The Poisson Test 


As mentioned above, we give this test for completeness. Because it provides 
a simple hypothesis directly, it is to be preferred in any instance where the 
Poisson distribution gives a good fit to the available data. This distribution 
has the density function 


ev p—-x 
i) ==, (16) 
where yw is the mean value. Taking my and py, for the hypothesized parameter 


values and proceeding in the same manner as with the n.b. we find the loga- 
rithm of the i-th factor in the likelihood ratio to be 


ae lu 
2; = x, log 7 (Ha — Mo) - (17) 
Substituting the analogy of equation (7), i.e., b = fy/@o, and summing over 
t= 152, ..., m gives 


d= (logd) S's, — 2S} mm. (18) 
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As with the n.b. the setting of (18) in turn equal to log A and logB gives 
the acceptance and rejection lines in the form of a pair of equations like (11). 
From these the slope S and the intercepts h, and h, are readily found to be 


/ 


Seal dae fegeorne (19) 
i] 
_ logB 
i ere: (20) 
_ logA 
h, — ‘logb * (21) ' 


These values substituted into (11) give the critical lines. 
Although the 7? for the excess of the variance over the mean was significant 
_ for the 1888-1947 series, it was not for the 1925-1947 series. Tests of goodness 
of fit applied to each series fitted by the n.b. and Poisson distributions indi- 
cated ‘good’ fits for both distributions, i.e., neither of the series departed 
significantly in their frequencies from either distribution. This seems to be 
due to the fact that as k becomes large the n.b. approaches the Poisson and 
since we are dealing with a fairly large k the Poisson distribution fits well. 
On this basis it appeared reasonable to illustrate the application of the Poisson 
test on the same series of total hail-days. 

The procedure is very similar to that for the n.b. Since «and f are the 
same, log A and logB remain unchanged, s4; was chosen to be 4-78 hail-days, 
the value for the 1925-1947 period, and @ is again 2, i.e., Mg is one half of p44. 
Substituting these values in equations (19), (20), and (21) we obtain the slope 
and intercepts which are in turn substituted in (15) to give the critical lines 


m 


Dix; = 3-4483 m — 23070, 3x; = 3-4483 m + 6-3223 . (22) 


The critical lines and test results are shown in Figure 2. As was to be 
expected, the results are quite similar to those for the n.b., hence there is little 
further to add in respect to the interpretation of the experimental results. One 
difference may be noted, however, this is that the critical lines are somewhat 
closer together for the Poisson test indicating greater sensitivity. In contrast 
to the n.b. test, this test, performed with 0-10 risk of accepting the decrease 
hypothesis when it is false, gave decisions identical with those obtained with 
the 0-20 risk test for both periods. The Poisson distribution also has the ad- 
vantage of having no parameter other than the test parameter and therefore 
is to be preferred where applicable. We shall have a better insight into the 
sensitivity of sequential tests after developing the operating characteristic and 
average sample number curves in the next two sections. 
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8. Operating Characteristics 


We are clearly more interested in accepting the hypothesis H, that hail 
suppression is effective than in accepting the hypothesis H, that it is not 
effective. This does not mean that we would like to accept H, entirely falsely, 
and so we have guarded against this by choosing probabilities of acceptance 
and rejection at fy and mu, according to our preassigned estimate of their im- 
portance. This defines the capabilities of the sequential test at these two 
points. To know the capabilities of the test completely, we need to give its 
operating properties at other values of y. In view of our first interest in ac- 
cepting H, an appropriate measure of the operating properties would be the 
probability of accepting H, when the true but unknown hail frequency parameter 
is w. Hence we are interested in a function L(u) which is the probability of 
accepting Hy when w is the true value or power. This is appropriately called 
the operating characteristic (OC), c.f. WALD [8], p. 24. 

From our discussion in the latter part of section 5 it is clear that when 
= fy, L(y) = B; when [6 = My, L(y) =1—a. Since we have chosen 
a= 0-01, L(t) = 0-99 and as u becomes small, L(u) approaches one, or the ability 
to reject H, becomes very high. This is what we need in an exploratory experi- 
ment where we wish to give a phenomenon a reasonable chance to show itself. 
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Apart from L(0) = 1 and L(co) = 0 which are not of great interest here, we 
_have only two points on the OC curve which is not sufficient for our purposes. 
Other points on the OC curve must be obtained by more complex methods. 
WALD [8], pp. 49-50, gives the required analysis which we shall apply to the 
Poisson and n.b. distributions. 
For a discrete variable such as we are dealing with here, WALD, by an 
ingenious analysis, found the basic equations 


= (4,Ha) * Ls 
Xs lien] = 1 ce 
and 
i 
L(u) s- oe a (24) 


such that h/ is a function of w and h + 0. Note that equation (23) expresses the 
fact that the expected value of the factor immediately after the summation is 
unity. These two equations determine L(y). 

We recall that the density function for the Poisson distribution is given by 
equation (16). Substituting this with the parameter values jz) and 1; into (23), 
after some simplification we have 


= (a) 
[exp (2 (uo — oa) — #)] XY = 1. (25) 


Since the summation term is the power series expansion of exp (uw u/mp), we 
find making this substitution in (25) that 
ato 
exp[h (Wo — M1) a] EXP = 1 (26) 


Taking logarithms and solving for gives finally 


2. a el 27 
BS (yl Mo)” — 1 : : ( ) 
In terms of 0 this becomes 
h a Al 
a eg ema (28) 


The OC function may now be drawn from a series of points obtained from 
equations (27) or (28) and (24). The easiest way to do this is to choose a con- 
venient set of values of # and make a table (see Table 2). The last three rows 
have been added to facilitate later computations. The calculation is straight 
forward except for 4 = 0 where equation (27) is indeterminate. For this value 
a simple application of the |’Hospital rule gives the limit. Thus differentiating 
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numerator and denominator of (27) with respect to h gives 


Sa B= Boy Pa Bo 20) 
: a (H1/Mo)*—1 log (Ha / Mo) * sek 
The values from Table 2 were used to draw the dashed curve in Figure 3. 


Table 2 
Poisson OC Function 


h | Ub | L(u) 

=2 6-37 0-041 

ae 4-78 0-200 

2 4-08 0-420 

0 3-45 0-733 

2 2-89 0-935 

1 2-39 0-990 

2 1:59 1-000 

21/2 1-28 1-000 

3 1-02 1-000 
1:00 
0-90 
0-80 
0-70 
0-60 
0-50 
0-40 
0-30 
@ 0-20 

z= 


30 “48 
(Poisson) ——= 


0:20 0:25 0-30 0:35 0-40 
p(negative binomial) ——— 


Figure 3 
Power curves, Poisson and negative binomial distributions. 
== ——— Poisson: (ty = 4:78; Ug = 2:39; w = 0-01; B = 0-20. 
Negative binomial: Py = 0°3875; pp = 0:1938; k = 12-338; x = 0-01; B = 0-20. 


The OC for the n.b. test is found in a similar manner. Substituting the 
density-function equation (6) with the parameter values p; and 4, into (23) 
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after some simplification we. find 


ae ieee AUD et Meee aS Le 60 


i a), Po 1+ py Ik py! ol aed) CG) hs pete 
If this is multiplied by [f9(1 + bo) |**, the result is the negative binomial form 
[Po (1 + Po)l”* {(1 + P) [Po (1+ Pidl* — B [hs (1 + Do) l"}* - (31) 
Raising both sides to the k-th power and solving for # gives 
T+ Po)F 
Sen) eae 
rile st 
es ( (lL + Pe |’ ; oe 
Po (1 + Pr) 
Substitution of b and c for the proper terms gives the simplified form _ 
(1/c)®—1 
P= T= Oe ~ 


Proceeding in the same.manner as with the Poisson for convenient values 
of h, we obtain a similar table from equations (24) and (33). We can save com- 
puting L(p) by choosing the same values of has for Table 2; whence L(p) = L(u) 
for a given / since A and B are the same as before. 


Table 3 
Negative Binomial OC Function 


jh p | Lie) 

Be 0-530 0-041 
ay 0-387 0-200 
~1/2 0-329 0-420 
eis: 0-277 0-733 
1/2 0-232 0-935 

1 0-194 0-990 

2 0-132 1-000 
21/2 0-109 1-000 

3 0-089 1-000 


Again there is a difficulty at 4 = 0 which can be resolved as before by taking 
the limit of equation (33). This gives 
Pee a ae ee 34 
p= lim @e® ~ Tog (je) ° oH) 
Values of Table 3 have been used to draw the OC for the n.b. shown as the 
solid line in Figure 3. 
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The OC curves for the Poisson and the n.b. are remarkably similar and 
close together so there is little reason to differentiate between them. The ~ 
curves show the rather obvious fact that it is much easier to detect big effects — 
of hail suppression than it is small effects. At 2-3 total hail-days for the three | 
stations the ability to detect the effect is nearly perfect whereas near the - 
natural mean for the 1925-47 period yw, the ability is low. Thus if hail sup- 
pression operations could reduce the total number of hail-days to two at the 
three stations in the Magadino, there would seem to be little chance of not 
finding the effect. Of more direct interest in evaluation by sequential testing 
is the length of time an operation would have to run to reach a decision. This — 
comes under a part of the theory concerning average sample number which 
we will now consider briefly. 


9. Average Sample Number 


If a sequential test were applied to a number of sample sequences from the 
same population, the tests would terminate at various sample numbers m. 
These m’s would comprise a random variable and would consequently form 
a frequency distribution and have a mean value. This mean or average sample 
number (ASN), as it is commonly called, would seem to be very helpful since 
it would tell us on the average how long our operations must continue in order 
to come to a decision. 

Wa Lp [8], p. 53, has developed the theory necessary to determine the ASN. 
The ASN as-we might have expected depends not only on Mo and j,, consider- 
ing these again as general parameters, but on the true value wu. The formula 
given by WALD for the expected sample size or ASN is approximate but gives 
values very close to the true ASN. It is 


Eis) og ll Eheivgs 


ee E(e| #) 


(35) 
Here E(m | u) is read the expected value of m for the parameter value yw. All 
other quantities have been met in previous discussions except the expected 
value of z. Since we wish to draw the ASN asa function for all values of Lb, 
there is also a difficulty at h = 0 as with the OC function; equation (35) be- 
comes indeterminate at its maximum value which corresponds to h = 0. Here 
let us assume that w has the value i’. WALD has surmounted this difficulty 
for us by providing an approximate expression for the ASN at its maximum. 
This is 
; — log B log A 
Em | we etapa (36) 


Expressions (35) and (36) require that E(z | fm) and E(z?| uw’) be evaluated. 
This we now proceed to do for the Poisson test. 
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Taking the expected value of (17) and substituting the expression for b 
we have 


b—1 
3 E(z) = E(x) logb — (P35) ly - (37) 
Since E(x) = w, substitution gives ; 
b= 
E(z| p) = plogs — ( 5 ) ma (38) 


Hence by substitution in (35) we have finally the ASN 
Dy logB+[1—L log A 
Pe Sy ee 2 tar Ne (39) 
tw logb — (—} 1 


To evaluate E(z? | wu’) we begin by squaring and taking the expected value 
of equation (17) and substituting b which yields 
b—1 b—1)\2 9» 
E(22) = E(x?) log2b — 2 E(x) my ( a ) logb + ( a =) we. (40) 
For w= mw, E(x) =p’ and E(x?) = w' + w’ by the ordinary power moment 
relationships. Hence 


r i F f b—1 b—1\2 4» 
Bet |p) aw (+e) logtb— 2H! ms (25) tows + (252) at. on 


Since this expression is rather long and is only computed for one value, it is 
easier to find the numerical value of (41) and substitute into (36) rather than 
write one expression of the maximum value of the ASN. 
The values of L(u) of Table 2 may now be employed in equation (39) with 
yw’ = 3-45 in (41) and (36) to obtain Table 4. 
: Table 4 
Poisson ASN Function 


h | L(u) | le E(m| 1) 
wg 0-041 6:37 2-04 
ai 0-200 4-78 3-45 
ANZ, 0-420 4-08 4-27 

0 0-733 3-45 4-23 

1/2 0-935 "2-89 3-10 

i 0-990 2-39 2-10 

2 1-000 1-59 1-24 

21/2 1-000 1-28 1-07 

3 1-000 1:02 0-95 

1-000 0 0-67 


(fee Sse See SE 
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Note that the approximate relationship does not give a true maximum in this 
instance. These values were plotted against the w-scale in Figure 4 and used 


for drawing the lower curve. 


Negative Binomial 
Poisson 


Average Sample Number 


DTK 04 05 0-6 
p(negative binomial) ——— 


H( Poisson) 


Figure 4 
Average sample number curves. 


For the n.b. distribution the ASN function is developed in a similar fashion 
to that for the Poisson. From equation (9) we see after taking expected values 
that 

1 
E(z) = E(x) logb + [k + E(x)] log— : (42) 


Since E(x) = k p we find after substitution 
E(z| p) = kp logd + k (1+ ) log. (43) 
Hence substituting in (35) we find the ASN function 


_, Ep) logB +'t1 — L(y) log 4 
© (m |?) os ~ 7p iogs — OB) loge 


Proceeding to the ASN at 4 =0 we obtain squaring (9) and taking the 
expected value 


E(z?) = log2b E(x®) + 2[k E(x) + E(x?)] log log 
| (45) 


(44) 


+ [B® + 2B E(x) + E(x¥)] log? 
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For the n.b. [9] with parameter p’, E(x) = k p’ and the variance o? = k p’(1— #’). 
From the moment formulae it is clear that E(x?) = o? + E*(x). Using these 
relationships we readily find 


E(x) = hp’ (1+ p' +h’). (46) 


Substituting E(x) and E(x?) in equation (45) would give E(z? | #’). Since this 
would make an awkward expression and is only needed for one value of the 
ASN, it is simpler to compute E(x) and E(x?) separately for p = ~’, substitute 
in equation (45) to find F(z? | #’), and finally substitute this value in (36) to 
find the maximum value of the ASN. 

Using the required values of Table 2 in equation (44) together with p’ = 
0-2773 in (46), (45), and (36) in succession we prepared Table 5. 


Table 5 
Negative Binomial ASN Function 


Pie ee) p | E(m|p) 
==2 0-041 0-530 2-45 
1 0-200 0-387 4.31 
='y0 0-420 0-329 5:39 

0 0-733 0-277. | 5:44 

1/2 0-935 0-232 4-04 

1 0-990 0-194 2:76 

2 1-000 0-132 1-65 
21/2 1-000 0-109 1-42 
3 1-000 0-089 1-27 

1-000 0 0:86 


Using the values from Table 5 on the p-scale of Figure 4, we drew the ASN 
curve for the n.b. test shown as the upper curve. 

The ASN curves are quite different, the Poisson requiring somewhat less 
sample size on the average to come to a decision. Both distributions show rather 
favorable ASN’s, however, and since ASN is number of years on the average 
for our tests to reach a decision, with appreciable suppression of hail frequency 
by seeding this would be detected in about three years. Near the center of the 
interval of indifference as to which hypothesis is accepted the ASN’s have 
maxima. If we consider fy and py or pp and A, as the true values, the ASN’s 
for these values may be read on the lines labelled H, and H,. Had we chosen 
the conventional risks « = 0-05 and f= 0-05 the ASN’s would have been 
2.88 for either hypothesis on the Poisson and 3.64 for either hypothesis on 
the n.b. Note that these are higher for Hy and lower for H, because we have 
made H, harder and H, easier to accept than the test chosen and applied pre- 
viously. If a test with « and f equal 0-05 had been. selected, the decision for 
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the first experimental period would have been unchanged. The test for the : 
second experimental period provided no decision at the third year indicating 
that at least another year of experimental data must be obtained. 


10. Discussion of Results 


As we were careful to point out previously, we consider our results first of 
all a development of methodology, and the application to the Magadino Plain 
operations only an example of the use of this methodology. The development 
we put forth is not in any sense meant to be a final finished product; it is 
simply the beginning of one line of development of evaluation methods. There 
are many ways to develop and adapt these to the individual experiment, some 
of which would certainly involve even the use of other observations and fre- 
quency distributions. If, however, there is one result to be emphasized, it 
would appear to be the futility of using areal controls. One way to compensate 
partly for the lack of a control is to take advantage of the favorable ASN of 
the sequential test. 

On the experimental aspects we must reiterate that much is yet to be 
learned about the physics of hail suppression and to emphasize that our tests 
prove nothing conclusively about the suppression effects of AgI ground gener- 
ators. A phenomenon which was suspected might exist was given a preassigned 
favorable chance to exhibit itself and did, in fact, show itself in the results of 
our tests. This would appear to be exploratory information of some value in 
planning future experimentation. 

Further consideration we feel should be directed toward both the theoretical 
and experimental aspects of artificial suppression. At present the theoretical 
basis of hail suppression by artificial nuclei as well as qualitative experimental 
evidence obtained by eye and radar seems to be much weaker than that on the 
stimulation of rainfall at cloud level. It would appear therefore that the results 
presented above as well as other clues of recent years on hail suppression now 
more than ever need confirmation or refutation by theory and case study by 
radar and other means. It would appear that this project should be continued 
so that experimental results might be tested under more stringent probability 
levels on the hypothesis that suppression is effective. Concurrently, the physical 
model, on which such testing is based, should be strengthened wherever pos- 
sible by design and additional observation. 

It should be emphasized that experiment IT, for which the above analysis 
has been made, was not a statistically designed experiment, and therefore 
our analysis had to be designed for data which had already been collected 
according to procedures beyond our control. In this sense our analysis had to 
be planned in much the same way as for commercial operations. It seems un- 
necessary for a statistician to have to say that data from an experiment, for 
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which the observational program and analysis have been designed in advance, 

would be much preferred to the data which were available. We understand that 

_the operators of this experiment have also recognized this and have now de- 
‘signed a new experiment. 
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Zusammenfassung 


ee 


Ein Sequenztest-Verfahren wurde fiir die Priifung der Annahme entwickelt, - 
dass Hagelverhinderungsversuche die Frequenz der Hageltage verkleinern. Das 
Verfahren ist auf die Angaben itiber den Grossversuch II zur Unterdriickung — 
des Hagels in der Magadinoebene (Schweiz) angewandt worden. Funktionen fiir 
die Versuchskenngréssen und Stichproben-Mittelwerte sind ebenfalls entwickelt 
und am selben Beispiel zur Anwendung gebracht worden. Schliesslich werden die 
Folgerungen, die sich hinsichtlich der Versuchsverfahren und Testergebnisse 
abzeichnen, er6rtert. 


(Received: July 20, 1957.) 
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An Asymptotic Method for Free Convection Past a Vertical Plate’) 
By Jor R. Foote, West Lafayette, Indiana, USA?) 


In a recent note FinstTon [1]%) observed that the transformation given by 
E. POHLHAUSEN, used to reduce the free convection equations to ordinary dif- 
ferential equations, belongs to a family of similarity transformations which allow 
the plate temperature to be a function of ¥ rather than a constant: 


O9(*) = 1+ a x h(0). (1) 


Besides permitting this function, the generalization may also permit bringing /’, 
hence u, into closer agreement with experiments made on large plates, without 
much change in agreement with 0,. Conceivably w would be chosen differently 
for each plate size. 

It is the purpose of this note to give two terms of the solution of the equations 
in addition to the obvious first terms of the asymptotic expansions. This solution 
will hold for arbitrary w, o, «, where « is used for the B of [1]. Details of this 
theoretical method were given by Ferris [2], who gave as a concluding remark 
the first term, beyond the obvious first term, for the usual case wm = 0. 

The streamfunction and temperature ratio are respectively : 


p(x, y) =A xm f() PSI (2) 


IRGS 
O(x, 9) = SP) a + a 9 hy) (3) 
foe} 
When these functions are used in the convection equations, it is found that A is 
inessential and that equation (1) can be satisfied if « is left arbitrary and h(0) = 1. 


1) Written under the sponsorship of the National Science Foundation, Contract NSF-G3050 
New York University. Reproduction in whole or in part permitted for any purpose of the United 
States Government. 

2) Purdue University, Division of Engineering Sciences. 

3) Numbers in brackets refer to References, page 67. 
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Also 3 = 
4 m= @ Sie : n= ae (4) 


-and with o as the Prandtl number the equations with boundary conditions are 


"4 


"+ mff"—(m+n)f?+ah=0, (5) 
h"+o(mfh’—-ofh)=0, (6) 
y=O f=f=0, h=1; n>oo: f'>0, h>d. (7) 


It should be noted that to reach the equations of POHLHAUSEN we would use 
foi O ct Den 4 pee AL Te 

The method uses a series for f(7), h(n), beginning with f and 0 respectively, 
since we assume 


a ee One. 


and the boundary conditions give 


Lime (97) 40 


™4—>0o 


Then assume a series in powers of some parameter p: 


Hn) = B+ Palm) + P* fal) +o, h(n) = 0+ Phan) + DP aly) Ho (8) 


Comparing equation (8) to the boundary conditions (7), it is found that it is 
sufficient to take as boundary conditions: 


f4(0) per, MOF MO=% 722, 


fy(0o) = 9, hy(oo) = 9 O0= B+ ph (0) + p? f.(0) +:--- 


The last equation of (9) is /(0) = 0, and will be used to determine f. Thus, with B 
and 0 to care for the conditions at oo, and f,(7), 4,(7) to care for the remaining 
conditions, the terms f,(7), 4,(), # = 2, serve to bring the expansions (8) into 
agreement with the solution. 
Substituting equation (8) into equations (5) and (6), and collecting terms by 


powers of p: 


0, f,(co) = 0 


fa tmBoh=0, | | es 
“lil +m Bh s—aebr, 
4 hi tm Bohs =o(ofta— mfr h) | ts 
‘ Sure tt = tm + 7 vo—mpfi — «he, 


3. ht +m Bohs =o[o (fi he + feta) — ™ (fifa + fe ei) (12) 
‘ CAP aa Erna) 1) — «hg - 


Clearly the equations for /,(m), y(n) could also be written, and each pair of 
equations has the same homogeneous part as the first pair. Since the equations 
have constant coefficients, it is merely a matter of algebra to carry k to as lorge 
an integer as we please. It is‘found that f, and h, each contain the factor 1 [p*; 
hence, we take p = 1 in equation (8). For o + 2-", y any non-negative integer, 


ZAMP IX/5 
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the first three pairs of functions are 
ieee ei), 
fo(n) = 6, 
SS eels DP 
oe ere S 


(m B)* 0? (1 — o) 


—mBon — —me in 3n __ 4—mBo 
h(n) a Be (1 aa a) >; +o (1 a q: Bn) = 262 (1 é ")| ; 


2027 (w—mo)— 3n(1+o0) —mBn —mBon 
fo(n) = ot | a B)? o we 0)? (1 + 0) (oe —e ) 
o3 (1 — o) (w— mo) — 2na (1+) +m (1 — oa)? (1 4+ 0) 
ks (m By? of (1 — 0)? (1 + 0) 
CPS ot ee) ee 


ees —mBn _ ,—2mBo 
+ 8 (m B)? «4 (26 — 1) (1 — o)? (20e e 5 n) 
n sara baby 
7 “Aim phot (1 — a)? (2 — e-™8n) ¢ "|. 


a 


By using these functions in equation (8) it may be seen that the proper form _ 


of the series expansion is 


An) = B+ a Fy(n) + a? F(m) +--+, (mn) = Ayn) + & He(m) +--- 
There is no difficulty in using a limiting process on the above functions if the 
solution for o = 1 is desired. 
To determine the constant £ from /(0) = 0, we have an eighth-degree biquad- 
ratic equation: 
pe — a B4 i ce (o — mo) — 2no(1+o) + m(1—o)? (1 +0) 
m’? o§ (1 — o)? (1 + o)8 
207 (w— mo) — 3n(1 +o) n(29° 3 3e= 1) 6 
2m? o4 (1 — o%) art aa ; 


To check this equation we may put w= 0, a=1, m= 3/4, n= —1/4, and 


B= 44, knowing that PontHauseNn gave x = 0-6 (see, e. g., [3], p. 449) when 
o = 0-733 for air: 


ee oe bs 4 3 (1 — a)? (1 + 0) +2o(1+4 0) — 3 afl — o) 
21x 64.3 16 ¢? "© 48 X 27% 63 (1 — o)® (1 + 0/8 i | 
tas 
Ate = 26% 1—30— 20? _— fe 
82% 272 o%(1— of) * 24xi6xa7of—aje 


The appropriate root turns out to be x = 0-612. With only two terms of the 
Series FEettis found + = 0-513, so that the next term has given significant 
improvement. With patience one could easily solve also equation (12), and likely 
would find that four terms yield all that experimental data would require. As 


seen in equation (13), even for 7 = 0 the higher-order terms are quite small, the 
constant term being about 0-2 x 10-8 for air 


-_ 
- 
' 


eo, 


f 
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Since a proof that this method converges is lacking, it is worth noting that the 
solutions (9) become power series in «, so that convergence is assured for suffi- 
ciently small «. However, an iterative numerical procedure, using perturbations 
from f, and h,, can be based on the prediction of three terms, and that method 
‘might well be the more economical one. 
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Zusammenfassung 


: Das Problem der freien Str6mung an einer ungleichmassig erhitzten Wand wird 
auf die Berechnung einer Reihe von Funktionen zuriickgefiihrt, von denen jede 
einer linearen Differentialgleichung mit konstanten Koeffizienten geniigt. 


(Received: March 5, 1957.) 


Remarks on the Paper by M. Finston: 
‘Free Convection Past a Vertical Plate’ 


By FRANK NIUMAN and KARL POHLHAUSEN, Dayton, Ohio, USA?) 


Recently M. Fryston has shown in a ‘Brief Report’ [ZAMP 7, 527 (1956)] 
that under the assumptions of the boundary layer theory the problem of free 
convection past a vertical plate has an exact solution for a plate temperature 
which is proportional to a power of the distance from the leading edge of the plate~ 
The well-known solution of a constant plate temperature is included as a special 
case. 

By a similarity hypothesis Frnston derives the following system of differential 
equations for the temperature and velocity distribution (the notations are those 
of FINSTON): 


o+1 
2 


t wo +3 " m “ w +3 “ 5 1 " 
ae ott = IN + Bh, of h Gol om. ls 


where primes denote the differentiation with respect to the argument 7 of the 


functions f and h 
n= yxe-nit 


The boundary conditions are 
for 7=0: f=f=9, eles 
for n= co: f=h=0. 


1) Aeronautical Research Laboratory, Wright-Patterson Air Force Base, USA. 
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Thermal and dynamical profiles and stream function for Pr = 0-733 and for parametric values of w. 


Figure 1b 
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For a given Prandtl number the solutions form a two-parameter family with 


and B as parameters. A ome-parametric problem is obtained by the following 
substitutions: 


H(n) = BI4 FH), h(n) = H(m) with 7= By. 
One then finds 
; wo +3 1 


aoa 2 o+3 
F’ 7 BW PL LD H / ie ee 
2 4 Seah se Tah 3 FH aes 


Here the differentiation is to be carried out with respet to 7. These equations 
were evaluated numerically for air (Pv = 0-733) using the Bendix Digital Dif- 
ferential Analyzer. The figure shows the quantities H, F and F’ as functions of 7 
for m = 0, 0:5, 1 and 2. The values of H’ and F” at the plate, which determine 
the heat transfer and the shearing force as well as the asymptotic values of F 
are given in the following table. 


eee eS 


H’(0) |—0-359 |— 0-386 |—0-410 |— 0.450 |— 0-483 |—0-512 |— 0-536 |— 0-630 
F’ (0) 0-953 | 0-926} 0-902] 0-861 0-828 | 0-800] 0-777] 0-649 
F (oo) 1-684 | 1-610] 1-550 1-412) 1:327 Ihe 1-169) | 0-957 


For w = 0-2 the heat transfer along the surface of the plate is constant. The case 
may be of some technical interest. 


Zusammenfassung 


Die von M. Finston abgeleitete exakte Lésung fiir die nattirliche Konvektion 
an einer vertikalen Platte, deren Oberflachentemperatur nach einem Potenz- 
gesetz veranderlich ist, wird fiir gegebene Prandtl-Zahl auf eine einparametrische 
Schar zuriickgefiihrt. Einige Lésungskurven fiir Luft (Pr = 0,733) werden mit- 
geteilt. 


(Received: July 6, 1957.) 


An Inequality Relation Between Yule’s Characteristic K and 
Shannon’s Entropy H 


By Gustav Herpan, Bristol, England?) 


If the entropy is used in linguistics in the form and with the interpretation 
given to it by the communication engineers [1]*), as so and so many ‘bits of 
information’, it requires much additional explanation, and even then confusion 
is not always avoided. For one thing, it may be now regarded as established that 
the entropy H measures not the information, but rather the lack of it, and that, 


1) University of Bristol. F 
2) Numbers in brackets refer to References, page 73. 


70 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 


therefore, the negative entropy — the negentropy — should be used for measuring 
information [2]; secondly H is a logarithmic quantity, and as such not readily 
interpretable in non-logarithmic terms; more especially the ‘bits of information’ | 
have not really a linguistic meaning. : : 

The transformation which we shall propose takes these shortcomings of the 
traditional expression of entropy into account, and by doing so changes it into an © 
easily interpretable statistical quantity of general applicability and of interest to — 
linguists — the repeat vate. It is based upon an inequality relation between the two © 
parameters which are indigenous to literary statistics: SHANNON’s entropy 
H = — J)’ p; log,p;, where p; is the probability of occurence of the symbol i, and 
YULE’s characteristic (apart from a factor 104 to avoid decimals) 


Sa ‘ 


where 


N =D %,. and)S., = Sir? Ny, (1) 


and where u, is the number of different words occuring with the frequency r¢ [3, 4]. 


The Transformation 


The negative entropy —H = Y p; log.p; is the arithmetic mean of the loga- 
rithmic variable log,p. Its antilog is, therefore, the geometric mean G of p;- In 
symbols: 

Antilog (-H) = G. (2) 


The arithmetic mean A of p, as a weighted average, or the expectation of p, is 


A=)"; 2; - (3) 


Dy b; is the repeat rate, i. e., the probability that two particular linguistic events, 


say words, selected at random from the text under consideration should turn out 
to be the same word. 


Using the inequality between arithmetic and geometric mean 
GS Ao, (4) 


we obtain the following inequality between the entropy and the repeat rate of 
particular events: 


Antilog (—H) = G 5 A = }"p?. (5) 


On the other hand, it can be shown [5] that YULE’s characteristic K is an 
unbiased estimate of the sum of the repeat rates for particular linguistic events 


Oo ns 


Proof 


We denote ie Cm 4 measure of heterogeneity of the statistical population 
ions that c,,= Pi + py +++ +p”, where p, 1s the probability of occurence of 
the w-th word. Then the expectation of that parameter for which we will write 


ul 


4 
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- ¢ is given by: 


= 


4 1 1 . “1 (N 
ae [A = r a 
fm Nil "3 Apom Nim a. de 2, () Py (? a : (6) 
¥. = 
: 1 N 
= [m] (ay a Nyries N-r 
a2 nin” () Pat ey) 


_ : CN m)! r—m —m—(r—m) ».m 
2 2 em mom te Pe 


wal 


In other words @,,, is an unbiased estimate of c,,,. 
For our case m = 2, and we write 


1 bee 
eEceE ie ay Cia eet: (8) 


The Inequality 


Yutr’s characteristic K can be shown to equal ¢, (1 — 1/N): 


esi aN te ED (9) 
and thus tends to p? + $3 +--- as N oo. It appears thus to be the repeat 
rate of the linguistic events in question. 

We write therefore: 


, (1-3) S,—N ( S,—WN 


Antilog (-H) SK. (10) 

In words: 
The vepeat vate ov the characteristic is at least equal to the antilog of the negen- 
tropy, or the antilog of the negentropy is a lower limit of the repeat vate or of Yule’s K. 
The relation between the interpretation of K as the repeat rate and as the 


coefficient of variation of the mean, given in a former paper [4], is as follows: 
In formula (1) both and , are frequencies, of word occurences and vocabu- 


lary respectively, and may by division through the corresponding totals be 
transformed into relative frequencies. For great samples, these relative frequen- 
cies are sensibly stable and may be regarded as estimates of the corresponding 
probabilities, p, = 7/N and a, = Np|n (N= Ny + Ng + ++ Ns). 

Using the first transformation, we have 


S 2 v2 - 2 
Kt = 2 =D) age DP Or 
r=1 r=1 


and summing over the individual p, instead of over their groups ”,, we may 
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K* - YP, 
p=l1 


which shows K to be the repeat rate of particular words?). 
For using the second transformation, we write 
Ke it Sis 1 nm Rv? ny — (Xv m,)* 


[Sine ew Ec n( Xv n,|* 


write: 


Dividing numerator and denominator by m3, the expression on the right becomes 


{Yrialn—(Zra]2}n _ 9,/n 
[Sv a} Wp 


which shows K, for great 1, to be equal to the coefficient of variation of the mean 
of the variable 7. 

Relation (10) is of great importance both practically and theoretically. From 
the practical angle, it enables us to obtain an estimate in terms of a lower limit 
of the repeat rate if K cannot be calculated because the variable is a qualitative 
one. Furthermore since, as a rule, H is easier calculated than K, it provides a 
method for obtaining quickly a lower limit for the latter. From the theoretical 
angle, it provides the sought-for correct linguistic interpretation of the entropy 
in terms of the repeat rate of linguistic forms.-« 


Application (Uni-Variate) 


Using the inequality (10), we can express, for instance, the difference in the 
entropy for length of English words in terms of number of letters H = 2-628 and 
that for word length in terms of sounds (phonemes) H = 2-274 [6] by saying that 
the repeat rate for specified word lengths in English is at least 1 in 6-3 words, if 
word length is measured in terms of letter number, and 1 in 4:8 words, if word 
length is measured in phoneme number. In other words, picking two words:at 
random from a very long text and examining them for word length, we shall find 


that word length repeats itself slightly more often if measured in terms of pho- 
nemes. 


Application (Bi-Variate) 


Language translation may be regarded as bi-variate coding of messages [6]. 

If the entropy is calculated for a given quantitative characteristic, such as word 
length‘), for each of the two languages, original and translation, independently, 
the values so obtained will differ, as a rule, among themselves and also from the 
value of the bi-variate entropy, that is the entropy calculated from the word 
length of translational equivalents (words denoting the same concept in the two 
languages). More precisely, the conditional entropy, i. e., the entropy for word 
length in translational equivalents of words in the original of a specified length is 
less than the unconditional entropy. Thus for translation from French into 
English we find the former as 0-923 bits and the latter as 1-410; in translation 
from German into English the former is 1-162 bits and the latter 1-648; in 


3) K* approximates to K as N->oo. 
4) In terms of syllabie number. 


? 
YY 


Vol. IX,1958 — Kurze Mitteilungen — Brief Reports - Communications bréves 7s} 


‘translation from Russian into English the former is 1-153 bits and the latter 
2-100 [6]. 

These relations are not too clear as long as they are expressed by means of 
the entropy in its conventional form. The inequality relation (10) enables us to 
“make the difference between the conditional and the unconditional entropy more 
‘readily understandable. In terms of that relation we express the differences as 
follows: for French-English translation the repeat rate of words in the French 
original is at least 1 in 2-6, or 1 in 3 approximately, and rises in the translational 
equivalents. to at least 1 in 2; for German-English translation the repeat rate of 
word length in the original is also 1 in 3 approximately, and in the translational 
equivalents rises again to at least 1 in 2 approximately; for Russian-English 
translation the repeat rate of word length in the original is 1 in 4 approximately, 
and rises in the translational equivalents to 1 in 2 (2-4 precisely). 
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Zusammenfassung 


Nach einer einleitenden Erérterung der Unzulanglichkeit der bisherigen 
Interpretation von SHANNONS Entropie H fiir Zwecke der strukturellen Linguistik 
wird gezeigt, dass zwischen SHANNoNs Entropie H und YutEs Charakteristik Kk, 
als einem Ausdruck fiir die Wiederholungshaufigkeit p? von Sprachformen, eine 

-Ungleichheitsrelation wie folgt besteht: 


Antilog (-H) < K = J ??, 


die sich als eine Beziehung von ziemlichem praktischem und theoretischem Inter- 
esse erweist. 

Was die praktische Seite betrifft, so erméglicht sie bei Kenntnis von H die 
schnelle Berechnung einer unteren Grenze von kK, wenn diese Grésse nicht direkt 
berechnet werden kann oder nur mit grésserem Arbeitsaufwand. 

In theoretischer Hinsicht erweist sich die Relation als die lange gesuchte 
linguistisch relevante Formulierung der Entropie, im Sinne der Anwendung 
dieses Ausdruckes in der Informationstheorie, als eine untere Grenze der Wieder- 
holungshaufigkeit von Sprachformen. » 

Die Anwendung der Ungleichheitsrelation, und die Vorteile der implizierten 
Auslegung von H gegeniiber der friiheren Auslegung dieser Grésse, werden an- 


hand von Beispielen gezeigt. 


(Received: July 25, 1957.) 
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A New Approach to the Theory of an Oscillating Cylinder Viscometer_ 
by Dimensional Analysis 


By Att ABpEL KeErim IBRAHIM, Alexandria, Egypt?) 


The Theory 


The experimental variables which influence the deflection 6, of the inner 
cylinder are listed below: 


Table 1 
Name of Quantity Symbol Dimensional oe 
Moment of inertia of inner cylinder i ML? 
Angular frequency of outer cylinder o lige 
Torsion constant of wire G MEP ae 
Viscosity . n iM Ea 


Other variables such as the radius of the inner cylinder, the radius of the 
outer cylinder and the length of inner cylinder immersed in the liquid may also 
appear in the final solution, but since the dimensional formula of each is L, 
therefore we shall consider an apparatus constant K which include them and has 
the dimensional formula L’, where / is an unknown constant number. 

Now since J w? and t have the same dimensions and each is acting opposite 
the other, therefore (I w? — t) may be given as a single variable. 

The dimensional formula of the several variables in terms of the unit chosen 
as fundamental (namely, mass, length, and time) are given in Table 1. The most 
general form of the product of powers of the variables is 


a 
(<*) RK? (I w? — 1) we xf , 
Po 


where a, b, c, d, e are unknown constant numbers and Po is the deflection of the 
outer cylinder. 
Substituting the dimensional formula for each of the variables into the 


previous expression and equating to zero the indices of the three fundamental 
units, we have the following three equations: 


Cpe ome, Fibre! os Cimn Oy aes Cio wie al a) < (1) 


Solving equations (1) for d, e and b in terms of the other index c, we obtain 
for the general dimensionless product 


0)\% (1 wo? —7e 
Po ( Kit ws 
Since the combination of variables within each bracket is raised to an arbitrary 
power, each such combination must form a subsidiary dimensionless product. The 


1) Physics Department, Faculty of Science, The University. 
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general solution of the problem may then be written 
M\. (lo=r 
(3°) - (  @ Rar) (2) 
Erhere gy is an unknown function of the dimensionless product whose form can 


be determined either from theoretical or experimental [1, 2]?) data. 


For a particular apparatus K is constant, we may, therefore, conclude that 
under such circumstances: 


ise () 


The Single General Curve 


It therefore follows from equation (3) that a plot of (g,/0)) as Y against 
[(I w? — t)/n wm) as X gives the single general curve. This curve does not depend 
on either the viscosity or the density of the liquid but does depend on the dimen- 
sions of the apparatus used. 

Figure 2 represents the single general curve for our apparatus and Figure 4 
that of OLDROYD’s apparatus. This curve is based on either our theory [2] or an 
experimental data of a liquid of known viscosity. Also Figures 1 and 3 represent 
the variation of [(I m2 — t)/m] versus n (n being the frequency where w = 2 an) 

for our apparatus and OLDROYD’s apparatus respectively. 


Viscosity Determination 


The value of (q,/0,) for a certain liquid is measured for a certain frequency n. 

Then it follows that: 

(1) From the (x [(I m? — t)/@]) curve we get the value of [(J w*? — t)/@] corre- 
sponding to the value of used. 

(2) From the {[(I w? — t)/n J, (o/9)} curve we get the value of (I w? — t)/7 wo] 
corresponding to the value of (9/09) measured. 

(3) From dividing the previously obtained values in step (1) and step (2) the 
value of 7 is determined. 


Experimental Measurements 


Experiments have been carried out on parraffin oil, glycerine and olive oil. 
The work was carried out at room temperature, which remained constant during 
the half hour required to take any one set of observations. 

Tables 2 and 3 below give full details of the readings taken and results obtained. 

This offers a consistent and convenient method for measuring 7. 


Table 2 
efi Temperature n Po 
Liquid oC é Js 6, 
ee es nS 
atciiimrOUMeane tN al «= so sass 20 aS 2-085 
Gy CCTGs knee fa pore us sp 20:5 3-0 1-746 
Oliverol 7. ! - Ys : 0-85 1-162 
Paraffin oil [3] (Gusev s sresults) 


2) Numbers in brackets refer to References, page 78. 
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SKe 


Figure 1 


=.) curve for our apparatus. 
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Figure 2 


Single general curve for our apparatus. 
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Figure 3 


Iw? —Tt t 
n— —J curve for OLDROYD’s apparatus. 


Figure 4 
Single general curve for OLDROYD'S apparatus. 
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Table 3 

att. Io?—t Ia?—t n 

Liquid ( oe =) ae ; Heike 
Paratfinvoll "2s. wee 1-8 0-94 1-91 
Glycerine”. -. 9.5 3 ao eg ee 4:37 0-75 5-83 
OUVe Ol eh a. fe ree 0-3 0-3 1-0 | 
Paraffin oil [3] (OLDROYD’s results) 0-049 0-0295 1-662 | 
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Zusammenfassung 


Es wurde ein Verfahren entwickelt, um den Viskositaétskoeffizienten 7 einer 
Newtonschen Fliissigkeit mit Hilfe eines schwingenden Zylinderviskosimeters zu 
messen. Es wird gezeigt, dass die Ablesungen von (6)/g,) fiir verschiedene Fliis- 
sigkeiten zu einer einzigen Kurve fiihren, wenn sie als dimensionsloses Produkt 
der verschiedenen Variablen aufgetragen werden. Diese einfache Kurve hangt 
nicht von der Viskositat und Dichte der Fliissigkeit ab, hingegen von den Dimen- 
sionen der Apparatur. 

Experimente wurden an Paraffinél, Glyzerin und Olivenél ausgefiihrt. Die 
Ergebnisse scheinen sowohl die Theorie wie das Verfahren zu stiitzen. 


(Received: March 28, 1957.) 


A Note on the Second Fundamental Problem of Elasticity Applied 
to a Plane Circular Ring 


By FRANTISEK Kroupa, Prague!) 


A brief report by W. A. Gross [1]2) from 1957 gives the general solution of 
the second boundary problem for annular regions by the method of MusKHELIsH- 
VILI and as an example treats the mutual displacement of both circles. An inter- 
esting new point in this contribution is that attention is drawn to an error in the 
older paper by REISSNER. 

I should like to point out that a detailed general solution of the second boun- 
dary problem for annular regions, together with specialization for an identical 
case of mutual displacement of the circles, also given by using expansion in 
series in a complex variable by the method of MUSKHELISHVILI, was published in 
1952 [2]. The same results are obtained here as in paper [1], but the procedure in 

1) Institute of Physics, Czechoslovakian Academy of Sciences. 

2) Numbers in brackets refer to References, page 79. 
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calculation differs in details. Specialization for the mutual displacement of circles 
is also discussed in this paper with respect to technical applications and it is 
shown that this solution can be used as the first approximation in calculating a 
rubber damper in the shape of a hollow cylinder. 

' Another general solution of this problem, particularly suitable for numerical 
calculation, is given in the monograph [3] published in 1955, devoted to making 
more accurate and further developing the method of MUSKHELISHVILI. The solution 
here, on pages 308 to 322, is given as a combination of the Muskhelishvili method 
of the function of a complex variable and the variational method, and here too is 
re ‘detailed calculation of the same special case of mutual displacement of the 
circles. . = 
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(Received: July 29, 1957.) 


Comment on the Preceding Note 


The author regrets that he was unaware of the publication by F. Kroupa, 
The aim of the present paper was to correct the classic solution of H. REISSNER. 

It is of passing interest that the author was able to solve a plane circular 
ring problem with mixed boundary conditions of unknown extent. Solutions 
for the first, second, and mixed boundary value problems were programmed 
for digital computer evaluation using fifteen Fourier coefficients. The final 
result was obtained by successive improvements. W. A. GRoss 


(Received: November 28, 1957.) 


Nachtrag zu «Ein Verfahren zur Stabilitatsfrage 
: bei Matrizen-Eigenwertproblemen» 


Von Hans-Rupo.ir Scuwarz, Ziirich’) 


Nach Erscheinen meiner Veréffentlichung [1]*) haben wir erfahren, dass das 
yon mir dort auf Seite 474 zitierte Kettenbruchkriterium zuerst von H. S. WALL 
nur fiir Polynome mit reellen Koeffizienten aufgestellt und bewiesen wurde Wale 
wiahrend die Verallgemeinerung auf Polynome mit komplexen Koeffizienten auf 
E. Franx zuriickgeht [3, 4]. Fiir diesen wertvollen Hinweis von E. Frank danken 


wir bestens. 


1) Institut fiir angewandte Mathematik der ETH. wa 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 80. 
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The Flow of a Non-Newtonian Liquid Near a Stagnation Point 


By ABINASH CHANDRA SRIVASTAVA, Kharagpur, India?) 


1. Introduction 


FROEsSSLING?) has obtained an exact solution for the flow of a viscous liquid © 
near a stagnation point. A stagnation point occurs when a liquid stream impinges 
on a wall at right angles to it and flows away radially in all directions. 

In this paper equations of motion for the flow of non-Newtonian liquid near 
a stagnation point have been set up in cylindrical coordinates and integrated by 
Karman-Pohlhausen method. It is observed that the cross-viscous effect depends 
on a non-dimensional parameter K = pu, a/m (where B, Be are coefficient off 
viscosity and. cross- -viscosity respectively and a is a constant of dimension 7-1, | 
depending on the velocity in potential motion). The boundary layer thickness 
decreases when K varies, from K = 0 to K ~ 0-15 and then increases with K. 
The shearing at the wall decreases with the increase of K, while the normal stress 
at the wall is independent of both viscosity and cross-viscosity. 


2. Equations of Motion 


The equations of steady axisymmetric motion in cylindrical coordinates, 
when the velocity component perpendicular to the meridian plane is zero, are 


Ou Ou OPry OP,» Pry im Poo 
@ (wu pe Re )= greet age ¢ ; oa | 


ow Ow Op, OP. , Pye ; 
a(n sp +e 7 aia ear (2) 


where wu and v are velocities in directions of y and z respectively. The equation | 
of continuity is 


“Or y Oz i (3) 


» Matianages Department, Indian Institute of Technology. 
2) H. Scuricurine, Boundary Layer Theory (Pergamon Press Ltd., London 1955), p. 73-75. 


. 
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_ The stress matrix || p;; || and the rate of strain matrix || e, ; || in the case of 
non-Newtonian liquids are related as) 


ine Pre Pox 100 Cnn org Ore Cnr v6 ng Cry 20 ve 


Por Poo Poz || = —P || 919 |) + wI| 6, 290 202 || + Me || Cor 00 Coz Cor "66 “bz 


|| Pee Pr9 Pes 001 Cer org Oxe Cur er exe Cur e270 exz 
where 

Ou u Ou ow Ow 
oe, 475 Dy oct Barra. 2,9 = 9, = 0, Cag FE a Mee ne tar vaass (5) 


A stream of non-Newtonian liquid impinges on a wall z = 0, and flows away 
radially in all directions. The stagnation point is at the origin and the flow is 
in the direction of the negative z-axis. We shall denote the radial and axial com- 
ponents of the velocity in frictionless flow by U and W respectively, whereas 
those in viscous flow will be denoted by u = u(r, z) and w = w(r, z). Theboundary 
conditions are uw = 0, w= Oatz= Oandu= U at z= 6*, where o* is the boun- 
dary layer thickness. 

For frictionless case we can write, U = ar, W = — 2.42, where a isa constant. 
It is seen that such a solution satisfies the equation of continuity. Denoting the 
stagnation pressure by », we find the pressure in ideal flow, 


p— y= © (U4 W2) = Far? + 428), 


The frictionless case suggests that we should take the following form of the 
solution for the velocity in the boundary layer, 
u=rf(z), w=—2f(2), (6) 


the solution of this form satisfies the equation of continuity identically. Substi- 
tuting from (6) into (1) and (2), we get 


vy (f/? — 2ff" — y fe + a = -—.$ +r, f"?, ; (7) 
0 
Aff 4:29 f"— 28», ft” — Ly fee = -— ee. (8) 


The condition of integrability is 
8 (2¢f'=f24 vf" — aS A es + vy, f"?) = ERO TE we 

Integrating, 
g Sfpe af of 2a =H f+ So, (9) 


where C is an integration constant. In view of the boundary conditions 
(or) =a, f"(68) = f"(O*) = 0, C= —-a?. 
1 op 


meee rs: 


From (7) and (9), 


_8) I. Braun and M. REINER, Quart. J. Math. appl. Math. 5, 43 (1952). 
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(4) 
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Integrating, | 

p= ~~ a? [v2 + 4F(z)])+m,7? f", (10) 

where F(z) is an arbitrary function of z, to be determined from (8). Substituting 
for p in (7) and (8), we have 

PRE | pm Ome Rf ON Se op eds (11) 

2) He of ta ve pe. (12) 

The boundary conditions are | 

iar Oman Cl Or aie cia Ol i a ee 


The constants a? and vy can be removed from (11) by a similarity transforma- 
tion 7 = ie i ath and f(z) = /a v y(n, 6), where 6 = |/a/v 6*. The non-dimensional 
form of (11 


6p” + 209 yp" — 6 gp + 63] = K (p"2? $2 g’ o"), (13) 


where a prime denotes differentiation with respect to 7. Now, boundary condi- 
tions are 


(0) = y(0) = 0, g’(1) = 6, p"(1) = 0 and 4d[y”(0) + 63] = K y"%(0). (14) 


3. Solution of Equations by Karman-Pohlhausen Method 


Assume 
(1) = a + ag + ag? + bP + c Mf. (15) 
Applying. the boundary conditions (14), we get 
6+ 36 j 
e(n) = Aen Dt by eiaegawt, (16) 
and 
9K b*—6bd6(3K+4)+ 62(9 K — 462) =0. (17) 
Integrating equation (13) with respect to » between the limits 0 to 1, we get 
1 1 
" \ / é, 9 « K 
[p"]o + 2 Ole w’]s — 3 a p*() dn + 6% + 5 fe) dn=0. (18) 
0 


Peat the value of (7) from (16) in equation (18) and putting 62 = 
b/d = p’, (11) becomes 


25 eee | 


9 39 i 
(ove a) pe (U “+ SK +15) p’— a e+ 4Gg K+3=0 (19) 
and equation (17) reduces to 
9K p?— (18K + 24) p'+9K—4¥=0, (20) 


Elvneeee p’ between (19) and (20), we have 


9 v4 + (132 — 624 K) *8 + (10609 K® — 3916 K — 2736) x? 
+ (26523 K* + 15634-5 K® + 32235 K + 7560) « + 15876 K4 Sat 
+ 40131 K8 4+ 54684 K® — 101430.K.0, 
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Putting K = 0 in (21), we have the solution corresponding to the flow of a 
viscous liquid. In the reduced equation, either 7 = 0, which is not possible or 


3 «84 44 x2 912 ¥ 4 2520=0. (22) 


' Equation (22) possesses two positive and one negative roots. We need not 
consider the negative root, since 6 cannot be imaginary. Two positive roots are 
3-491 151 and 8-894 263. After plotting the velocity profile corresponding to 
these values of 6?, it is found that corresponding to 62 = 8-894 263, u/U becomes 
greater than 1 with in the boundary layer region which is absurd. The velocity 

profile corresponding to 6% = 3-491 151 agrees closely with that given by FRoEss- 
LING?). For K = 0-02, 0-04, 0-06, ... , 0-30, we find that root of the equation (21) 
which is nearest to the smaller root of (22), namely, 62 = 3-491 151. The shearing 
stress at the wall (f,,),_) is found to be wa vValy y"(0) (table). The normal 
stress at the wall (p,,),_)= @4@?¥7?/2, which is independent of both viscosity 
and cross-viscosity. The stream function is y = 7*,f/(z). The stream lines are drawn 
for K = 0, 0-1, 0-2, 0-3 (Figure 2). For the same value of y stream lines are drawn 
closer to the wall as well as to the axis with the increase of K. 


a 
0-00 1-8684 1-2699 0-16 1-5630 1-2523 
0-02 1-8086 1-2670 0-18 1-6250 1-2208 
0-04 1-7567 1-2649 0-20 1-7141 1-1836 
0-06 1-7107 1-2635 0-22 1-8061 1-1506 
0-08 1-6688 1-2628 0-24 1-8978 1-1216 
0-10 1-6307 1-2625 0-26 1-9887 1:0961 
0-12 1-5963 1-2625 0-28 2.0790 1-0735 
0-14 1:5752 1-2587 0-30 2:1688 1-0534 

2:25 1:30 
2:00 1:20 
1:75 110 

oe iS 

= ise 
te 0-10 0-20 0:30 
K 


Figure 1 
The graph of 6 and p”(0). 
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Figure 2 
The graph of w. 


4. Conclusion 


The effect of, cross-viscosity in the flow of non-Newtonian liquid with a 
stagnation point depends upon a non-dimensional parameter K = yp, a/mu. The 
boundary layer thickness decreases when K varies, from K = 0 to K = 0-15 and 
then increases more rapidly (Figure 1). The shearing stress at the wall decreases 
slowly upto K = 0-15 and then rapidly upto K = 0-30 (beyond which calculations 
have not been made, Figure 1). It is interesting to see (Figure 1) that 6 has the 
same value when K = 0 and K ~ 0:22 but the flow characteristics in the two 
cases are different. In the latter case, the difference arises due to the cross-vis- 
cosity. The normal stress at the wall is independent of both viscosity and cross- 
viscosity. The stream function is given by y = rv? f(z). We find that for the same 
value of y the stream lines are drawn closer to the wall as well as to the axis with 
the increases of Kk. 

I thank Dr. S. D. Nrcam for his help in the preparation of this paper. 


Zusammenfassung 


Die Bewegungsgleichungen fiir die Str6mung einer nicht-Newtonschen Fliissig- 
keit nahe einem Staupunkt sind integriert worden nach der Karman-Pohlhausen- 
Methode. Es ist gefunden worden, dass der Effekt der Querviskositat von einem 
dimensionslosen Parameter K = ,a/ju abhangt. Die Grenzschichtdicke nimmt 
ab von K=0 zu K= 0,15 und wachst dann mit K. Die Wand-Scherbean- 
spruchung nimmt ab mit zunehmendem K, wahrend die Normalkraft an der Wand 
unabhangig ist von Viskositat und Querviskositat. 


(Received: August 12, 1957.) 
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Herbsttagung der Schweizerischen Physikalischen Gesellschaft 
vom 22. September 1957 in Neuchatel 


Berichte tiber die Sitzungen fiir angewandte Physik und Mathematik 


Ein magnetisch abstimmbares Filter im Mikrowellengebiet, unter Ver- 
wendung der paramagnetischen Resonanz. Von E. BaLDINGER und B. STapD- - 
LER, Basel?), 

Um die Méglichkeit eines magnetisch abstimmbaren Filters im Mikrowellen- 
gebiet zu untersuchen, haben wir auf der Grundlage der paramagnetischen Reso- 
‘nanz eine Briickenanordnung mit gekreuzten rechteckigen Hohlleitern aufgebaut. 
Im Gegensatz zu den bekannten Methoden zum Nachweis der paramagnetischen 
Resonanzabsorption?) verwenden wir in unserem Experiment den Induktions- 
effekt in Anlehnung an das Blochsche Prinzip der Kerninduktion. 
Figur 1 zeigt zwei gekreuzte, rechteckige Wellenleiter, deren Seiten und 
Achsen paarweise aufeinander senkrecht stehen. Diese Anordnung besitzt die 
Eigenschaft einer Briickenschaltung. Bei vollkommener Symmetrie dieser Anord- 
nung sind die beiden Hohlleiter elektrisch iiber ein breites Frequenzband ent- 
koppelt. Um dies zu zeigen, betrachten wir die elektrischen und magnetischen 
Felder eines in seinem Grundzustand angeregten, rechteckigen Hohlleiters (Hy - 
Welle). Der elektrische Feldvektor liegt im Hohlleiter 4 parallel zur y-Richtung, 


Paramagnetische 
Substanz 
H, i 
h | fingang 
° 1 
me) 2 O'G \e! 
e Oo e' ! 
fels © fOJOJOO ele! z 
| fe} [e) Oo @!,! 
Keven eat. Ue Bled 
Sf cn Die eas = 2 
TSSSSRSSSSSASATSS SSS 
{uso Hy-Welle 
Figur 1 Figur 2 
Prinzip des magnetisch abstimmbaren Filters. H9-Welle im rechteckigen Hohlleiter. 
Feldlinien: ————> e; ————> fA. 


1) Physikalisches Institut der Universitat, Abteilung fiir angewandte Physik. 
2) Siehe zum Beispiel G. Feuer, Bell Syst. techn. J. 36, 449 (1957). 
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I 
die zugehérigen magnetischen Feldlinien verlaufen in der (¥%, z)-Ebene. Setzt man 
den in seiner (x, y)-Ebene geschnittenen Hohlleiter A symmetrisch auf die Breit- 
seite eines zweiten Hohlleiters B gemass Figur 1, so tritt dort durch die Offnung 
ein hochfrequentes Magnetfeld, welches aus Symmetriegriinden keine gemein- | 
samen Komponenten mit einem H,,)-Zustand im Hohlleiter B besitzt. Daher wird 
im Hohlleiter B keine Welle im Grundzustand angeregt, sondern nur hdhere 
Schwingungszustande, fiir welche durch die Wahl der Abmessungen die Rand- 
bedingungen nicht erfullt sind. 

Bringt man in die Verzweigung eine paramagnetische Substanz, so erzeugt 
ein ausseres statisches Magnetfeld H, = H, einen selektiven Durchlassbereich in 
unserer Briicke. Der Durchlassbereich liegt bei der Frequenz 

ee oe 
wobei y das gyromagnetische Verhaltnis bedeutet. 

Dieser selektive Durchlassbereich kommt wie folgt zustande: Vom Wellen- 
leiter A her wirkt in der Verzweigung die h,-Komponente auf die paramagnetische 
Substanz. Diese h,-Komponente lasst sich in eine links- und eine rechtszirkulare 
Komponente aufspalten. Ein in der Verzweigung befindliches magnetisches 
Moment nimmt dann im Resonanzfall Energie von jener Komponente auf, 
welche mit dem Drehsinn seiner Larmorprazession im dusseren statischen Magnet- 
feld H, iibereinstimmt. Auf diese Weise wird im Hohlleiter B die Symmetrie 
gestort und damit ein magnetisches Wechselfeld in der (*, y)-Ebene angeregt, 
welches zu einer existenzfahigen H,)-Welle im Hohlleiter B Anlass gibt. 

Mit einer derartigen Briicke ist es méglich, ein magnetisch abstimmbares 

Filter aufzubauen. Als Testsubstanz haben wir ein organisches freies Radikal, 
Diphenyl-Pikryl-Hydrazyl (DPPH), in Form eines kristallinen Pulvers verwendet. 
y hat den Wert 2,8 MHz/Oe, die Linienbreite (und somit im wesentlichen auch der 
Durchlassbereich des Filters) betragt etwa 5 MHz. Bei einer vorlaufigen Messung 
ergab sich im Resonanzfall ein Ubertragungsfaktor von —30 db (Pulver-Fiull- 
faktor rund 0,3) im Bereich 8-10 kMHz. Der Abgleich der Briicke ausserhalb der 
Resonanz war etwa —50 db. Der verwendete Magnet hatte einen Poldurchmesser 
von 25 cm und einen Polabstand von 4,5 cm. Die Feldinhomogenitat betrug im 
Bereich der Briicke weniger als 0,5 Oe/cm. 

Als Unterlage fiir die theoretische Abschatzung der Briickeneigenschaften 

haben wir folgende Daten an DPPH gemessen: 
IDichves) bein 02(C amr sq as et. hn o = 1,40 + 0,03 g/cm3, 
Dk. von DPPH-Pulver, Fiillfaktor 0,4. é= 2,01 + 2%, 
NSTIUSEWITEL Me tee e Agee See iret tgd = 0,02 + 10%. 
Da uns zurzeit keine Kristalle geniigender Griésse zur Verfiigung standen, 


haben wir die an Pulver gemessenen Gréssen auf den festen Korper umgerechnet. 
Wir erhalten nach BérrcHEr?) 


Bey = 47035 ted 


X 5) 


4 0,05. 


fes 


Bei F, und M. Brutn, Physica 22, 130, wird die Dichte von DPPH mit 1,2 + 108 kg/m® 
angegeben, ein Wert, der uns zu klein erscheint. Unser Wert ist ein Mittelwert aus 3 verschiedenen, 
sorgfaltig durchgefiihrten Messungen. Bei 20°C ergab:; 


1. Schwebemethode, mit gut entgastem Pulver. . ........ - @ = 1,42 + 0,05 g/cm3, 
2. Pyknometermethode, Kristall von rund 200 mg unter Wasser entgast 
(durch Anlegen eines leichten WWM HTIS SE BY GS Roa MSG oe Pe @ = 1,37 + 0,05 g/cm%, 


-. . @= 1,40 + 0,05 g/em’. 
C. J. F. Borrcuer, Theory of Electric Polarisation (Elsevier Publishing Company, Amsterdam 
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Die Dk.-Messungen wurden mit Kristallpulver nach der Methode von A. v. Hip- 
PEL’) mittels kurzgeschlossenem Hohlleiter ausgefiihrt®). 

Fiir die Unterstiitzung unserer Arbeiten sind wir der Schweizerischen Kom- 
mission fiir Arbeitsbeschaffung zu grossem Dank verpflichtet. Herr Professor Dr. 
P. HuBER hat uns freundlicherweise den Magnet zur Verfiigung gestellt. Ferner 
moéchten wir an dieser Stelle den Herren B. Kratzer und H. D. KiRNER (Uni- 
versitat Basel) fiir die Darstellung des freien Radikals DPPH bestens danken. 


Utilisation des transistors dans les dispositifs de résonance nucléaire. 
Par P. Cornaz et R. MERcIER, Lausanne?). 


Introduction 


Pour étudier les possibilités qu’offrent actuellement les transistors dans la 
construction d’appareils demandant un faible bruit propre et une bonne stabilité, 
nous avons réalisé un dispositif complet de résonance nucléaire n’utilisant que des 
transistors. Il couvre la gamme de fréquence de 0,5 a 1,5 MHz. Actuellement 
V’observation se fait sur l’oscillographe, mais nous ajouterons prochainement un 
détecteur synchrone («lock-in») qui permettra de passer sur un enregistreur. 

Le principe du montage est la mesure de l’absorption ou de la dispersion 
nucléaire a l’aide d’un pont en T shunté (figure 1) décrit par GRIVET, SOUTIF et 
Buys?) (méthode de PURCELL). 


200 
-9V 


ny) é 
o-——= récepteur 


EB echantillon 


Figure 1 
Oscillateur et Pont. 


Oscillateur 


Un dosage soigné des couplages M, et M, (figure 1) permet d’obtenir une 
stabilité de fréquence trés grande. La variation relative de la fréquence due a une 
variation de la température du transistor a pour valeur A//f}, = 3+ 10-5/°C. Celle 
due & une variation de la tension d’alimentation est Af/f) = 4-10~°/V (alimen- 


5) A. v. HippE., Dielectric Materials and Applications (Wiley and Sons, New York 1954). 

6) Es sei noch darauf hingewiesen, dass eine ahnliche Anordnung zum Bau eines Spektro- 
graphen fiir paramagnetische Resonanzuntersuchungen verwendet werden kann. Zu diesem Zweck 
muss die Verzweigung in Figur 1 zu einem Resonator umgestaltet werden, in welchem zwel 
orthogonale Schwingungszustande gleicher Frequenz méglich sind (zum Beispiel: Zylindrischer 
TMy41- oder TM,19-Schwingtopf). Derartige Untersuchungen sind im Gange. 

1) Laboratoire de Physique de l’EPUL. , 

2) P, Grivet, M. SouriF et M. Buyze, C. R. Acad. Sci., Paris 229, 113 (1949). 
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tation normale 9 V). Les variations de la charge ont une influence plus importante 
sur la fréquence; mais, la puissance de sortie utilisée étant faible, nous avons pu 
limiter cet effet a A//f,) = 3 - 10-4 entre sortie court-circuitée et sortie ouverte, par | 
l’adjonction d’une charge fixe en série avec la charge utile. Pour des puissances de 
sortie supérieures a 0,5 mW, il est cependant nécessaire d’ajouter un étage de 
sortie neutrodyné et faiblement couplé a l’oscillateur. 


Récepteur 


C’est un récepteur hétérodyne composé d’un étage changeur de fréquence, un 
étage moyenne fréquence (80 kHz), un étage démodulateur et deux étages basse 
fréquence (figure 2). Le tout prend place dans une boite de 25 x 17 x 17 em: 

Le facteur de bruit de cet amplificateur est de 3,5 environ (en puissance). Ce 
chiffre pourra sans doute étre amélioré, en sélectionnant le premier transistor et 
en remplacant le changeur de fréquence par un étage HF. 


-9V 
27k 68k 
15k [5000 
56k $5k Op 
=e LPS 
Os 
16) 0C71 
E 1 1 
; *H00p 47k 100} 
o0 
Figure 2 
Récepteur. 
Conclusion 


Les paramétres qui décrivent le fonctionnement des transistors sont complexes 
et dépendent de la fréquence. Ceci rend le calcul des schémas long et laborieux. 
La dépendance des impédances d’entrée et de sortie de la fréquence, d’une part, 
et des impédances de charge et de source, d’autre part, est souvent trés génante. 
A ces inconvénients nous pouvons opposer les avantages suivants: 

Facilité d’alimentation sur batteries tout a fait courantes et absence d’échauf- 
fement. Ces deux avantages entrainent les suivants: une stabilité exceptionnelle 
des performances et l’absence totale des perturbations de 50 Hz dues au réseau. 
Dans les oscillateurs pilotes, les transistors nous semblent préférables aux lampes 
pour des fréquences inférieures A 5 MHz (50 MHz pour drift-transistor), A condition 
de bien découpler le transistor du circuit oscillant. 


Eine neue Untersetzerschaltung mit Transistoren. Von E. BALDINGER und 
P. SANTSCHI, Basel), 

In Figur 1 und 2 ist eine neuartige Untersetzerschaltung dargestellt, die mit 
vier Transistoren bestiickt ist. Die Arbeitspunkte der Transistoren sind so ge- 
wahlt, dass der dem sperrenden oberen Transistor zugeordnete untere stark iiber- 


1) Physikalisches Institut der Universitat, Abteilung fiir angewandte Physik. 
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_steuert ist, der dem oberen leitenden zugeordnete dagegen nur schwach oder gar 
‘nicht. Die oberen Transistoren (im folgenden als « Haupttransistoren » bezeichnet) 
haben die gleiche Funktion wie die in den tiblichen Schaltungen mit zwei stabilen 
Zustanden?) *) und sind nicht tibersteuert. Die zwei unteren Transistoren (« Hilfs- 

_ transistoren ») werden als « Gedachtnis» der Schaltung verwendet, das heisst, sie 
iibernehmen die Funktion der Kondensatoren C, in der Schaltung von BaLpIN- 
GER und NicoLer?) bzw. der Induktivitat L, in der Schaltung von CHAPLIN und 
OwEns’'). 


Figur 1 Figur 2 
Kapazitaten in pF. 


Der am Eingang eines Untersetzers angelegte Impuls hat die Aufgabe, die 
Schaltung von einem in den anderen stabilen Zustand iiberzufiihren. Es ist des- 
halb wiinschenswert, dass der Eingangsimpuls auf den vorher leitenden Haupt- 
transistor in Sperrichtung und gleichzeitig auf den vorher sperrenden Haupt- 
transistor in Leitrichtung wirkt. In der vorliegenden Schaltung haben die beiden 
Hilfstransistoren die Aufgabe, den Eingangsimpuls, welcher in Sperrichtung auf 
ihre beiden Basisanschliisse wirkt, in der gewiinschten Weise zu transformieren. 
Da der Hilfstransistor der leitenden Seite des Untersetzers nicht iibersteuert ist, 
wird sein Kollektorstrom und damit der Strom durch die leitende Seite des Unter- 
setzers unmittelbar gesperrt. Der Hilfstransistor der gesperrten Seite dagegen ist 
iibersteuert, das heisst, sein Emitterstrom ist gleich seinem Basisstrom. Die Uber- 
steuerung hat zur Folge, dass in der Basis dieses Transistors zahlreiche Ladungs- 
trager gespeichert werden‘). Wird nun an der Basis des iibersteuerten Transistors 
ein Impuls in Sperrichtung angelegt, so bringt man gleichzeitig seinen zugeord- 
neten oberen Transistor zum Leiten (Figur 3). 


2) E. BALDINGER und M. NICOLE, Untersetzer mit Flachentransistoren, ZAMP 6, 503-506 (1955). 

3) G. B. B. Cuapiin und A. R. Owens, A Junction-Transistor-Scaling Circuit with 2 ps 
Resolution, Proc. Inst. elect. Engrs. [B] 103, 510-515 (1956). ; 

4) J. L. Mott, Large-Signal Transient Response of Junction Transistor, Proc. IRE 42, 1773-1784 


(1954). 
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fo 
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Eingang 


Figur 3 
Serieschaltung eines tibersteuerten und eines gesperrten Transistors. 
Die Basisspannung des Transistors T, ist so gewahlt, dass er im Ruhezustand sperrt. Die Basis- 
spannung des T, dagegen so, dass ein betradchtlicher Basisstrom fliesst. Ein positiver Eingangs- 
impuls, der also den Transistor T, sperrt, bringt T, kurzzeitig zum Leiten. 


Die in Figur 1 dargestellte Schaltung zeichnet sich durch ein gutes Auflésungs- 
vermodgen aus. Mit dem Typ SB-100 wurde bei einem Spannungssprung am Kol- 
lektor von 0,5 V®) ein Auflésungsvermégen von 50 mus erreicht. 

In diesem Zusammenhang sei noch auf die Schaltung von CHaPLiIn und 
OweEns*) hingewiesen, die mit Niederfrequenztransistoren etwa das gleiche Auf- 
lo6sungsvermogen aufweist wie unsere Variante 1. Diese Schaltung’) hat jedoch den 
Nachteil, dass sie fiir Hochfrequenztransistoren des Typs SB-100 ungeeignet ist. 

Gegeniiber der Schaltung nach Figur 1 hat die in Figur 2 dargestellte Variante 
zusatzlich je einen Widerstand zwischen Kollektor des oberen und Basis des unte- 
ren Transistors. Sie zeichnet sich durch gute Toleranzen beziiglich der Variation 
der Schaltelemente aus. 

Das gleiche Prinzip der Kombination von iibersteuerten und nicht tibersteuer- 
ten Transistoren wurde mit Erfolg auch bei Ringschaltungen angewendet. 

Fir die finanzielle Unterstiitzung dieser Untersuchungen sind wir der Hasler- 


Stiftung und der Schweizerischen Kommission fiir Arbeitsbeschaffung zu Dank 
verpflichtet. 


Couples de pivotement d’origine hydraulique et centrifuge agissant sur 
les pales de roues Kaplan. Par B. CHarx, Genéve!), 

Les couples de pivotement agissant sur une aube de turbine hydraulique 
Kaplan ont la méme origine que le moment de torsion auquel est soumise une pale 
d’hélice d’avion. A part les efforts résultant de la répartition de pression a la 
surface de l’aube, les forces centrifuges jouent un réle — particuliérement bien 
connu en aéronautique. 

Le couple d'origine centrifuge unitaire, tel qu’il agirait sur l’aube d’un modéle 
de turbine ayant pour diamétre D, = 1 m et pour chute H, = 1 m, peut s’écrire: 


M > i L \2 ieee ae o 
a iene. Gling. 4 ( ie ) doi Sey apie [e a 
13 Wie 


a H, DE Di 
Dieser Spannungssprung reicht nicht aus, um die nachfolgende Stufe direkt anzukoppeln. 
Ausserdem sollte diese hochauflésende Stufe am Ausgang méglichst wenig belastet werden. Zur 
Kopplung der nachfolgenden Untersetzerstufe wurde deshalb vorlaufig ein zweistufiger R6hrenver- 
starker mit anschliessender Kathodenstufe verwendet. 
1) Ateliers des Charmilles S.A. 


By 
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- Dans le terme de droite: 
le 1e? facteur, avec le poids G de l’aube, dépend surtout de la chute; 
le 2° facteur contient le coefficient de chute ou de pression usuel dans les turbo- 
machines; 
‘le 3¢ facteur, avec la valeur relative de la corde moyenne L/D,, dépend avant tout 
du nombre des aubes; 
le 4°, qui rapporte les moments d’inertie volumétriques principaux de la pale a 
son volume V et a L?, varie remarquablement peu d’une roue a l'autre, gé- 
néralement entre 0,05 et 0,06; 
le 5¢ est fonction de langle de calage «, repéré par la direction principale de l’aube. 
Ce couple de pivotement tend toujours a la fermeture. A l?emballement, il 
peut atteindre 20 kg/m*. Cette valeur est alors, multipliée par le nombre des 
aubes, du méme ordre de grandeur que le couple moteur unitaire de la turbine. 
Pour une Kaplan tripale, le couple de pivotement d’origine hydraulique a 
. été déterminé au Laboratoire de mécanique des fluides des Ateliers des Charmilles. 
L’angle de calage «, variant entre 0 et 30°; la vitesse de rotation passant de zéro 
a la vitesse d’emballement, et au dela; l'appareil directeur étant fixe. Les résul- 
‘tats de mesure ont été présentés d’une maniére apparentée au coefficient de couple 
aérodynamique des profils isolés: 


ou S est la surface de l’aube et C,, la vitesse débitante moyenne. 

Dans tout le domaine de mesures, ce coefficient restait dans les imites com- 
prises entre — 0,1 et + 0,3. La tendance prédominante a la fermeture a été 
obtenue en déplacant l’axe du tourillon jusqu’aux 50% de la corde des profils. 

A vitesse normale, le rapport entre les maxima respectifs des couples centri- 
fuges et hydrauliques est de 0,05 a 0,3 suivant l’installation. 

Une publication intégrale aura lieu aprés que le travail ait été complété. 


Ein Gittergerat fir Untersuchungen des Raman-Effektes. Von E. LUSCHER 
und F. Troyon, Lausanne’). 

In der Raman-Spektroskopie werden an die Spektrographen besonders hohe 
Anforderungen an Lichtstarke gestellt, da die Raman-Strahlung sehr schwach 
ist. Die Wahl eines Prismengerates drangt sich in vielen Fallen auf, da man damit 
ein sehr giinstiges Offnungsverhaltnis erzielen kann. Andererseits hat man an den 
vielen Grenzflachen Glas—Luft unerwiinschte Schwachung des Lichtes durch 
diffuse Streuung. Da Untersuchungen sowohl mit der Resonanzlinie Hg 2536 A 
als auch bis hinauf zu Hg 6908 A auf unserem Programm sind und auch Struktur- 
untersuchungen an Raman-Linien angestellt werden sollen, entschlossen wir uns, 
ein Gittergerat zu bauen. Damit die Streuverluste méglichst klein gehalten wer- 
den kénnen, wurde eine Montierung nach PASCHEN-RUNGE?) gewahlt. 

Der Durchmesser des Rowland-Kreises betragt 2m. Das Gitter (B + L) 
besitzt 600 Furchen pro Millimeter, und die gravierte Flache betragt 65 x 104 mm, 
was eine hohe Lichtstarke gewahrleistet. Das theoretische Aufl6severmogen ist 
62400 in erster Ordnung (praktisch erreichen wir 30000, da nicht die ganze 
gravierte Flache ausgeleuchtet wird). Zur Aufnahme der Spektren kann sowohl 


1) Laboratoire de Physique de VEPUL. Unterstiitzt durch den Schweizerischen Nationalfond 


fiir wissenschaftliche Forschung, Nr. 623. é 
2) C. K. Runce und F, PAscHEN, Abh. kénigl. Akad. Wiss. Berlin, Anhang I (1902). 
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eine Filmkamera als auch ein Photomultiplier verwendet werden, def hinter — 


einem Austrittsspalt montiert ist. Spalt mit Photomultiplier k6nnen mittels eines 
Antriebes langs des Rowland-Kreises das Spektrum kontinuierlich abtasten. 


Dieser Wagen wird vom Zentrum des Rowland-Kreises aus radial gefiihrt. Der — 


Gitterhalter wurde so konstruiert, dass er in 6 Freiheitsgraden genau justiert 
werden kann. 

Als Lichtquellen zur Anregung der Raman-Spektren beniitzen wir den innern 
Zylinder einer Westinghouse Mitteldruck-Hg-Lampe A-H12 fiir qualitative Unter- 
suchungen und fiir Feinstruktur-Messungen eine Toronto*)-Spirallampe aus 
Quarz,. hergestellt nach unseren Angaben von der Quarzlampengesellschaft in 
Hanau. Die Elektroden werden mit Wasser gekiihlt, und die Betriebsdaten dieses 
Brenners sind: 100 V, 15 A. 

Die Linien dieser Toronto-Lampe sind dusserst fein. Figur 1 zeigt schematisch 
die Anordnung des Brenners und des Spektrographen. 

Zur Herabsetzung des Streulichtes wurden eine Reihe von Blenden, die mit 
schwarzem Stoff tiiberzogen wurden, eingebaut. 

Unser Dank geht an Herrn Prof. MERCIER fiir viele wertvolle Diskussionen 
und an Herrn RIEBeEN fiir seine prazise Arbeit bei der Konstruktion. 


Figur 1 


if Torontobrenner ; RK Raman-Rohr; P Parabolspiegel; F, Eintrittsspalt; G Gitter; A drehbarer 
Arm fiir Kamera und Photomultiplier; F, Austrittsspalt; Ph Photomultiplier; V7 Vorverstarker. 


The Physics of Fluids 
Published by the American Institute of Physics 


Starting with January 1958, the American Institute of Physics will publish 
a new journal, The Physics of Fluids. 

The new journal will be devoted to original contributions to the physics of 
fluids covering kinetic theory, statistical mechanics, structure and general physics 
of gases, liquids and other fluids; as well as certain basic aspects of physics of 
fluids bordering geophysics, astrophysics, biophysics and other fields of science. 


3) J. W. Kemp, J. opt. Soc, Amer. 41, 868 (1951). 
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The scope of these fields of physics includes: magneto-fluid dynamics, ionized 
fluid and plasma physics, shock and detonation wave phenomena, hypersonic 
physics, rarefied gases and upper atmosphere phenomena, physical aeronomy, 
transport phenomena, hydrodynamics, dynamics of compressible fluids, bound- 
vary layer and turbulence phenomena, liquid state physics and superfluidity. 

The Board of Editors of the new journal consists of: J. W. Beams (University 
of Virginia), W. BLEAKNEyY (Princeton University), J. M. BurGers (University 
of Maryland), S. CHANDRASEKHAR (Yerkes Observatory), H. L. DRYDEN (Natio- 
nal Advisory Committee for Aeronautics), W. M. ELsSAssErR (Scripps Institution 
of Oceanography), R. J. Emricu (Lehigh University), J. C. Evvarp (Lewis Flight 
Propulsion Laboratory), J.O. HIRSCHFELDER (University of Wisconsin), J. KAPLAN 

“(University of California, Los Angeles), J. G. KirKwoop (Yale University), 
G. KuErtI (Case Institute of Technology), A. M. KUETHE (University of Michigan), 
M. S. Piesset (California Institute of Technology), 5S. A. ScHaar (University of 
California, Berkeley), G. B.ScHUBAUER (National Bureau of Standards), 
E. TELLER (University of California, Berkeley), P. P. WEGENER (Jet Propulsion 
Laboratory), and F. N. FRENKIEL, editor (Applied Physics Laboratory, The 
Johns Hopkins University). 

The new journal will consider appropriate for publication original papers 
containing sufficiently significant results not published in any other recognized 
publication. The Physics of Fluids will start as a bimonthly and will become a 
monthly as soon as it appears desirable. 

Correspondence on editorial matters should be addressed to: F. N. FRENKIEL, 
editor, The Physics of Fluids, Applied Physics Laboratory, The Johns Hopkins 
University, Silver Spring, Maryland. 

Subscription information can be obtained from the American Institute of 
Physics, 335 East 45th Street, New York 17, New York. 

American Institute of Physics 


Internationale Tagung fiir Elektronenmikroskopie 1958 in Berlin 


Der 4. Internationale Kongress fiir Elektronenmikroskopie 1958 findet in 
Berlin vom 10. bis 17. September statt. Das Organisationskomitee leitet Prof. 
Dr. Ernst Ruska, Prisident der International Federation of Electron Microscope 
Societies, Berlin-Dahlem, Faradayweg 4-6. Auskiinfte durch den geschafts- 
fiihrenden Vorsitzenden Dr. H. KEHLER, Farbwerke Hoechst AG, Frankfurt am 


Main-Hoechst. H. KEHLER 
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Spheroidal Wave Functions — Including Tables of Separation Constants 
and Coefficients. Von J. A. StRaTTON, P. M. Morse, IDeA (Ciehone AIR D, C. LITTLE 
und F. J. Corsato (The Technology Press of M. I. T. and John Wiley & Sons, 
New York 1956). 613 S.; $12.50. . 

Tafeln der Entwicklungskoeffizienten der Spharoidfunktionen (das heisst der 
Lésungen der in gestreckt- und abgeplattet-rotationselliptischen Koordinaten 
separierten Schwingungsgleichung Ag + kq@= 0 nach zugeordneten Kugel- 
funktionen und Bessel-Funktionen mit halbzahligem Index). Die Tabellen laufen 
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nach zwei Eigenwertindizes m und / sowie nach einem vom Koordinatensystem 
und von k abhangigen Parameter h bzw. g, und zwar in den Bereichen: 


m=0(1)8, 1=m(1)8, h bzw. g = 0(0,1) 1 (0,2) 8. 


Eine ausfiihrliche mathematische Einleitung ist dem Tabellenwerk beigegeben. — 
P. LAUCHLI 


Anleitung zum praktischen Gebrauch der Laplace-Transformation. 
Von Gustav DoetscH (R. Oldenbourg, Miinchen 1956). 198 S., 12 Abb.; Hln. 
DM 22.-. 

“Das Buch richtet sich vor allem an die Praktiker, welche die Laplace-Trans- 
formation als mathematisches Hilfsmittel benutzen. Alle Operationsregeln und © 
Satze sind mathematisch korrekt und unter Angabe aller der bekanntlich oft 
miihsamen Voraussetzungen formuliert; es werden jedoch keine Beweise darge- 
stellt. Die sehr in die Einzelheiten gehende Schilderung der Lésung gewohnlicher 
und partieller Differentialgleichungen macht das Buch zu einem sehr wertvollen 
Hilfsmittel, und auch speziellere Anwendungsgebiete (Differenzengleichungen, 
Integralgleichungen) werden leicht verstandlich und leicht anwendbar auseinan- 
dergesetzt. Besondere Erwahnung verdient der kurze, aber sehr schéne Ab- 
schnitt itiber Asymptotik. Wie in der klassischen Laplace-Theorie iiberhaupt, 
bleibt auch in diesem Buch die «Impulsfunktion» (Dirac-Funktion) ein Fremd- 
korper. 

Das Buch macht die Ergebnisse, die der Verfasser in seinem grossen mathe- 
matischen Werk tiber Laplace-Transformationen entwickelt hat, einem grésseren 
Kreis zuganglich und kann warmstens empfohlen werden. E. STIEFEL J 


Bibliography of Russian Mathematics Books. Von GrorGcE E. ForsyTHE 
(Chelsea, New York 1956). 104 S.; $3.95. é 

Diese verdienstvolle Schrift enthalt die bibliographischen Daten von iiber 
600 Bichern tiber reine und angewandte Mathematik in russischer Sprache und 
die Buchtitel in englischer Ubersetzung. Wie der Autor in der Einleitung richtig 
bemerkt, geht aus der russischen mathematischen Literatur hervor, dass Sowjet- 
Mathematiker die westlichen mathematischen Entwicklungen bedeutend genauer 
verfolgen als wir die russischen, Diese Unsymmetrie kénnte dazu fiihren, dass 
man in Russland tiberhaupt mehr iiber Mathematik weiss als bei uns. Die durch 


den Verfasser gebotene Erleichterung in der Beniitzung russischer Biicher ist 
daher sehr willkommen. E. STIEFEL 


Proceedings of Symposia in Applied Mathematics. Vol. VI; Numerical 
Analysis. Herausgegeben von John H. Curtiss (McGraw-Hill, New York 1956). 
303 S., 22 Abb.; s. 73.-/$9.75. 

Das sechste Symposium der amerikanischen mathematischen Gesellschaft 
iiber ein Thema der angewandten Mathematik fand im August 1953 statt und 
war der numerischen Analysis gewidmet. Der vorliegende Sitzungsbericht gibt 
einen guten Querschnitt durch die damals aktuellen Probleme. Die Vortrage des 
Symposiums betrafen hauptsachlich folgende Gebiete: Methoden sukzessiver 
Approximation zur Lésung linearer Gleichungen und Differentialgleichungen, 
Anwendung der Funktionalanalysis auf Approximationsprobleme, lineare Pla- 
nung und kombinatorische Aufgaben, Anwendung der numerischen Analysis auf 
die Zahlentheorie, konforme Abbildung. Der Bericht enthalt die vorztiglich redi- 
gierten Ausarbeitungen von 19 Vortragen. E. STIEFEL 


Dn 
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Funfzig Jahre Relativitatstheorie. Bern, 11. bis 18. Juli 1955. Helvetica 


‘Physica Acta, Supplement IV. Herausgegeben von A. MERCIER und M. KErR- 


VAIRE (Birkhauser Verlag, Basel 1956). 284 S.; Fr./DM 36.—. 

Seit einem halben Jahr liegen die Verhandlungen des Jubilaumskongresses 
vor, der als Geburtstagsfeier fiir die Relativitatstheorie EINsTEINs geplant war, 
der aber unerwartet zu einer Totenfeier fiir den Geehrten wurde. Diese Tatsache 
allein schon macht den vorliegenden Band zu einem Dokument von historischer 
Bedeutung. — In diesem Lichte interessieren besonders die persénlichen Erinne- 
tungen von L. KoLiros, vor allem aber die gut dokumentierte kritische Dar- 


_ stellung der entscheidenden Jahre der Entwicklung des Relativitaétsprinzips von 


Max Born. — Die Vielschichtigkeit der Referate tiber Spezialprobleme sei in ein 
paar Stichworten angedeutet: Experimentelle Ergebnisse tiber die Lichtablen- 
kung und Rotverschiebung und tiber die Expansion des Universums (hier ist das 
Fehlen des Referats von Baabe sehr zu bedauern), Kosmologie, Integration der 


_ Feldgleichungen, Gravitationswellen, Erweiterungen der Theorie, Deduktion der 


Bewegungsgleichungen, Schwierigkeiten der Quantisierung. Als Uberblick und 
Bereicherung der Diskussion bedarf das tiberlegene Schlusswort des Prasidenten 
W. Pauti besonderer Erwahnung. — Im Ganzen kann man den Herausgebern, 
vor allem A. MERCIER, zu diesem hochinteressanten Bande nur gratulieren. 

C. Enz 


> 


Elements of Partial Differential Equations. Von Jan N. SNEDDON 
(McGraw-Hill Book Company, New York, Toronto, endon S054) 6327o., 
50 Fig.; £2.16.6. 

Die vorliegende Einfiihrung in das Gebiet der partiellen Differentialgleichun- 
gen ist aus Vorlesungen hervorgegangen, die der Verfasser an der Universitat 
Glasgow und anderswo gehalten hatte. Dem gesteckten Ziel entsprechend, steht 
die Bestimmung der Lésungen einer vorgelegten Gleichung im Vordergrund, 
wofiir im allgemeinen mehrere (elementare) Lésungsmethoden angegeben werden, 
wiahrend die allgemeine Theorie nur immer so weit entwickelt und dargestellt 
wird, als dies fiir das Verstandnis der Methoden erforderlich ist. Das Buch gliedert 
sich in die folgenden Kapitel: 

I. Ordinary Differential Equations in More Than Two Variables. 11. Partial 
Differential Equations of the First Order. 111. Partial Differential Equations of the 
Second Order. IV. Laplace’s Equation. V. The Wave Equation. V1. The Diffusion 


Equation. 
Das Buch zeichnet sich durch einen sehr klaren und iibersichtlichen Aufbau 


aus. Besonders wertvoll fiir den Anfanger sind die vielen durchgerechneten 


Beispiele sowie die bemerkenswert grosse Zahl von Ubungsaufgaben, darunter 
viele angewandte Aufgaben aus dem Gebiet der Physik, der Technik und der 
Biologie. Als Ganzes gesehen, handelt es sich um eine ausgezeichnete und sehr 
empfehlenswerte Einfiihrung fiir die Studierenden der Mathematik, der Physik 
und der Ingenieurwissenschaften; das Buch bietet aber auch dem Praktiker, der 
gelegentlich mit partiellen Differentialgleichungen zu tun hat, viele niitzliche 
Hinweise. E. RotH-DESMEULES 


The Theory of Games and Linear Programming. Von S. VajDA 
(Methuen & Co., London 1956). 106 S., 14 Fig.; 8s. 6 d. 

Dieses in der Serie Methuen’s Monographs on Physical Subjects herausgege- 
bene Bandchen gibt einen guten Uberblick tiber ein neues Gebiet der angewandten 


4 


7 


Mathematik, das erst mit der Anwendungsméglichkeit von elektronischen Re-_ 
chenmaschinen seine volle Bedeutung erlangt hat. hobs “oe 
Ausgehend von der elementaren Theorie der Spiele, wie sie von J. v. NEv-. 
MANN im Jahre 1928 entwickelt wurde, versteht es der Verfasser, die Problem- ' 
stellung der « Linearen Programmierung» klar herauszuarbeiten. Anhand von | 
_vielen Beispielen aus den Gebieten der Technik, Statistik und Wirtschaft, von’ 
denen das bekannteste wohl das Transportproblem ist, wird der Leser mit Extre- 
_malaufgaben bekanntgemacht, die als Nebenbedingungen Systeme von linearen — 
Ungleichungen besitzen und deren Lésung mit Hilfe der Variationsrechnung. 
nicht mehr moglich ist. 

Ein kurzes Kapitel ist der graphischen Darstellung des Problemkreises ge- 
widmet, bevor — aufbauend auf elementaren algebraischen Kenntnissen — die 
Simplexmethode zur Lésung solcher Extremalprobleme hergeleitet wird. Die Ver- 
bindung mit der Theorie der Spiele, die sich als Spezialfall der Linearen Pro- 
grammierung ergibt, wird durch die Anwendung der Simplexmethode auf soge- 
nannte « Two-person games» hergestellt. 

Viele im Detail durchgerechnete Aufgaben erlauben es dem Leser, sich schnell 
in das Gebiet einzuarbeiten, wahrend die ausfiihrliche Bibliographie das weitere 
Studium erleichtert. TH. GINSBURG 
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Interpolation and Allied Tables. Prepared by H.M. Nautical Office. 
Published by order of The Lords Commissioners of the Admiralty (Her Majesty’s 
Stationery Office, London 1956). 80 S.; 5s. 

Dieses kleine Buch enthalt alles, was ein Rechner fiir die Interpolation in 
numerischen Tabellen braucht; es ist eine Ausgestaltung der Tabellen gleichen 
Titels, die erstmalig 1936 von L. J. Comrie herausgegeben und haufig abgedruckt ~ 
wurden. Das Buch enthalt zunachst sehr bequeme Tafeln fiir lineare und qua- 
dratische Interpolation und dann Tafeln der héheren Besselschen und Everett- 
schen Interpolationskoeffizienten. Weiter findet man eine Formelsammlung iiber 
Differenzenrechnung und numerische Integration. Erwahnt sei noch, dass ausgiebig 
der Comriesche « Riickwurf» der héheren Differenzen auf die niedrigeren benutzt 
wird und dass nun auch von der Tschebyscheffschen Abkiirzung zu Interpola- 
tionszwecken Gebrauch gemacht wird. E. STIEFEL 


Elektrische und magnetische Potentialfelder. Von H. Bucunorz 
(Springer-Verlag, Berlin 1957). 552 S., 202 Abb.; DM 72.-. 

In diesem Werk werden, ausgehend von den physikalischen Grundlagen, 
zunachst die mathematischen Hilfsmittel fiir die Integration der Gleichung von 
LaPLace gebracht. Diese werden dann auf ebene elektrische Felder, auf elektrische 
Felder von Hochfrequenzleitern, auf raumliche elektrische Felder, auf raumliche 
magnetische und hochfrequente Felder sowie auf ebene magnetische Felder ange- 
wandt. In jedem dieser Falle werden mehrere Beispiele durchgerechnet, welche 
dazu geeignet sind, die Anwendungen der allgemeinen Grundsiatze fiir die Bestim- 
mung der Felder zu zeigen. Es ist sehr erfreulich, dass viele Beispiele auch nume- 
risch durchgerechnet werden, wodurch der Leser ein klares Bild der angewandten 
' Verfahren erhalten kann. Historisch ist einzuwenden, dass zwar M. ABRAHAM um 
die Jahrhundertwende hochfrequente Felder mittels elektrostatischer Felder 
angenahert berechnet hat, dass aber Lord KELvin diese Analogie bereits mehr 
als 20 Jahre friiher ausfiihrlich begriindet und ver6ffentlicht hat. 

Das Buch erfordert ein eingehendes Studium, um den vollen Nutzen aus dem 
dargebotenen Stoff zu ziehen. Es ist jedem zu empfehlen, der sich mit der Bestim- 
mung von Potentialfeldern beschaftigen muss oder will. M. Strutr 
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Neuer scheinung 


VOM DENKEN IN BEGRIFFEN 


Mathematik als Experiment des reinen Denkens 
von 


ALEXANDER ISRAEL. WITTENBERG 


Professor an der Laval-Universitat in Quebec 


mit einem Geleitwort von PAUL BERNAYS, 
Professor an der Eidgenéssischen Technischen Hochschule in Ziirich 


, Sammlung «Wissenschaft und Kultur» — Band 12 
386 Seiten. In Ganzleinen gebunden Fr. 25.—- (DM 25.-) 


’ Nirgends wird die Macht und die innere Kohidrenz begrifflichen Denkens 
so schonungslos auf die Probe gestellt wie in der Mathematik. Im mathe- 
matischen Grundlagenproblem wird, wie dieses Werk in eingehender Dis- 
kussion zeigt, eine echte Problematik menschlichen Denkens offenbar. Die 
Untersuchung derselben fiihrt zu wesentlichen Einsichten in die Natur und 
die Grenzen menschlicher Erkenntnis, und damit tiber den Menschen 

schlechthin. 


Das Buch will neues Licht auf die mathematische Grundlagenforschung 

werfen. Zugleich wendet es sich an alle, insbesondere Wissenschaftler, 

Philosophen und Theologen, die sich iiber die erkenntnismassigen Gegeben- 
heiten des menschlichen Daseins Rechenschaft ablegen wollen. 
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